Search

We present a method to jointly invert surface wave dispersion data and gravity measurements for 3-D shear wave velocity and density models. We implemented a petrophysical approach to combine the kernels of both methodologies in a single process. The synthetic experiments show that jointly inverted models recover shear wave velocity and density better than separate inversions. In particular, density models benefit from the good vertical resolution of surface wave dispersion data, while shear velocity models benefit from the good lateral resolution of gravity data. We also proposed two methods to stabilize the solution when using high-grade polynomials. We applied the methodology to the Los Humeros Geothermal area to demonstrate its applicability in a complex geological scenario. Compared with separate inversion, the joint inversion contributes to enhancing key aspects of the geothermal system by (i) delimitating better the geometry of the caldera deposits in the first 0–2.8 km deep by increasing the vertical resolution in density, (ii) delimitating better the lateral borders of low- Vs bodies at different depths interpreted as a part of a complex magmatic chamber system and (iii) estimating the local shear wave velocity–density relationship that conforms to other known relationships for sedimentary and igneous rocks but with some differences that bring us additional information. [ABSTRACT FROM AUTHOR]

The combination of electrical resistivity and seismic refraction tomography is a common practice for the characterization of subsurface features. Presently, the cross-gradient inversion scheme stands out as one of the most robust joint approaches, and some authors modified it to manage complex topographies on unstructured meshes even if at the expense of introducing additional parameters in the inversion process. We propose in this work a cross-gradient algorithm for jointly inverting electrical and seismic tomographic data on structured meshes in cases with non-flat topography. The proposed approach preserves the benefit of the classical cross-gradient approach without the need to impose physical length scales, as for irregular meshes. The quality of the results is evaluated in comparison with independent inversion through a new standardized cross-gradient index and a fuzzy c -means analysis that provides an assessment of the reconstruction accuracy through the membership function. The proposed method was applied to both synthetic models and field-scale examples located in Central Italy, where an accurate geophysical reconstruction is needed for the rehabilitation of existing dams. For all cases, joint inversion yielded superior results compared to independent inversion, demonstrating better agreement with available borehole data. The effectiveness of the joint approach was also demonstrated by the post-inversion tools, where the new cross-gradient index highlighted changes in structural similarity whilst fuzzy c -means clustering allowed for a quantitative reconstruction (position and shape) of the main units at the sites, facilitating the detection of site layering modifications. [ABSTRACT FROM AUTHOR]

3-D shear wave velocity and density models are important for understanding the structures, material composition, tectonic deformations and dynamical mechanisms of the Earth's crust. Such models are usually based on surface wave tomography and gravity inversion, which provide high resolution in the vertical and horizontal directions, respectively. The joint inversion of surface wave and gravity methods can promote the mutual constraints and complementary advantages of seismic and gravity information, improving the imaging resolution and reducing the uncertainty in the individual methods. However, the traditional joint inversion methods directly construct only the shear wave velocity model, excluding density models. We present a joint inversion method of surface wave and gravity data that simultaneously constructs both crustal shear wave velocity and density models. Unlike the previous studies, we reconfigure the seismic kernels of surface wave tomography to preserve the seismic kernel of density. Moreover, the gravity kernel of density is combined with the seismic kernels to establish an objective function of simultaneous joint inversion. Consequently, the imaging resolution of density structure is improved. Our method is validated on the northeastern Tibetan plateau. The inversion results show that the Dingxi, Jiuzhaigou and Jishishan earthquakes occurred in the high–low-anomaly transition zones of shear wave velocity or density, implying that they were induced by accumulation of strain energy in the upper crust of the northeastern Tibetan plateau when the surrounding harder blocks extruded during the tectonic deformation process. [ABSTRACT FROM AUTHOR]

The Pannonian Basin, situated in Central Europe, is surrounded by the Alpine, Carpathian and Dinaric orogens. To understand its tectonic characteristics and evolution, we determine a shear wave velocity model of its crust, mantle lithosphere and asthenosphere consistently by jointly inverting Rayleigh wave phase velocities measured consistently from earthquake (EQ) and ambient noise (AN) data. For the AN data, continuous waveform data were collected from 1254 stations, covering an area within 9° from the centre of the Pannonian Basin during the time period from 2006 to 2018. This data set enabled the extraction of over 164 464 interstation Rayleigh phase-velocity curves, after applying a strict quality control workflow. For the EQ data set more than 2000 seismic events and about 1350 seismic stations were used in the broader Central and Eastern European region between the time-span of 1990 to 2015, allowing us to extract 139 987 quality controlled Rayleigh wave phase-velocity curve. Using the combined data set, a small period- and distance-dependent bias between ambient noise and earthquake measurements, mostly below 1 per cent but becoming larger towards longer periods has been found. After applying a period and distance dependent correction, we generated phase-velocity maps, spanning periods from 5 to 250 s. 33 981 local dispersion curves were extracted and a new approach is introduced to link their period-dependent roughness to the standard deviation. Using a non-linear stochastic particle swarm optimization, a consistent 3-D shear wave velocity model (PanREA2023) encompassing the crust and upper mantle down to 300 km depth was obtained with a lateral resolution reaching about 50 km at the centre of the study area for shorter periods. The crust beneath the Carpathian orogen exhibits a distinct low-velocity anomaly extending down to the Moho. It is referred to as Peri-Carpathian anomaly. Similar anomalies were observed in the Northern Apennines, while the Eastern Alps and Dinarides, as collisional orogens, generally demonstrate higher velocities in the upper crust. High crustal shear wave velocities are also evident in the Bohemian Massif and the East European Craton. The brittle upper crust of the Pannonian Basin is characterized by alternating NE–SW trending high- and low-velocity anomalies: the western and central Pannonian low-velocity anomalies and the Transdanubian and Apuseni high-velocity anomalies related to Miocene sedimentary basins and intervening intervening interbasinal highs exposing Pre-Cenozoic rocks including crystalline basement rocks. Beneath the Southeastern Carpathians, a NE-dipping slab was identified, extending to depths of at least 200 km, while a slab gap is evident beneath the Western Carpathians. A short south-dipping Eurasian slab was imaged beneath the Eastern Alps down to only 150–200 km depth. The Adriatic lithosphere is subducting near-vertically dipping beneath the Northern Apennines, and a slab gap was observed beneath the Central Apennines. In the Northern Dinarides, a short slab was evident, reaching depths of around 150 km. The Southern Dinarides featured a thinned but possibly incompletely detached slab. [ABSTRACT FROM AUTHOR]

