The Balmuccia peridotite exposes relatively fresh mantle rocks at the Earth’s surface, and as such it is of interest for geologists and geophysicists. The outcrop is a kilometre-scale feature, yet its extent at depth is insufficiently imaged. Our aim is to provide new constraints on the shape of the density anomaly this body represents, through 3D gravity modelling. In an effort to avoid personal or methodology bias, we hereby launch an invitation and call for participative modelling. We openly provide all the necessary input data: pre-processed gravity data, geological map, in situ rock densities, and digital elevation model. The expected inversion results will be compared and jointly analysed with all participants. This approach should allow us to conclude on the shape of the Balmuccia peridotite body and the associated uncertainty. This crowd effort will contribute to the site surveys preparing a scientific borehole in the area in frame of project DIVE. The full description, the dataset, as well as the tentative timeline can be found at https://zenodo.org/records/10390437

How to cite: Hetényi, G., Baron, L., Scarponi, M., Subedi, S., Michailos, K., Dal, F., Gerle, A., Petri, B., Zwahlen, J., Langone, A., Greenwood, A., Ziberna, L., Pistone, M., Zanetti, A., and Müntener, O.: Participative gravity-modelling challenge to constrain the Balmuccia peridotite body (Ivrea-Verbano Zone, Italy), EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6402, https://doi.org/10.5194/egusphere-egu24-6402, 2024.

The Western Alps constitute a complex and heterogeneous orogenic system, generated by the continental collision between the European and Adriatic tectonic plates. Three main tectonic domains can be identified across local to regional scales, based on geophysical and geological observations: the European and Adriatic domains, and the high-density, high-velocity anomaly known as Ivrea geophysical body (IGB). Despite being one of the best-studied collisional systems in the world, the 3D Western Alpine lithosphere and its along-arc compositional and structural variations are still subjects of investigations.

In this framework, we exploit the inherently-3D information provided by gravity data. In particular, we set up a 3D Bayesian joint inversion of Bouguer gravity anomaly and surface wave dispersion data, to obtain a new 3D ρ-v_S model of the Western Alpine lithosphere. We benefit from the Bouguer anomaly map by Zahorec et al. (2021), obtained by homogeneous processing of gravity data across the Alpine domain, and from seismic data recorded by permanent and temporary seismic networks: e.g. IvreaArray, AlpArray (Hetényi et al. 2017, 2018), CIFALPS I and II (e.g. Paul et al. 2022).

We perform 3D forward gravity modeling by discretizing the study area in unitary volumes of constant density (voxels), accounting for spherical Earth structure and surface topography. The gravity effect of each voxel is pre-computed, and then only needs to be scaled with density during the inversion. This significantly decreases the computational cost of the forward model, and thus allows us to explore the parameter space with Monte Carlo sampling. We use a Bayesian framework and implement a Markov chain Monte Carlo (McMC) algorithm. We test different types of  parameterizations to reduce the non-uniqueness of gravity inversion. We plan to jointly invert gravity with surface wave dispersion data, providing complementary information on v_S. Finally, existing receiver function studies (e.g. Monna et al. 2022, Paul et al. 2022, Michailos et al. 2023) provide prior information on crustal and lithospheric geometry.

We expect to obtain a new 3D ρ-v_S model for the Western Alpine crust and lithosphere. This will provide new information on the European-Adriatic collision boundary, together with the IGB structure, and their three-dimensional variation along the orogen. The new model will be also useful to constrain rock composition, upon comparison with the geological observations at the surface.

How to cite: Scarponi, M. and Bodin, T.: Bayesian Joint Inversion of Bouguer Gravity and Surface Wave data: application to the Western Alps, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17176, https://doi.org/10.5194/egusphere-egu24-17176, 2024.

