DOI: 10.5285/4c8cfe1b-f426-4cc7-819b-236d1058ea0d
title: 3D density and susceptibility distribution model derived from 2024 gravity and magnetic surveys of Marguerite Bay, West Antarctica

description: We present a 3D crustal model of density and susceptibility distribution in Marguerite Bay, west Antarctica. The inversion is based on airborne gravity and magnetic data collected with a Windracers Ultra UAV (2023/2024 Antarctic season) with a line spacing of 2 km and a ground clearness of 500m.
The inverted densities and susceptibilities allow us to contain and identify large scale 3D intrusion in Marguerite Bay.
We provide a CSV file, which contains the coordinates, the depth and the inverted density and susceptibility distribution. 

This work was funded through "Protecting environments with unmanned aerial vehicle swarms" project number 10023377 as part of "Innovate UK Future flight challenge phase 3"
description: Crustal density and susceptibility distribution is inferred from "Variation of Information" joint inversion of airborne gravity and magnetic data. The inversion is carried out with the academic software JIF3D (https://sourceforge.net/projects/jif3d/) developed by Max Moorkamp.
The density and susceptibility inversion model are discretised into meshes with equal horizontal cell sizes of 250m. The vertical cell size is 100m at the surface. The vertical cell size increases with depth by a factor 1.1 for each cell. Increasing vertical cell size with depth is introduced to account for decreasing resolution with increasing distance to the source in potential field applications. 
Additionally, a padding area of 20% around the study area is added to avoid edge effects. The resulting inversion mesh contains 130 cells east-west direction, 314 cells in north-south direction and 33 cells in the z direction (20 km below sea level). To constrain the model geometry, the bedrock topography and Moho depths are used as model boundaries.
All cells are set to a starting value of 0 kg m-3 and 0 SI for density and susceptibility, respectively. The density and susceptibility values are iteratively updated during the joint inversion process.

Python code for preparing input data, setting up the inversion and plotting the results is openly accessible through zenodo.org (https://doi.org/10.5281/zenodo.14016296).

This dataset consists of 1 CSV file with 7 columns:
longitude (degree_east): longitude in degrees
latitude (degree_north): latitude in degrees
projection_x_coordinate (m): polar stereographic coordinates in meter. (straight_vertical_longitude_from_pole: 0; latitude_of_projection_origin: -90; standard_parallel: -71; false_easting: 0; false_northing: 0)
projection_y_coordinate (m): polar stereographic coordinates in meter. (straight_vertical_longitude_from_pole: 0; latitude_of_projection_origin: -90; standard_parallel: -71; false_easting: 0; false_northing: 0)
depth (meter, down = positive): Depths in meters. Positive values are defined as down.
density (kg.m-3): Density values in kg m-3
susceptibility (SI): susceptibility in SI
description: The data preparation and the inversion set up is performed in python (version 3.10.11). The inversion itself is carried out in the JIF3D software using a Linux terminal.
description: The mean residual between observed and inverted gravity and magnetic field is 0.1 mGal and -0.2 nT, with a standard deviation of 2.5 mGal and 3.5 nT.

The data have a horizontal resolution of 250 m and a vertical resolution of 100 m, which increases with a factor of 1.1 with depth.