Direct geological information in Antarctica is limited to ice free regions along the coast, high mountain ranges, or isolated nunataks. Therefore, indirect methods are required to reveal subglacial geology and heterogeneities in crustal properties, which are critical steps toward interpreting geological history. We present a 3D crustal model of density and susceptibility distribution in the Wilkes Subglacial Basin (WSB) and the Transantarctic Mountains (TAM) based on joint inversion of airborne gravity and magnetic data. The applied "variation of information" technique enforces a coupling between inferred susceptibility and density, relating these quantities to the same gravity and magnetic sources to give an enhanced inversion result. Our model reveals a large body located in the interior of the WSB interpreted as a batholithic intrusive structure, as well as a linear dense body at the margin of the Terre Adélie Craton. Density and susceptibility relationships are used to inform the interpretation of petrophysical properties and the reconstruction of the origin of those crustal bodies. The petrophysical relationship indicates that the postulated batholitic intrusion is granitic, but independent from the Granite Harbor Igneous Complex described previously in the TAM area. Emplacement of a large volume of intrusive granites can potentially elevate local geothermal heat flow significantly. Finally, we present a new conceptual tectonic model based on the inversion results, which includes development of a passive continental margin with seaward dipping basalt horizons and magmatic underplating followed by two distinct intrusive events associated with the protracted Ross Orogen. Plain Language Summary: Most rocks in Antarctica are hidden beneath a thick ice sheet. Therefore, indirect techniques are required to reveal rock provinces within Earth's crust below the ice. Rocks simultaneously influence the gravity and magnetic fields through their physical properties (density and susceptibility). Here we use both the gravity and magnetic fields to reveal rock provinces beneath the ice and use the relationship between density and susceptibility of the rocks to interpret the distribution of granitic rocks of the Transantarctic Mountains and the Wilkes Subglacial Basin region in East Antarctica. Granitic rocks can lead to elevated heat flow due to radiogenic decay of minerals within the rock and thus influence the overlying ice sheet. Based on our subsurface model of rock provinces we speculate on the tectonic evolution of the region. Key Points: We present a new 3D crustal density and susceptibility distribution model based on joint inversion of gravity and magnetic dataDensity and susceptibility data are used to identify crustal level intrusions and the craton marginOur new conceptual tectonic model identifies a earlier intrusive event in the interior of the Wilkes Subglacial Basin, separate from those exposed in the Transantarctic Mountains [ABSTRACT FROM AUTHOR]

In the field of mineral resource exploration, accurate imaging of subsurface structures is key to discovering and assessing potential mineral deposits. Traditional single geophysical methods, limited by terrain variations and their own constraints, can lead to divergent solutions and structural inconsistencies, affecting the reliability of exploration outcomes. To address these challenges, this paper presents a joint inversion method for three-dimensional direct current (DC) resistivity and gravity data based on a deformed hexahedral mesh. The article begins by outlining the current state of development of the method under study and proposes a research plan, followed by a detailed explanation of the theoretical basis and algorithmic implementation of the proposed method. Model tests confirm the advantages of the deformed hexahedral mesh in reducing terrain impacts and enhancing model resolution, demonstrating the optimization and complementarity of the resolution between the two methods after joint inversion. Finally, applying this method to actual data from the Huaniu Mountain area shows that joint inversion not only improves the consistency of the ore belt structure but also provides a more precise analysis for the quantitative interpretation of the distribution of underground mineral resources. This confirms the method's effectiveness and potential in practical geological exploration. [ABSTRACT FROM AUTHOR]

We present a general framework for multiphysics joint inversion of any number of geophysical data sets. Its main feature is the use of the variable splitting approach: an auxiliary multiparameter model space is introduced in which minimization of the coupling and stabilizing functionals is carried out. The use of rediscretization and interpolation to map between this auxiliary space and the model spaces allows the coupled models to have completely different parametrizations. Joint inversion is decoupled into the individual inversion and the coupling-regularization subproblems, each of which can be solved by a different optimization algorithm. For each subproblem, the linking term controlling the distance between the model and the corresponding auxiliary variable takes the form of a quadratic regularization with a reference model. As a result, any existing inversion code supporting such regularization can be integrated without modifications into the developed framework. As a concrete example scheme, we consider an application of the framework to 3-D joint inversion of magnetotelluric, seismic refraction and gravity data. We discuss different coupling functionals, mainly those corresponding to the more universal structural constraints: joint total variation, joint minimum support, cross-gradient, one-way cross-gradient and their combinations for a general multimodel case. The use of coupling based on explicit 'petrophysical' relationship between the properties is also considered. Performance of the developed framework is studied on three synthetic cases: a time-lapse joint inversion of full-tensor gravity gradiometry and seismic data, a joint inversion of magnetotelluric, seismic and gravity data and a joint inversion for electrical resistivity tomography and audio-magnetotellurics. [ABSTRACT FROM AUTHOR]

The 2008 M w 7.9 Wenchuan earthquake ruptured the middle and northeastern segments of the Longmenshan Fault Zone (LMSFZ), and the 2013 M w 6.6 Lushan earthquake ruptured a 50-km-long fault in the southwestern segment. Subsequently, an M w 5.8 earthquake occurred approximately 10 km distant from the M w 6.6 Lushan earthquake. Therefore, the potential risk for larger earthquakes (> M w 6.6) on the southwestern section must be considered. This study collects the latest seismological and GPS data to construct an integrated seismotectonic model for the two neighbouring earthquake sequences. The model integrates the fault planes involved, the main shock rupture processes, the main shock-caused Coulomb stress perturbation, the aftershock distribution and the 3-D velocity structure of the source region, providing information for seismic risk evaluation. We find that three fault planes were involved, two of which were related to the main shocks, and the third was generated by the aftershocks following the first main shock. The main shocks were caused by nearly pure thrust faulting on the two planes with dip angles of approximately 45° and almost opposite dipping directions, thereby forming a conjugate angle of around 90°. The third plane was located between the two main shocks, approximately parallel to the second main shock's fault plane. Each of the main shocks primarily ruptured a single asperity, displaying simple time history. The Coulomb stress change of the first main shock facilitated the generation of the second main shock and the third fault plane, and the second main shock increased the stress on the first main shock's fault plane. The aftershocks were distributed within stratified materials by spatially varying interfaces and characterized by high Vp and Vs velocity, and a low Vp / Vs ratio. The atypical dip angles of approximately 45° for thrust faults and the conjugate angle of approximately 90° are indicative of high stress state. The single asperity rupture implies simple stress accumulation. The main shock-caused Coulomb stress change did not reduce the seismic risk in the source region. The varying interfaces are interpreted as a consequence of long-term horizontal compression. All of these characteristics suggest that the two earthquake sequences were generated by the breakage of three immature faults under strong compression by background stress, and the high stress state remains within the southwestern LMSFZ. [ABSTRACT FROM AUTHOR]

Bedrock geology from Antarctica remains largely unknown since it is hidden beneath thick ice sheets. Geophysical methods such as gravity and magnetic inverse modelling provide a framework to infer crustal rock properties indirectly in Antarctica. However, due to limited availability of rock samples, validation against direct geological information is challenging. We present a new rock property catalogue containing density and susceptibility measurements on 320 rock samples from northern Victoria Land. This catalogue is used to assess the reliability of local and regional scale inverse results, including a new local high resolution magnetic inversion in the Mesa Range region and a previously published regional scale joint inversion of gravity and magnetic data in northern Victoria Land and the Wilkes Subglacial Basin. We compare our density and susceptibility measurements to global and local measurements from the literature to access the correlation to rock types and geological units. Furthermore, the measured values are compared against inverted values. The close correspondence between inverted and measured rock properties allows us to predict locations of rock types where currently such information is missing. The utility of measured susceptibility and density relationships for interpreting inversion output provides a strong incentive to incorporate local rock samples into geophysical studies of subglacial geology across Antarctica. [ABSTRACT FROM AUTHOR]

SUMMARY: We present a novel strategy for performing joint inversion with guided fuzzy c-means (GFCM) clustering coupling and apply it to electrical resistivity tomography (ERT) and ambient noise surface wave (ANSW) data. To accurately extract a priori clustering information, we use density peak clustering (DPC) rather than fuzzy c-means (FCM). The number and centres of resistivity and shear-wave velocity a priori clusters are extracted by DPC and then used to guide the joint inversion with the GFCM clustering coupling of ERT and ANSW data. Synthetic and field data are used to evaluate the flow and algorithm of DPC-GFCM clustering joint inversion. The results of synthetic examples show that the models recovered by the DPC-GFCM clustering joint inversion are nearly the same as the true models and are more accurate than those inverted using individual inversion and FCM-GFCM clustering joint inversion. In the field case, the depths of the stratigraphic interfaces shown in the resistivity and shear-wave velocity models inverted by DPC-GFCM clustering joint inversion are nearly consistent with those from the drilling data. In contrast, the strata recovered by the individual inversion and FCM-GFCM clustering joint inversion significantly differ from the drilling results. Both the synthetic and field examples verify the effectiveness of the DPC-GFCM clustering coupling method used for the joint inversion of ERT and ANSW data acquired from the near surface with strong heterogeneity. This novel approach can also be applied to other types of geophysical data. [ABSTRACT FROM AUTHOR]