In Western Europe, the Variscan belt contacts Avalonia along the Rhenohercynian Suture, a result of Early Carboniferous continental collision. Moving east of the Harz Mts., the Rhenohercynian suture disappears beneath a thick sedimentary sequence of the Permian-Mesozoic basin. Its extension is either truncated by major NW-SE strike-slip faults like the Elbe, Odra, or Dolsk faults or bends under the cover of a thick sedimentary succession. The extension of Avalonia into Poland is challenging to determine, with the thinned margin of Baltica considered the substratum of the Permian-Mesozoic basin. Deep seismic soundings show that the thinned margin of Baltica reaches the NW-SE oriented Dolsk or Odra fault, potentially bringing the crust of Baltica into direct contact with the crust of the Variscan internides of the Bohemian Massif. Along the Dolsk fault, there is the two-layered, low-velocity Variscan crust in the SW that contacts the three-layered Baltica crust. The geometry of this contact remains unknown, but the lower, high-velocity crust of Baltica may extend southwest to the Odra fault. In the basement of the sedimentary sequence between the Dolsk and Odra faults, low-grade metamorphosed phyllites with a metamorphic age of approximately 360 Ma are found. They apparently represent a fragment of Variscan metamorphic nappes.

The Variscan front is oriented NE-SW in Western Europe, but in Poland, it bends by 90° to the NW-SE direction, continuing to the border of Ukraine. In southeastern Poland, the front enters the slope of the East European Platform, constituting an undisputed example of a direct contact between the Variscan belt and Baltica. If the geometry of the Variscan front reflects the structure of the orogen, the edge of Baltica must have initially played the role of a transform margin with a right-lateral displacement. NW-SE strike-slip faults, parallel to this margin, truncated the Rhenohercynian and other Variscan sutures from the NE. The following accretion event resulted in NE-SW shortening, either thin-skinned, leading to folding of the external fold-and-thrust belt, or thick-skinned, resulting in the emplacement of the Variscan nappe stack on the Baltica margin.

The last folding of external Variscides in Poland occurred around 305 Ma and was immediately followed by the emplacement of a large igneous province at the Carboniferous to Permian transition. The centre of magmatism was in NE Germany, the area of greatest crustal thinning. The origin of the igneous province was linked to plate boundary forces leading to extension and continental rifting. The latter produced the Mid-Polish trough, an elongated continental rift running NW-SE parallel to the Teisseyre-Tornquist zone. Permian rifting further attenuated the Baltica margin and, jointly with coeval magmatism, reshaped the margin of Baltica masking its contact with the Variscan belt. Toward the east, the continuity of the Variscan internides was disrupted by early Mesozoic rifting in the area of the present-day Carpathians.

How to cite: Mazur, S.: From Carboniferous convergence to Permian continental rifting – the interaction of Baltica with the Variscan belt of Europe at the time of the Pangaea assembly, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4560, https://doi.org/10.5194/egusphere-egu24-4560, 2024.

This interdisciplinary study describes the integration of gravity field analysis, curvature techniques and various spatial applications. The data are based on land-based Free Air and Bouguer gravity data from the AlpArray Gravity Research Group, complemented by recent satellite missions. New seismic and seismological data from the AlpArray initiative and the German MB-4D Priority Program were used as independent boundary conditions for the 3D modeling and inversion of the gravity data. Prior to this modeling, Euler deconvolution, terracing/clustering techniques, and advanced filtering methods were employed to reveal intricate details of the region's gravitational signatures. For example, a distinct zoning of gravity is observed in the central part of the Ligurian Sea, pointing to traces of past rifting processes. Analysis of various curvature parameters (e.g., dip-, min-, max- and shape-curvature) of the processed gravity fields, in particular gradients and residual fields support the identified zonation of the gravity fields, which reflect the geological structures in the crust. The final 3D modeling of the Ligurian Sea area is based on a previous density model of the entire Alpine region and includes density distribution of the upper mantle. These densities were derived from tomographic velocity models, accounting for petrology, temperature, and pressure. Additional information of the upper crust was obtained from the refraction seismic results of the LOBSTER project, offering a comprehensive understanding of spatial phenomena. Calculations of the gravitational potential energy (GPE) provide additional information on local stresses, facilitating a deeper understanding of the flexural rigidity in the area. By elucidating the relationship between processing techniques and 3D modeling, this work advances interdisciplinary interpretation crucial for geological studies in the Ligurian offshore area.

How to cite: Götze, H.-J., Strehlau, R., Anikiev, D., Dannowski, A., and Scheck-Wenderoth, M.: Harnessing Modern 3D Gravity Analysis Techniques: A Study of the Ligurian Offshore Area, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6228, https://doi.org/10.5194/egusphere-egu24-6228, 2024.

Through Stokes kernel function and gravity anomaly, it is possible to calculate the gravity gradient disturbance on the geoid and its external space. For this Stokes’ integral expression, we apply Laguerre wavelet numerical integration to improve the accuracy of its computational results. Meanwhile, compute unified device architecture (CUDA) is used to implement parallel computing on the Graphic Processing Unit (GPU) for speeding up. The full tensors of gravity gradient in the experimental ocean area with 3°×2°are computed. Compared to serial computing, its computing acceleration ratio can be more than 10 times faster. The results of the vertical gravity gradient are compared and validated from the public model from the University of California San Diego.

How to cite: Ye, Z., Bu, J., and Sneeuw, N.: Application of numerical integration and CUDA parallel in the calculation of ocean gravity gradient, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17112, https://doi.org/10.5194/egusphere-egu24-17112, 2024.

The South China Sea, situated at the convergence of the Tethys and Pacific tectonic domains, holds immense geological significance due to its interaction with multiple tectonic plates (Hall 2002; Hayes and Nissen 2005; Metcalfe 2011). With its abundant sedimentary basins, the region is of paramount importance for geological and structural studies, particularly in relation to its potential for oil and gas resources. In this study, we propose the utilization of satellite gravity data to analyze the tectonic structure of the South China Sea, focusing on three key areas:

1. High-resolution construction of gravity gradient anomalies and fault identification: By integrating Fast Fourier Transform algorithms with satellite gravity anomalies and high-resolution terrain elevation data, we obtaina comprehensive dataset of full tensor gravity gradient information. Through spatial analysis of this data, we successfully identify 17 significant and deep faults, as well as partition the study area into 9 distinct tectonic units characterized by well-defined geological structures.

2. Moho Depth Determination and Interpretation: Employing an improved regularization Bott's method, we determine the Moho depth using information obtained from sonar-buoy detection and submarine seismograph detection profiles. Regularization parameters are introduced to ensure the smoothness of the inversion results. By analyzing the distribution characteristics of the Moho and its relationship with tectonic units, we conduct a comprehensive analysis to comprehend the coupling between shallow and deep structures. The resultsreveal distinct regional characteristics in the depth distribution of the Moho surface in the South China Sea, shedding light on the distribution of continental crust, oceanic crust, and the ocean-continent transition zone.

3. Comprehensive Geophysical Analysis: We employ a combination of seismically constrained Moho undulation, gravity data, gravity gradient anomalies, and unconstrained 3D correlation imaging to investigate the crustal structure of the South China Sea. Integrating various geophysical datasets, we gain a deeper understanding of the distribution of continental crust, oceanic crust, and transitional crust within the region. Notably, the results shows that the trench-island arc-back arc basin systemplays a pivotal role in the active continental margin of the Western Pacific. This comprehensive analysis provides valuable insights into the tectonic dynamics and geological processes occurring in the South China Sea region.

*This study was supported by y the Basic Frontier Science Research Program of the Chinese Academy of Sciences (No. ZDBS-LY-DQC028).

Reference:

Hall, R. (2002). Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations, J. Asian Earth Sci. 20:353–431.