Subsurface characterization is a crucial aspect of geophysical exploration, enabling the identification and understanding of valuable geological bodies and resources. In this context, joint inversion of gravity and magnetic data has emerged as a powerful geophysical exploration technique, allowing for a more coherent and consistent interpretation of subsurface structures. The study focuses on understanding residual gravity and magnetic anomalies by employing the gradient descent-based joint inversion approach. A MATLAB program was developed to determine the inverse gravitational and magnetic anomalies using the gradient descent approach. We explored the potential of 2D rectangular prisms as a popular geometry to represent mineralized bodies and oil and gas structures. To overcome the non-uniqueness issues, we designed code for joint inversion of gravity and magnetic data. Synthetic data was inverted using the gradient descent technique and compared with the least-squares approach. Numerical simulations and real data application successfully reconstructed the geometry of the prisms. An illustrative example of a prism fault was used for further evaluation. Real data from the Oka complex in Quebec, Canada, was collected from the literature and subjected to joint and individual gravity and magnetic modelling. The results highlighted the influence of heterogeneous mass distribution on matching forward anomalies. The high gravity anomaly in the Oka complex was attributed to carbonatite and silicate rocks. The presence of two intrusive centres within the complex caused the magnetic high. This work demonstrates the effectiveness of the gradient descent approach as it consistently outperformed the conventional method, offering a robust solution for subsurface characterization in geophysical exploration. [ABSTRACT FROM AUTHOR]

There is great exploration potential for ultra-deep reservoirs in Junggar Basin， and reservoir prediction with seismic data is the main technical means of oil and gas exploration at present. However， due to the low signal-to-noise ratio of seismic data of ultra-deep reservoirs， unclear reservoir-seismic data correspondence， and few actual drilling wells， it is difficult to establish an effective initial model for seismic inversion， which restricts the accuracy and reliability of ultra-deep reservoir prediction. Gravity inversion， as an important quantitative interpretation method， can obtain the characteristics of underground density distribution and provide support for geological interpretation. According to the density model， a relatively reliable low-frequency model can be established for seismic inversion， which can overcome the difficulty of applying seismic data in ultra-deep reservoirs to a certain extent. Meanwhile， the acquisition of gravity data is economical and convenient compared with that of seismic data， and it is easier to be applied in practice. Therefore， a new technique applying gravity inversion to seismic data-based reservoir prediction was developed in this article. Firstly， a quasi-neural network gravity inversion technique based on Gaussian radial basis functions was proposed to solve the gravity inversion problem and improve the resolution and reliability of gravity inversion. Then， the density body obtained by gravity inversion was used as training data， and a neural network was trained together with seismic and logging data to establish the initial model of seismic inversion. Finally， seismic inversion was carried out with the initial model constraints. This technique broke through the application bottleneck of single seismic data in ultra-deep reservoir prediction and overcame the limitations of logging constraints， providing a reliable initial model for seismic inversion. The application of this technique to the prediction of ultra-deep clastic rock reservoirs in Junggar Basin was consistent with the existing geological knowledge， indicating that this method had high practical value and application potential for ultra-deep reservoir prediction and could provide technical support for ultra-deep reservoir exploration.

The precise relocation of seismic events is critical for many engineering projects. Swarms of minor or micro earthquakes typically reveal stress concentration and spots of greater seismic hazards. Particularly in the context of deep underground mining, advanced techniques that can accurately relocate microseismicities are urgently in demand. Here, we developed a neural network-based modeling training method that can precisely relocate seismicities and invert for velocities at the same time, with preconfigured receiver network locations. Our model can be iteratively improved with field recorded data. We showed that, with roughly eight iterations, we can reasonably resolve for the earthquake locations for both clusters of events, namely spatially distributed with linear pattern or randomly scattered. Our initially trained model, which only focused on events that had a linear distribution pattern, was used as the base for the training of the subsequent models which could better resolve for randomly scattered event locations. Although we stopped at the eighth iteration, the process reported here can be continued, as the model will have a better performance with more iterations. [ABSTRACT FROM AUTHOR]

SUMMARY: For accurate modelling of groundwater flow and transport processes within an aquifer, precise knowledge about hydraulic conductivity K and its small-scale heterogeneities is fundamental. Methods based on pumping tests, such as hydraulic tomography (HT), allow for retrieving reliable K-estimates, but are limited in their ability to image structural features with high resolution, since the data from time-consuming hydraulic tests are commonly sparse. In contrast, geophysical methods like induced polarization (IP) can potentially yield structural images of much higher resolution, but depend on empirical petrophysical laws that may introduce significant uncertainties to the K-estimation. Therefore, this paper presents a joint inversion procedure for both HT and IP data, which allows for combining the complementary abilities of both methods. Within this approach, a traveltime inversion is applied to the HT data, while the IP inversion is based on a full-decay time-domain forward response, as well as a reparametrization of the Cole–Cole model to invert for K directly. The joint inversion is tested on a synthetic model mimicking horizontally layered sediments, and the results are compared with the individual HT and IP inversions. It is shown that jointly inverting both data sets consistently improves the results by combining the complementary sensitivities of the two methods, and that the inversion is more robust against changes in the experimental setups. Furthermore, we illustrate how a joint inversion approach can correct biases within the petrophysical laws by including reliable K-information from hydraulic tests and still preserving the high-resolution structural information from IP. The different inversion results are compared based on the structural similarity index (SSIM), which underlines the robustness of the joint inversion compared to using the data individually. Hence, the combined application of HT and IP within field surveys and a subsequent joint inversion of both data sets may improve our understanding of hydraulically relevant subsurface structures, and thus the reliability of groundwater modelling results. [ABSTRACT FROM AUTHOR]

SUMMARY: Reconciling rock unit boundary geometry is crucial for geological and geophysical studies aiming to achieve a comprehensive 3-D subsurface model. To create a unified 3-D parametrization suitable for both geological modelling and geophysical inversion, an integrated approach utilizing implicit modelling is essential. However, a key challenge lies in encapsulating all pertinent information within the 3-D model, ensuring compatibility with the utilized data sets and existing constraints. In this study, we present a workflow that enables the generation of an integrated 3-D subsurface model primarily using gravity and reflection seismic data sets. Our approach involves a cooperative geophysical inversion workflow, which incorporates the inverted model from the reflection seismic data while leveraging sparse petrophysical information. Despite advances in integrated modelling, the incorporation of implicit modelling approaches in cooperative inversion workflows remains unexplored. In our gravity inversion process, we use a generalized level set method to refine the boundaries of rock units in the prior model. We integrate the inverted model, derived from seismic and other sparse petrophysical data sets, to create a comprehensive 3-D prior model. To enhance the integration of reflection seismic data sets in the level set inversion, we introduce a weighting uncertainty matrix containing constraint terms. This step refines the model's accuracy and ensures greater consistency. Finally, we search for any missing rock units within inverted model through nucleation investigations. The introduced methodology has undergone successful testing in the Boulia region (Southern Mount Isa, Queensland), utilizing two 2-D reflection seismic profiles and regional gravity data sets. This study primarily aims to reconstruct the geometry of major structures within the basement units and the basin at a regional scale. By combining seismic profiles and gravity data sets with constraining information, we are able to create a 3-D model of the area that accurately represents distinct rock units and their boundary geometries. Additionally, relevant legacy data sets and prior modelling results from the region have been incorporated and refined, ensuring that the final model aligns with all available knowledge about the area. [ABSTRACT FROM AUTHOR]