Hayes, D.E., Nissen, S.S. (2005). The South China Sea margins: implications for rifting contrasts, Earth Planet Sci. Lett., 237: 601–616.

Metcalfe, I. (2011). Tectonic framework and phanerozoic evolution of Sundaland, Gondwana Res., 19 (1): 3–21.

How to cite: Guo, D.: Constrained Gravity Inversion for the Moho Depth and Tectonic Patterns in South China Sea, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4845, https://doi.org/10.5194/egusphere-egu24-4845, 2024.

Coseismic gravity changes provide significant information for the study of the mechanisms of large earthquakes and for developing fault models (Sun, 2012). In this research,  coseismic gravity changes of the 2008 Ms8.0 Wenchuan earthquake in China were studied by using gravity observation data and simulation based on a fault model.

Firstly, a fine processing of relative and absolute gravity data from the Longmenshan Gravimetric Network was carried out and observed gravity change of 22 stations near this earthquake were obtained; Secondly ,simulation of coseismic gravity changes was conducted based on half-space dislocation theory using the fault model obtained by Wang et al(2008) through inversion with multiple types of geodetic survey data, including GPS, INSAR, and leveling, and the results were compared with the observations..

It was found that the observed and simulated results are basically consistent, showing that the significant changes are mainly concentrated in the near-rupture zone in the hanging wall of the Yingxiu–Beichuan fault and that the changes decrease rapidly away from the rupture zone. The changes exhibit a positive to negative trend from east to west in the footwall of the Yingxiu–Beichuan fault and have a distribution characterized by alternate positive and negative changes in the hanging wall of the fault. This demonstrates the reliability of the observed results and the reasonableness of the fault model used in this paper.

In the near-rupture zone on the west and east sides of the Yingxiu–Beichuan fault, there are still some differences between the observed and simulated results. The trends in the spatial distribution of these differences exhibit a deviation similar to “phase delay”; in other words, an observed result deviates from the corresponding simulated result in terms of spatial position, which is speculated to be caused by errors in the geometric parameters and in the slip distribution of the fault model. After the slip distribution of  the Pengguan fault model was modified based on the actual surface rupture distribution, the simulated result at the Hongjiawan station near the eastern boundary of the fault model showed greater consistency with the observed result. This indicates that the observed gravity change results in this paper can provide an important reference for further detailed study of the fault model.         

            Fig1.Schematic of the Chengdu Gravimetric Network                        Fig2.Spatial distribution of observed gravity changes and simulated results

How to cite: Hao, H. and Hu, M.: Coseismic gravity changes of the 2008 Wenchuan earthquake in China observed by surface gravimetric data, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7561, https://doi.org/10.5194/egusphere-egu24-7561, 2024.

Multiple magnetic surveys, including fixed-wing, helicopter-borne, uncrewed aerial vehicle (UAV), and ground magnetics have been acquired over parts of the Bjerkreim-Sokndal layered intrusion (BKS) in Rogaland, Norway. The Proterozoic 230 km² Bjerkreim-Sokndal layered complex intrudes into anorthosites and hosts recurrent megacyclic units (MCU) with varying cumulus and critical minerals. Some MCUs are associated with strong magnetic remanence, resulting in Koenignsberger ratios (Q ratio) over 5 and anomalies of 12 000 nT below background.  A comparative analysis of these surveys over the Bjerkreim Lobe provide insights into what features can be mapped at different scales. Here we focus on new geological details provided by UAV, ground, and mineral scale surveys.  A UAV can typically operate at a maximum altitude of 150 m above terrain to a minimum of tens- of centimeters in ideal conditions. Thus, UAV magnetic surveys are optimal for understanding the change of a magnetic anomaly with varying source-separation through multiple flight altitudes. Survey altitudes by UAVs overlap with the source-sensor separation of ground and low-altitude crewed flights, therefore allowing a comparative analysis.

In 2023, the Norwegian University of Science and Technology and Helmholtz Institute Freiberg acquired coincident magnetic survey grids by UAV and ground magnetometer over key sites in the Bjerkreim lobe. Here we compare results of crewed-, uncrewed-, and ground-based data collected over the eastern margin of the Bjerkreim lobe and assess how these impact subsequent geologic interpretations. The petrophysical database for the survey area also contains > 1500 previously collected samples in combination with surface geometry information. This database in combination with the extensive lateral magnetic survey data at various sensor heights and other available complementary geophysics, including gravity, provide excellent parameters and constraints for forward modelling and inversions.

In the Bjerkreim lobe two MCU have significant magnetic remanence where anomalies are several thousand nT below background due to natural remanent magnetizations that are typically > 15 A/m. Therefore, an additional focus is on understanding the nature of the magnetic mineralogy using high-resolution scanning magnetic microscopy. These large amplitude magnetic anomalies may also cause logistical challenges for both airborne- and ground magnetic surveying. UAVs employ onboard magnetometer for navigation and attitude corrections which can be impacted by these large magnetic gradients. Similarly, significant noise or sensor drop-outs when the sensor’s dead zone is aligned with these large, often steep, magnetic gradients.

How to cite: McEnroe, S., Lee, M., Madriz, Y., Gloaguen, R., Pastore, Z., Lelièvre, P., and Church, N.: What can we learn from magnetic surveys at different scales? Geological insight from mineral to airborne surveys in the Bjerkreim-Sokndal Layer Intrusion, Norway, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8934, https://doi.org/10.5194/egusphere-egu24-8934, 2024.

High-precision observations of gravity plays a central role in modern approaches for active volcano monitoring. Time-lapse observations over a network of benchmarks are frequently used to detect underground mass redistribution in the plumbing system of active volcanoes. Such approach however does not allow to retrieve small mass variations occurring over short terms (i.e. few hours or days). To fill this gap, continuous gravity monitoring at a fixed station may be employed. In January 2023 the relative gravimeter gPhoneX#116 was installed at the WWF Nature Reserve of Astroni volcano, in the Campi Flegrei caldera, to further complement time-lapse observations periodically performed over a network of benchmarks. During the 1-year of recordings, the gPhone has continuously recorded the relative gravity changes, only shortly interrupted by a few technical issues. The purpose of the observations is to monitor continuously the short-term gravity signals in one of the world's highest-risk volcanoes; to pursue this objective targeted and meticulous corrections need to be applied to remove the effect of several other geophysical effects, such as tides and atmospheric effects, which may superpose on the signals of interest. Special effort was devoted to the study of instrumental drift, which can mask actual gravity changes due to mass variations occurring in the volcanic and geothermal systems. In this contribution we report the various processing steps and analysis performed to obtain reliable parameters of the Earth tides, non-tidal corrections and gravity residuals. The retrieved Earth tide model is then used to properly reduce tidal effects in high-precision relative and absolute gravity measurements.

How to cite: Fedele, A., Riccardi, U., Pivetta, T., Carlino, S., and Ricciardi, G.: A year-long gravity record at Astroni, Campi Flegrei (Southern Italy): some considerations on data processing for volcano monitoring and precise gravity tides, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17476, https://doi.org/10.5194/egusphere-egu24-17476, 2024.

Quantum sensors have gained increased attention in the last years, both from the applied but also from the manufacturer perspective. Most instruments are still limited to operation within a dedicated lab. With the Absolute Quantum Gravimeter (AQG), one device has already proven mobile capabilities and was used in several research studies. This application perspective is an important topic for new instruments, in order to meet scientific requirements in terms of usability and usefulness for various research interests.

One of these research interests is hydrology. From the monitoring perspective, hydrological observations in the field traditionally rely on point measurements, often in form of invasive sensor installations. These spatially-limited observations sometimes complicate natural hydrological process investigations. An advantage is provided when using the hydrogravimetric method, with its integral nature of monitoring water mass changes as a whole.