SUMMARY: The misalignment of the observation and predicted waveforms in regional moment tensor inversion is mainly due to seismic models' incomplete representation of the Earth's heterogeneities. Current moment tensor inversion techniques, allowing station-specific time-shifts to account for the model error, are computationally expensive. Here, we propose a gradient-based method to jointly invert moment-tensor parameters, centroid depth and unknown station-specific time-shifts utilizing the modern functionalities in deep learning frameworks. A $L_2^2$ misfit function between predicted synthetic and time-shifted observed seismograms is defined in the spectral domain, which is differentiable to all unknowns. The inverse problem is solved by minimizing the misfit function with a gradient descent algorithm. The method's feasibility, robustness and scalability are demonstrated using synthetic experiments and real earthquake data in the Long Valley Caldera, California. This work presents an example of fresh opportunities to apply advanced computational infrastructures developed in deep learning to geophysical problems. [ABSTRACT FROM AUTHOR]

SUMMARY: Greenland's tectonic history is complex, and the resulting lithospheric structure is, although extensively researched, not well constrained. In this study, we model the lithospheric structure of Greenland in a consistent, integrated framework with three steps. First, we build a lithospheric background model by forward modelling, adjusted to gravity gradient data and shear wave velocities from a regional tomography model. Subsequently, we jointly invert for the upper crustal density and susceptibility structure by minimizing the gravity residuals and magnetic total field anomaly misfit. The last modelling step searches for upper crustal thermal parameters to fit our model to the most recent geothermal heat flow predictions for Greenland. Finally, we present 3-D models of the density, temperature and velocity structure for the lithosphere as well as thermal parameters and susceptibilities for the upper crust. Our model also includes the depth of the Moho and LAB in Greenland. A comparison between inverted crustal parameters and surface geology shows a clear correlation. The novelty of our model is that all these results are consistent with each other and simultaneously explain a wide range of observed data. [ABSTRACT FROM AUTHOR]

Hydraulic tomography (HT) is a promising technique for high-resolution imaging of subsurface heterogeneity, which addresses the limitations of traditional methods, such as borehole drilling and geophysical surveys. This study focuses on the application of HT to detect and characterize fractured zones in bedrock and addresses the gap in the understanding of the role of distributed flux data in the joint inversion of hydraulic head and flux data. By conducting synthetic injection tests and using sequential successive linear estimators for inversion, the study explores the effectiveness of combining limited head data with distributed temperature sensing (A-DTS)-derived flux data. The findings highlight the fact that integrating flux data significantly enhances the accuracy of identifying fracture permeability characteristics, even when head data is sparse. This approach not only improves the resolution of hydraulic conductivity fields but also offers a cost-effective strategy for practical field applications. The results underscore the potential of HT to enhance our understanding of groundwater flow and contaminant transport in fractured media, which has important implications for carbon capture, enhanced geothermal systems, and radioactive waste disposal. [ABSTRACT FROM AUTHOR]

The Tanlu Fault Zone (TLFZ) and Chifeng‐Kaiyuan Fault (CKF) serve as tectonic boundaries between the North China Craton (NCC) and the Central Asian Orogenic Belt (CAOB). Clarifying the refined structure of these tectonic boundaries is crucial for understanding the relationships between the tectonic units and the heterogeneity in the destruction of the NCC. In this study, two linear seismic arrays were deployed across these tectonic boundaries. Based on the phase velocity dispersion and receiver functions extracted from the seismic arrays, the Hamiltonian Monte Carlo algorithm was employed for the joint inversion of the S‐wave velocity (Vs) in the crust and uppermost mantle. The Vs model was then used to correct the time differences in common conversion point (CCP) stacking. The CCP stacking results indicate that the boundary faults TLFZ and CKF are both whole‐crustal faults that separate the NCC and CAOB. The Vs structure showed a significant low‐velocity anomaly in the mantle beneath the NCC, with intense seismic activity within the crust. This suggests that the NCC was affected by the subduction of the Western Pacific, leading to crustal and mantle destruction. In contrast, the CAOB exhibited a clear high‐velocity anomaly with relatively stable crustal structures. We believe that the NCC and CAOB have undergone structural modification and destruction due to the closure of the Paleo‐Asian Ocean and the activities of the TLFZ since the Late Mesozoic. During the Cenozoic, the region east of the TLFZ experienced more significant destruction in the NCC than the other adjacent tectonic units. Plain Language Summary: The joint inversion of Rayleigh wave phase velocities and receiver functions provides a reliable method for obtaining high‐precision S‐wave velocity (Vs) models. Here, we applied a joint inversion based on the Hamiltonian Monte Carlo algorithm to acquire detailed crust‐mantle Vs models in Northeast China and utilized this velocity structure to constrain the common conversion point (CCP) stacking to obtain a more reliable Moho depth. A comprehensive high‐resolution seismological model was established by integrating the Vs model, Moho depth, Vp/Vs ratio, and seismic activity. We then revealed the differential structural features and seismic heterogeneity on both sides of the tectonic boundaries between the Central Asian Orogenic Belt (CAOB) and the North China Craton (NCC). Our research results offer new constraints on the crust and upper mantle structure in Northeast China, contributing to the understanding of the geodynamic context underlying the differential structures between the CAOB and the NCC and the characteristics of their boundaries. Key Points: We obtained the fine crust‐mantle structure of tectonic boundaries in NE China using Rayleigh wave phase velocities and receiver functionsThe tectonic boundaries Chifeng‐Kaiyuan Fault and Tan‐Lu Fault are both whole‐crustal faults with strong seismic activityThe upper mantle of North China Craton exhibits low S‐wave velocity with strong seismic activity [ABSTRACT FROM AUTHOR]

Accurately determining the seismic structure of the continental deep crust is crucial for understanding its geological evolution and continental dynamics in general. However, traditional tools such as surface waves often face challenges in solving the trade‐offs between elastic parameters and discontinuities. In this work, we present a new approach that combines two established inversion techniques, receiver function H‐κ stacking and joint inversion of surface wave dispersion and receiver function waveforms, within a Bayesian Monte Carlo (MC) framework to address these challenges. Demonstrated by synthetic tests, the new method greatly reduces trade‐offs between critical parameters, such as the deep crustal Vs, Moho depth, and crustal Vp/Vs ratio. This eliminates the need for assumptions regarding crustal Vp/Vs ratios in joint inversion, leading to a more accurate outcome. Furthermore, it improves the precision of the upper mantle velocity structure by reducing its trade‐off with Moho depth. Additional notes on the sources of bias in the results are also included. Application of the new approach to USArray stations in the Northwestern US reveals consistency with previous studies and identifies new features. Notably, we find elevated Vp/Vs ratios in the crystalline crust of regions such as coastal Oregon, suggesting potential mafic composition or fluid presence. Shallower Moho depth in the Basin and Range indicates reduced crustal support to the elevation. The uppermost mantle Vs, averaging 5 km below Moho, aligns well with the Pn‐derived Moho temperature variations, offering the potential of using Vs as an additional constraint to Moho temperature and crustal thermal properties. Plain Language Summary: Knowing the seismic structure of the deep crust helps us understand Earth's geological history and how continents evolve. However, traditional methods of studying the deep crust face challenges due to tradeoffs that can impact accuracies of the results. In this paper, we present a new approach that combines two existing techniques intending to measure the deep crust more accurately. We tested this method using both synthetic and real data and learned that it works better than previous methods. We applied this method to the Northwestern US and found that the results are aligned with the area's geology, suggesting that the new method is feasible to be applied on a regional scale. The new method provides a more accurate way to study the deep crust and improves the mapping of the uppermost mantle. Key Points: We have developed a new joint inversion approach that incorporates stacking of receiver function multiple phases with multiple data setsThe new approach reduces the trade‐offs and improves the determination of deep crustal shear velocity, Moho, and Poisson's ratioApplication of the new method to the northwestern US produces a more accurate model that exhibits geologically coherent structures [ABSTRACT FROM AUTHOR]