In this contribution, we address the above-mentioned important topics for a first feasibility study of an emerging instrument: the Differential Quantum Gravimeter (DQG). Developed by Exail Quantum Sensors, the DQG measures the acceleration due to gravity and the vertical gravity gradient simultaneously. It is an industry-grade demonstrator that has been operational for three years now and has achieved state-of-the-art sensitivity on the gradient of about 60E/sqrt(tau) and a long-term stability on the gradient around 1E. For gravity measurements the performances are on par or better than the AQG with a sensitivity of 600nm/s²/sqrt(tau) and a stability down to 5nm/s².

In preparation for first field measurements, we were interested in its performance for resolving hydrological dynamics and processes. We set up different hydrological modeling scenarios to forward model gravity responses and their DQG-related gravity gradients from water mass changes. Scenarios for obtaining these water mass changes consisted of vertical 1D models using the software Hydrus. Developing scenarios from very simple to more complex soil layer setups, we tested different forcing types (precipitation, evapotranspiration) with varying magnitudes and durations. The overall objective was to simulate resulting gravity gradients at different locations as the DQG would monitor them. Varying the theoretical placement of the DQG within the model domain enabled us to investigate its sensitivity to the simulated hydrological processes with respect to its location and distance to different magnitudes of water mass changes. Forward modeled data was averaged at different time periods and combined with realistically expected noise of the instrument. The study helps to evaluate the capabilities of the instrument as a tool to observe water fluxes in the soil, as well as optimal implementation of the DQG for planning first field measurements.

How to cite: Reich, M., Janvier, C., and Güntner, A.: Investigating the feasibility of resolving hydrological processes with a Differential Quantum Gravimeter by hydrological scenario modeling, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15133, https://doi.org/10.5194/egusphere-egu24-15133, 2024.

Asteroids and moons are promising targets for physical space exploration. The use of physically-based simulations within a virtual environment for (deep) space missions can significantly benefit the testing and validation of guidance, navigation, and control algorithms. This approach offers advantages in terms of cost and time efficiency. Especially for orbit propagation and landing maneuvers, information about the gravitational field is crucial. However, several factors contribute to the complexity of this task, such as limited information available about the inner structure of celestial bodies. The lack of detailed knowledge about their shapes further adds to the challenge.

This study presents a voxel-based mass concentration (MASCON) method to model detailed and realistic density distributions, enabling accurate gravity field determinations. We chose a cube with constant density as first case due to the perfect shape reconstruction and the availability of an analytical solution for its gravity field. To validate our results, we calculated the surface gravity and compared it with the analytical solution, ensuring the accuracy of our calculations. Furthermore, the surface gravity is derived for different resolutions and compared against other state-of-the-art methods like the polyhedral method that provides a closed-form analytical solution of the gravity field for homogeneous density. The other two methods for validation also use a MASCON approach, one utilizing polydisperse sphere packing and another with MASCON represented in spherical coordinates. The relative errors of the gravitational acceleration between the four methods will be evaluated for a cube and sphere, with homogeneous density.

The second aspect of this study was to create a tool that generates realistic density distributions. We are able to successfully reproduce natural environments by placing body-specific restrictions on three-dimensional Perlin noise with additional normalization. The simulator can add the following structural features to the density distribution: an arbitrary number of centralized or decentralized shells, with varying thickness and densities, anomalies of arbitrary size and shape, only restricted by its maximum permille of the body's volume. Furthermore, we implemented different normalization techniques to keep the mass of all generated bodies fixed. Our results show that the tool can generate realistic density distributions and calculate the corresponding gravitational field correctly. The data generated here is used to train Machine Learning and Deep Learning algorithms for gravity inversion.

How to cite: Haser, B., Andert, T., and Förstner, R.: Voxel-based Density Models for Accurate Gravitational Field Computation, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-1294, https://doi.org/10.5194/egusphere-egu24-1294, 2024.