Copyright of Petroleum Geology & Recovery Efficiency is the property of Petroleum Geology & Recovery Efficiency and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Abstract—A new technique for simultaneous construction of the "gravity" and "magnetic" images of Mercury from satellite data based on the regional version of the S approximation is proposed. The mathematical statement of the inverse problem of finding the images of a planet from the data on the potential fields recorded at different instants with different degrees of accuracy is reduced to solving ill-conditioned systems of linear algebraic equations (SLAE) with approximate right-hand sides. Based on the analytical approximations of the mercury magnetic and gravity fields determined from the solution of the ill-conditioned SLAE, the distributions of the equivalent sources on the spheres are determined. The results of the mathematical experiment on constructing the magnetic image of Mercury from the radial component of the magnetic induction vector continued analytically to the field sources are presented. [ABSTRACT FROM AUTHOR]

We present a kinematic slip model and a simulation of the ensuing tsunami for the 2020 M w 7.0 Néon Karlovásion (Samos, Eastern Aegean Sea) earthquake, generated from a joint inversion of high-rate GNSS, strong ground motion and InSAR data. From the inversion, we find that the source time function has a total duration of ∼20 s with three peaks at ∼4, 7.5 and 15 s corresponding to the development of three asperities. Most of the slip occurs at the west of the hypocentre and ends at the northwest downdip edge. The peak slip is ∼3.3 m, and the inverted rake angles indicate predominantly normal faulting motion. Compared with previous studies, these slip patterns have essentially similar asperity location, rupture dimension and anticorrelation with aftershocks. Consistent with our study, most published papers show the source duration of ∼20 s with three episodes of increased moment releases. For the ensuing tsunami, the eight available gauge records indicate that the tsunami waves last ∼18–30 hr depending on location, and the response period of tsunami is ∼10–35 min. The initial waves in the observed records and synthetic simulations show good agreement, which indirectly validates the performance of the inverted slip model. However, the synthetic waveforms struggle to generate long-duration tsunami behaviour in simulations. Our tests suggest that the resolution of the bathymetry may be a potential factor affecting the simulated tsunami duration and amplitude. It should be noted that the maximum wave height in the records may occur after the decay of synthetic wave amplitudes. This implies that the inability to model long-duration tsunamis could result in underestimation in future tsunami hazard assessments. [ABSTRACT FROM AUTHOR]

Three-dimensional inversion based on source growing uses systematic searches. Compared with the regularization inversion, this method has lower computational requirements and faster processing speed. The criteria determining the source growth is crucial for the quality of the results. This study proposes an inversion based on source growing with full-tensor gravity gradiometry data to improve vertical inversion effectiveness. First, a depth weighting function is introduced for the criteria to optimize the determination of source growing at different depths. Second, the weights of different data are adjusted based on the inversion results of single-component gradient data, and a joint inversion method is established. Finally, matrix compression reduces memory occupation and improves computational efficiency. By the tests of synthetic data and real data from Vinton Dome, it is demonstrated that the proposed method can effectively guide source growing, providing stronger ability for distinguishing deep targets and being suitable for the inversion of complex-shaped targets. Furthermore, the method has high computational efficiency and anti-noise ability. [ABSTRACT FROM AUTHOR]

The Mw 6.8 Adassil earthquake that occurred in the High Atlas on September 8, 2023, was a catastrophic event that provided a rare opportunity to study the mechanics of deep crustal seismicity. This research aimed to decipher the rupture characteristics of the Adassil earthquake by analyzing teleseismic waveform data in conjunction with interferometric synthetic aperture radar (InSAR) observations from both ascending and descending orbits. Our analysis revealed a reverse fault mechanism with a centroid depth of approximately 28 km, exceeding the typical range for crustal earthquakes. This result suggests the presence of cooler temperatures in the lower crust, which facilitates the accumulation of tectonic stress. The earthquake exhibited a steep reverse mechanism, dipping at 70°, accompanied by minor strike-slip motion. Within the geotectonic framework of the High Atlas, known for its volcanic legacy and resulting thermal irregularities, we investigated the potential contributions of these factors to the initiation of the Adassil earthquake. Deep seismicity within the lower crust, away from plate boundaries, calls for extensive research to elucidate its implications for regional seismic hazard assessment. Our findings highlight the critical importance of studying and preparing for significant seismic events in similar geological settings, which would provide valuable insights into regional seismic hazard assessments and geodynamic paradigms.

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the ChinArray II temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method (Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vS structures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vS ratios. While, lower velocities and higher vP/vS ratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane (SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust. Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones (LVZs) in the SPGZ. The crustal thickness, vS, and vP/vS ratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Because of its complex tectonic setting, the Philippines experiences vigorous seismic activity. On December 15, 2019, an Mw 6.8 earthquake struck Davao del Sur, the Philippines, severely damaging infrastructure and adversely affecting the population. In this study, we used Interferometry Synthetic Aperture Radar (InSAR) and teleseismic broadband data to explore the earthquake’s source parameters, rupture process, and tectonic characteristics through a Bayesian modeling approach and Coulomb stress analysis. An analysis of InSAR suggested that the mainshock occurred on an NW–SE oriented single plane within the Cotabato fault system, with the Makilala-Malungon fault being the causative fault. The Joint inversion revealed that the rupture initiated in the middle crust (depth: approximately 16.5 km), exhibiting a blind, left-lateral strike-slip fault rupture propagation with minimal reverse slip and no surface rupture. The slip was concentrated in a distinct asperity spanning a depth of 3–24 km, with the peak slip reaching approximately 1.8 m. The source time function of this event indicated a single pulse with a total source duration of approximately 17 s. The static stress drop value was low, which may be related to the brittle condition of the crust. Coulomb stress changes calculations suggested earthquake risk in the future.

The phase of secular evolution of continents is manifested as the degree of compositional differentiation, modification and maturation of continental crusts, which is vital in understanding the mechanism of continental evolution but is difficult to quantify. Here we use integrated passive- and active-source seismic profiling to conduct joint analysis and inversion and derive Vs and Vp/Vs section models across the North China Craton (NCC) to southeastern Altaids boundary zone that bears a tectonic transition from a reworked ancient craton margin to a Phanerozoic accretionary orogen. We systematically exploited the imaged multiple physical properties as precise and delicate proxies to constrain the compositional architecture in the crust across this important tectonic transition subject to various crustal evolutional phases. Our Vs and Vp/Vs imaging, together with the existing isotopic data, characterizes the Yin Shan–Yan Shan belt as the northern NCC margin with layered homogeneous compositions that point to an evolved crust. However, the lower-crustal low-Vs/high-Vp/Vs signature that overlaps the shallowly dipping to horizontal reflective fabrics suggests that the crust of the northern NCC margin has undergone considerable reworking through lower-crustal-stretching-assisted melt migration and mixing since the late Paleozoic to Mesozoic eras. The process probably involved crust–mantle interaction and thus resulted in a compositionally modified ancient crustal basement. On the contrary, the southeastern Altaids domain manifests crustal complexity in compositions and structures inferred to be indicative of a juvenile crust of the Phanerozoic accretionary orogen. Our results provide deep physical-property constraints that shed new light on the crustal evolution of a complex craton margin. [ABSTRACT FROM AUTHOR]

More accurate inversion of source fault geometry and slip parameters under the constraint of the Bayesian algorithm has become a research hotspot in the field of geodetic inversion in recent years. In nonlinear inversion, the determination of the weight ratio of the joint inversion of multisource data is more complicated. In this context, this paper proposes a simple and easily generalized weighting method for inversion of source fault parameters by joint geodetic multisource data under the Bayesian framework. This method determines the relative weight ratio of multisource data by root mean square error (RMSE) value and can be extended to other nonlinear search algorithms. To verify the validity of the method in this paper, this paper first sets up four sets of simulated seismic experiment schemes. The inversion results show that the joint inversion weighting method proposed in this paper has a significant decrease in the large residual value compared with the equal weight joint inversion and the single data source joint inversion method. The east–west deformation RMSE is 0.1458 mm, the north–south deformation RMSE is 0.2119 mm and the vertical deformation RMSE is 0.2756 mm. The RMSEs of the three directions are lower than those of other schemes, indicating that the proposed method is suitable for the joint inversion of source parameters under Bayesian algorithm. To further verify the applicability of the proposed method in complex earthquakes, the source parameters of the Maduo earthquake were inverted using the method of this paper. The focal depth of the inversion results in this paper is closer to the focal depth released by the GCMT agency. In terms of strike angle and dip angle, the joint inversion in this paper is also more inclined to the GCMT results. The joint inversion results generally conform to the characteristics of left-lateral strike-slip, which shows the adaptability of this method in complex earthquakes. [ABSTRACT FROM AUTHOR]

In near-surface surveys, shallow-seismic and ground-penetrating radar (GPR) full-waveform inversions (FWIs) have received increasing attention because of their ability to reconstruct high-resolution subsurface models. However, they have different sensitivities to the same targets and thus may yield conflicting geophysical parameter models. To solve this issue, we have developed an indirect joint petrophysical inversion (JPI) integrating shallow-seismic and multi-offset GPR data. These data are used to reconstruct porosity and saturation whereby we use only strong sensitivities between petrophysical and geophysical parameters. To promote its field application, we proposed an input strategy to avoid measuring rock matrix parameters and make indirect JPI more robust. We apply indirect JPI to the field data acquired in Rheinstetten, Germany and find that it reveals the mechanical, electrical and petrophysical properties more reliably than individual inversions. The reconstructed models are assessed by direct-push technology, borehole sample measurements and migrated GPR image. Indirect JPI can fit seismic and GPR observed data simultaneously and provide consistent multiparameter models, which are hard to achieve by FWIs and individual petrophysical inversions. We also find that the method is robust when there are uncertainties in petrophysical a priori information. Overall, the field example proves the great potential of using indirect JPI to solve real-world problems. [ABSTRACT FROM AUTHOR]

The diffusion field theory has been widely used to interpret ambient noise wave fields. Based on this theory, 1-D subsurface velocity structure inversion method is developed. However, few studies have referred to the noise horizontal to vertical (NHV) spectral ratio inversion of 3-D subsurface velocity structures, and almost no effective 3-D NHV inversion tools have been developed. To develop a useful tool for obtaining 3-D soil layer velocity structures, we combined the NHV forward calculation formula derived from diffusion field theory with the guided Monte Carlo algorithm and then extended the single-point NHV inversion to multipoint joint inversion through a joint objective function. Subsequently, a new 3-D soil layer velocity structure inversion method was proposed. Subsequently, a synthetic 2-D case was used to verify the proposed method. Finally, the proposed method was applied to the Xiangtang Array in Tangshan, China, to identify the 3-D velocity structures of the site based on noise observations. The results show that the proposed multipoint joint 3-D inversion method is effective for identifying 3-D underground velocity structures. [ABSTRACT FROM AUTHOR]

Intracrustal low-velocity zones (LVZs) indicate a mechanically weak crust and are widely observed in the southeast margin of the Tibetan Plateau. However, their spatial distributions and formation mechanisms remain controversial. To investigate their distribution and detailed morphology of the LVZs in the southeastern Tibetan Plateau, here we used teleseismic events and continuous waveform data recorded by 40 broadband seismic stations newly deployed in the Sichuan-Yunnan region from December 2018 to October 2020. A total of 12,924 high-quality P-wave receiver functions and 5–40 s fundamental Rayleigh surface wave phase velocity dispersion curves from ambient noise cross-correlation functions were obtained. The S-wave velocity model at a depth interval of 0–100 km in the study area was inverted by using the trans-dimensional Markov chain Monte Carlo strategy to jointly invert the complementary data of the receiver function waveform and Rayleigh surface wave phase velocity dispersion. Our results show that there are two separate LVZs (∼3.5 km/s) surrounding the rigid Daliangshan subblock at crustal depths of approximately 30–40 km, providing new constraints on the geometry of the LVZs in our study region. The two LVZs obtained in this study may represent the middle crustal flow channels, through which the material in the center of the Tibetan Plateau extrudes to its southeast margin. Blocked by the rigid Sichuan Basin and the spindle-like Daliangshan subblock, the material continues to flow southward through the mechanically weak middle crustal channels surrounding the Daliangshan subblock. In addition, the existence of thin LVZs in the middle crust plays an important role in understanding the decoupling between the upper and lower crust in the study area. It also provides new constraint on the complex tectonic deformation process of the southeastern margin of the Tibetan Plateau caused by the collision and compression of the Indian and the Eurasian plates. [ABSTRACT FROM AUTHOR]

In order to explore the constraint degree of GNSS and InSAR data on the Maduo earthquake under the layered earth model structure, and to understand the cause of the Maduo earthquake and the danger of the surrounding area, this paper uses D-InSAR technology to obtain the InSAR co-seismic deformation field of the Maduo earthquake on May 22,2021. Based on GNSS and InSAR data, the co-seismic slip model and fault plane stress distribution of the Maduo earthquake are jointly inverted. This paper calculates the Coulomb stress changes caused by 14 M ≥ 7 strong earthquakes co-seismic, post-seismic viscoelastic relaxation and inter-seismic tectonic stress loading of 19 fault segments in the Bayan Har block research area (96°E-106°E, 29°N-36°N) since 1900. The results show that: (1) The maximum Line Of Sight (LOS) deformation obtained based on the ascending and descending InSAR data was about 0.9 m. The maximum slip amount of GNSS + InSAR joint inversion is 4.87 m, which is located in the Dongcao Along Lake section (∼98.6°E). Combined with the stress drop distribution of the fault surface, it is found that the Maduo earthquake is dominated by strike-slip. (2) The historical strong earthquakes delayed the Maduo earthquake by nearly 103 years. The 1937 Huashixia earthquake, the 1947 Dari earthquake, the 1963 Dulan earthquake, and the 1973 Luhuo earthquake had the greatest stress unloading effect, which was equivalent to delaying the Maduo earthquake by 53 years, 21 years, 11 years, and 17 years, respectively. The cumulative Coulomb stress at the source location of the Maduo earthquake reached 0.0112 MPa. Most area of the fault sections such as the East Kunlun Fault, the Maduo-Gande Fault, the western and eastern sections of the Dari Fault, the western and eastern sections of the Qingshuihe Fault, the Ganzi-Yushu Fault, and the eastern section of the Kunlunshankou-Jiangcuo Fault are in the Coulomb stress loading area, which should be paid attention to. [ABSTRACT FROM AUTHOR]

Copyright of Journal of Geodesy & Geodynamics (1671-5942) is the property of Editorial Board Journal of Geodesy & Geodynamics and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Enhancing the reliability of inversion results has always been a prominent issue in the field of geophysics. In recent years, data-driven inversion methods leveraging deep neural networks (DNNs) have gained prominence for their ability to address non-uniqueness issues and reduce computational costs compared to traditional physically model-driven methods. In this study, we propose a GMNet machine learning method, i.e., a CNN-based inversion method for gravity and magnetic field data. This method relies more on data-driven training, and in the prediction phase after the model is trained, it does not heavily depend on a priori assumptions, unlike traditional methods. By forward modeling gravity and magnetic fields, we obtain a substantial dataset to train the CNN model, enabling the direct mapping from field data to subsurface property distribution. Applying this method to synthetic data and one-field data yields promising inversion results. [ABSTRACT FROM AUTHOR]

Joint inversion is a crucial approach to mitigate the non-uniqueness problem in geophysical inversion. Nevertheless, existing joint inversion methods fall short of meeting the stringent requirements of high-precision exploration, necessitating the development of new techniques. In this paper, we introduce the Structural Similarity Index (SSIM) as a novel structural consistency constraint for the joint inversion of gravity and magnetic data. Compared with the results of cross-gradient inversion, our method demonstrates outstanding performance and stability. SSIM inversion not only introduces a new class of joint inversion with structural constraints but also enhances the consistency of inversion results in the distribution of physical attribute values. The structural constraints of SSIM inversion are more comprehensive and robust, significantly improving the reliability of the inversion. Both synthetic and real data applications demonstrate that the proposed method can effectively handle both synthetic and real data, yielding outstanding results. [ABSTRACT FROM AUTHOR]

Copyright of Fly Ash Comprehensive Utilization is the property of Hebei Fly Ash Comprehensive Utilization Magazine Co., Ltd. and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The 1 January 2024 Noto-Hanto (Noto Peninsula) earthquake (MJMA 7.6) generated strong ground motion, large crustal deformation and tsunamis that caused significant damage in the region. Around Noto Peninsula, both offshore submarine and partially inland active faults have been identified by previous projects: Ministry of Land, Infrastructure, Transport and Tourism (MLIT) and Japan Sea Earthquake and Tsunami Research Project (JSPJ). We inverted the tsunami waveforms recorded on 6 wave gauges and 12 tide gauges around Sea of Japan and the GNSS data recorded at 53 stations in Noto Peninsula to estimate the slip amount and seismic moment on each of active faults. The results show that the 2024 coseismic slips were 3.5 m, 3.2 m, and 3.2 m on subfaults NT4, NT5 and NT6 of the JSPJ model, located on the northern coast of Noto Peninsula and dipping toward southeast. A smaller slip, 1.0 m, estimated on NT8 on the southwestern end of the 2024 rupture, may be attributed to its previous rupture during the 2007 Noto earthquake. The total length of these four faults is ~ 100 km, and the seismic moment is 1.90 × 1020 Nm (Mw = 7.5). Almost no slip was estimated on the northeastern subfaults NT2 and NT3, which dip northwestward, opposite to NT4–NT5–NT6, and western subfault NT8. Aftershocks including the MJMA 6.1 event occurred in the NT2–NT3 region, but they are smaller than the potential magnitude (Mw 7.1) those faults can release in a tsunamigenic earthquake. Similar features are also found for the MLIT model; the 2024 slip was only on F43 along the northern coast of Noto Peninsula, and northeastern F42 did not rupture, leaving potential for future event. [ABSTRACT FROM AUTHOR]

The thermochemical structure of the lithosphere controls melting mechanisms in the mantle, as well as the location of volcanism and ore deposits. Obtaining reliable images of the lithosphere structure, and its complex interactions with the asthenosphere, requires the joint inversion of multiple data sets and their associated uncertainties. In particular, the combination of seismic velocity and electrical conductivity, along with proxies for bulk composition and elusive minor phases, represents a crucial step toward fully understanding large‐scale lithospheric structure and melting processes. We apply a novel probabilistic approach for joint inversions of 3D magnetotelluric and seismic data to image the lithosphere beneath southeast Australia. The results show a highly heterogeneous lithosphere with deep conductivity anomalies that correlate with the location of Cenozoic volcanism. In regions where the conductivities have been at odds with sub‐lithospheric temperatures and seismic velocities, we observe that the joint inversion provides conductivity values consistent with other observations. The results reveal a strong relationship between metasomatized regions in the mantle and (a) boundaries of geological provinces, elucidating the subduction‐accretion process in the region; (b) distribution of leucitite and basaltic magmatism; (c) independent geochemical data, and (d) a series of lithospheric steps which constitute areas prone to generating small‐scale instabilities in the asthenosphere. This scenario suggests that shear‐driven upwelling and edge‐driven convection are the primary mechanisms for melting in eastern Australia, contrary to the conventional notion of mantle plume activity. Our study presents an integrated lithospheric model for southeastern Australia and provides valuable insight into the mechanisms driving surface geological processes. Plain Language Summary: The lithosphere is the Earth's outermost rigid layer and the focus of important geological processes, such as earthquakes, volcanism, and mineralization. The location of these processes often coincides with deep discontinuities in the lithosphere. Imaging the structure of the lithosphere using geophysical techniques is crucial to fully understand the nature of these processes. The most reliable images of the lithosphere are obtained by joint analysis of two or more geophysical data sets. In particular, the combination of magnetotellurics (an electromagnetic technique) and seismic data holds great potential because of their complementary sensitivity to the Earth's properties. Furthermore, our understanding of the variability of the lithospheric structure is enhanced by combining this joint analysis with a probabilistic approach. This is because probabilistic methods provide a large number of models that can explain the data. Given the substantial data coverage in southeast Australia, we use a new probabilistic approach for the joint analysis of magnetotelluric and seismic data to image the lithospheric structure beneath this region. Our results show a complex lithosphere in line with the location of volcanism and the tectonic history of the region. The lithospheric composition derived from the models provides significant insights into melt production in the area. Key Points: We apply a novel approach for joint probabilistic inversions of 3D magnetotelluric and seismic dataWe use the new method to image the lithosphere‐asthenosphere system beneath southeastern of AustraliaThe imaged lithospheric structure provides insights on metasomatism and melt production in the region [ABSTRACT FROM AUTHOR]

The seismological data recorded at our four semi-permanent broadband seismographs installed by the CSIR-NGRI, Hyderabad, Telangana in the Aravalli mobile belt (AMB) and Rajasthan craton in the northwest Indian shield during the period 2014–2016 are used in the study. The 16 selected regional Indian earthquakes of Mw 5.5–7.8 that occurred in the North India region recorded at the four seismic stations are used to compute the surface wave (both Rayleigh and Love waves) group velocity dispersion characteristics and average 1-D regional shear wave velocity (Vs) structure beneath the North India which covers the northern part of India between the whole Himalayan region and Rajasthan craton including AMB. First, we compute Rayleigh waves (at 7–87 s) dispersion curves and Love waves (at 7–82 s) group velocity dispersion curves by utilizing the standard frequency time analysis (FTAN) of 3-component broadband data recorded at these four 3-component broadband seismic stations. After that, using the logarithmic stacking technique, we compute the final average stacked fundamental mode dispersion curves of surface waves (both Love and Rayleigh waves) for our study area. Then, we finally invert the final average dispersion curves of surface waves (both Love and Rayleigh waves) to compute the average one-dimensional regional shear velocity structure beneath North India. Our modelling result in North India suggests a two-layered crustal structure with a crustal thickness of 39 km. The thickness of the first crustal layer (upper crust) is 20 km with a shear wave velocity (Vs) of 3.30 km/s. On the other hand, the second crustal layer (lower crust) marks a 19-km thickness with a Vs of 3.70 km/s. Our modelling result suggests that the Vs of the upper mantle is found to be 4.35 km/s below the study region. Research highlights: Modelled average one-dimensional regional Vs structure below northern part of India. A 2-layered crustal velocity model below North India. Upper crust at 20 km thick with Vs of 3.30 km/s. And, lower crust at 19 km thick with Vs of 3.70 km/s. The Moho depth is 39 km and Vs below the Moho depth is found to be 4.35 km. [ABSTRACT FROM AUTHOR]

The Debrecen area, as part of the Great Hungarian Plain (GHP), is associated with a multi-aquifer system that is overly exploited to fulfill the development plans. This research aims to jointly interpret and model gravity and magnetic data to map the subsurface geology and structures that govern groundwater occurrence. Various potential field techniques, including spectral analysis, anomaly derivatives, analytical signal, and Euler deconvolution were employed to map the distribution and depth of the buried geological structures. The combination of the potential field techniques enabled the construction of a detailed lineament map, providing valuable insights into the distribution of the subsurface structural features. It was indicated that the main structural trend is NW–SE and NE–SW, that coincides with the main structural trends in Hungary. Subsequently, a lineament density map is derived, indicating that the eastern, central, and northwestern parts of the area form the most promising zones for groundwater prospection. The joint inversion of gravity and magnetic data has further enhanced the understanding of subsurface geology. The depth to the basement rock varied between 1.18 and 2.2 km. The highest depth to the basement meets with thick sedimentary sequences bounded by normal faults forming graben and horst structures. Moreover, the distribution of these sediments is investigated using lithological logs indicating the thickness of the main hydrostratigraphic units in the Debrecen area. These units include Nagyalföld Aquifer, Algyő Aquitard, Endrőd Aquitard, and Miocene Badenian Aquifer units, which mainly consist of sand, silt, marl, and gravel. The recent study demonstrated the effectiveness of the joint interpretation in enhancing the knowledge of lithology and geological structures. However, a detailed geophysical survey is recommended to characterize the hydrostratigraphic units in the Debrecen area. [ABSTRACT FROM AUTHOR]

Copyright of Remote Sensing for Natural Resources is the property of Remote Sensing for Natural Resources Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Induced polarization datasets, which assess the overvoltage characteristics of the subsurface, are becoming increasingly prevalent. These datasets are not limited to mineral resource exploration but are also gaining traction in various engineering applications. We present the result of joint inversion of direct current resistivity and induced polarization data. We have obtained two subsurface properties: electrical resistivity and chargeability, by minimizing the single objective function. For this purpose, we have computed the distorted apparent resistivity using the secondary voltage measure at each time gate. We compared our results with the conventional inversion of induced polarization on synthetic data with both high and low chargeable mineralized ore bodies. We noted that our approach provides a more reasonable model for high chargeable ore bodies. However, both joint inversion and nonlinear approaches give similar results for low-chargeable ore bodies. Nonlinear inversion of induced polarization was accomplished by solving two objective functions: One for resistivity and another one for induced polarization. Thus, error in the resistivity inversion propagates in the chargeability computation due to consequences. Further, the developed joint inversion algorithm is also applied to a field data set measured at the Beldih mine located in the West Bengal, India. We compared our final chargeability models with the local geology and geologic models obtained via other geophysical methods on the same profiles and found an excellent correlation in both models. [ABSTRACT FROM AUTHOR]

The leakage of leachate is a crucial but often overlooked hydrological process in the landfill water cycle. The concealed leakage of leachate not only leads to soil and groundwater pollution in surrounding areas but also affects the distribution of water and the metabolism of organic matter within the landfill. To accurately quantify this concealed hydrological process, we propose a method for detecting leachate leakage based on joint inversion of multi-source geophysical exploration data. By integrating multiple geophysical exploration data (resistivity and self-potential information), we reconcile the spatial differences in different exploration data to improve the accuracy of leachate and its pollution plume imaging. Additionally, we introduce an efficient alternating iteration technique within the joint inversion framework to ensure convergence of separately inverted models toward similar spatial structures. Simulation results indicate that: (1) for the early detection of small-scale high-concentration leachate pollution, the proposed confidence-induced joint inversion framework (CI-JIF) improves precision by 15.6 % compared to separate inversions. (2) for the detection of leachate long-term diffusion, CI-JIF accurately delineates the distribution of widespread high-concentration leachate, improving precision by 17.4 % compared to separate inversions. Further on-site experiments demonstrate that CI-JIF can more accurately reconstruct the distribution of leachate.

Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson's ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth's structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M 9 event, or distributed deformation along plate boundaries as constrained from InSAR. [ABSTRACT FROM AUTHOR]

There are abundant igneous gas reservoirs in the South China Sea with significant value of research, and lithology classification, mineral analysis and porosity inversion are important links in reservoir evaluation. However, affected by the diverse lithology, complicated mineral and widespread alteration, conventional logging lithology classification and mineral inversion become considerably difficult. At the same time, owing to the limitation of the wireline log response equation, the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations. To overcome those issues, this study takes the South China Sea as an example, and combines multi-scale data such as micro rock slices, petrophysical experiments, wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks. Specifically, we define the lithology and mineral characteristics through core slices and mineral data, and establish an igneous multi-mineral volumetric model. Then we determine element cutting log correction method based on core element data, and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity. However, it is always difficult to determine the elemental eigenvalues of different minerals in inversion. This paper uses multiple linear regression methods to solve this problem. Finally, an integrated inversion technique for altered igneous formations was developed. The results show that the corrected element cutting log are in good agreement with the core element data, and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction. The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations. [ABSTRACT FROM AUTHOR]

We study the 2016 January 21 (⁠|${{{M}}}_{\rm{w}}$|  5.9) and 2022 January 8 (⁠|${{{M}}}_{\rm{w}}$|  6.7) earthquake sequence that struck the Menyuan region in northwest China's Qinghai province. These two earthquakes are destructive events that occurred around/on the Lenglongling fault (LLLF). Here, we derive the line-of-sight displacement fields of the two earthquakes using Interferometric Synthetic Aperture Radar (InSAR) measurements of Sentinel-1 SAR data, and map the range and horizontal offset fields of the 2022 event using Sentinel-1 amplitude images and Planet-Lab optical images. Based on the offset maps, we determine the detailed surface rupture trace of the 2022 event. We perform slip inversions for the two earthquakes on triangle fault patches whose size increases with depth. Results show that the western branch segment of the 2022 event has a ∼0.5-m normal dip-slip motion. This result contradicts previous inferences on dip-slip sense of this branch segment. We identify a left-stepping fault structure with a ∼5-km step width in the transition zone between the Tuolaishan fault (TLSF) and LLLF, which may serve as a kinematic barrier to prevent further propagation of seismic rupture along the TLSF. Stress calculation shows that a stress drop of ∼0.4 bar produced by the 2016 event on a ∼5-km long LLLF segment may act as a negative stress barrier to suppress rupture propagation of the 2022 event toward the southeast of the LLLF. [ABSTRACT FROM AUTHOR]

Copyright of Acta Geoscientica Sinica is the property of Acta Geoscientica Sinica Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)