Understanding the communication pathways of deep aquifers in the Tunisian Oriental Atlas along the southern Mediterranean margin, particularly within the coastal Sahel region, is of utmost importance for designing effective well drilling strategies and reducing risks for groundwater drilling.

In this study, we employed Gravity, Time Domain Electromagnetic (TDEM) methods and the variation of the piezometric level to investigate the structural setting and aquifer characteristic. Gravity and TDEM are two geophysical methods that provide insights into the density variation of underground bodies and reveal resistivity distribution at different depths, respectively. By integrating these methods, we aim to unravel the intricate hydrogeological system in the study area.

Our findings highlight a major fault line with a significant water level discrepancy, which is crucial information for groundwater exploration and exploitation. This study provides insights into the hydrogeological dynamics of the coastal Sahel region, facilitating the design of new drilling strategies. The gained knowledge supports informed decision-making in selecting optimal target production zones, ultimately minimizing drilling risks and promoting sustainable groundwater management in the Tunisian Oriental Atlas.

How to cite: Charrek, K., Walraevens, K., Hermans, T., and Gabtni, H.: Unravelling Deep Aquifer Communication in the Coastal Sahel Region: Insights from Geophysical Methods in the Tunisian Oriental Atlas, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16919, https://doi.org/10.5194/egusphere-egu24-16919, 2024.

Regional gravity and magnetic surveys are essential sources of information about the structure and geodynamics of the lithosphere. However, geologically meaningful inversion of gravity and magnetic data usually requires integration with other geophysical methods. We have developed a 3-D joint inversion framework that has the flexibility of using independent inversion codes and model discretizations for each of the included methods, is easily expandable and supports a wide range of the coupling constraints. Here we show its application to the regional geophysical datasets available in Ireland. We present the results of joint inversion of long-period magnetotelluric data, seismic traveltimes, and land gravity – a multiparameter geophysical model of the crust and uppermost mantle of the whole Ireland. On a smaller scale, we present the results of joint inversion of gravity, airborne magnetic and magnetotelluric data for the Limerick Basin, focusing on imaging of a Carboniferous volcanic structure.  The main aim is to better understand the Pb-Zn mineral systems which are controlled by the tectonics of the basement and lower crust. Exploration-scale geophysical surveys and geothermal exploration will also benefit from the regional 3-D geophysical models.

How to cite: Molodtsov, D., Kiyan, D., and Bean, C.: 3D joint inversion of regional magnetotelluric, seismic, gravity and magnetic datasets to image lithospheric structure of Ireland, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12072, https://doi.org/10.5194/egusphere-egu24-12072, 2024.

Sub-ice geology significantly influences the dynamics and future evolution of the Antarctic Ice Sheet, but largely inaccessible for direct sampling. Here, we present an approach, where we use a Self-Organizing Map (SOM) to describe sub-glacial properties. Based on attributes derived from gravity, magnetics and radar data from the NASA Operation Ice Bridge dataset in East Antarctica, we train a SOM, where attributes are selected to best represent sub-glacial conditions. Therefore, we study the trade-offs between these data sets helping to identify for which properties these are most sensitive.
The trained SOM identifies the outlines of the main geological structures beneath the ice and supplements models based on inverse and forward modelling. In contrast to such often regional interpretations, the SOM captures small-scale structures at the ice bed, as we illustrate with case examples, and highlights areas with inconsistencies in existing geological interpretations. The SOM can furthermore be used as input for inverse modelling of the physical properties of the sub-glacial geology in Antarctica.

How to cite: Ebbing, J., Liebsch, J., and Matsuoka, K.: Enhancing sub-ice geology in East Antarctica with Self-Organizing maps based on gravity, magnetic and radar data, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8059, https://doi.org/10.5194/egusphere-egu24-8059, 2024.