The antbgm from ll-geo

antbgm's Introduction

大家好！ Hello! 😉

I'm a postdoc and a geophysicist at the University of Western Australia, modelling the lithospheric structure of the North Australian Craton and trying to link petrophysics with geophysics. 🦘

I completed my PhD in Australia as well, 🐧, and used geophysics to understand the interaction between the solid earth and the cryosphere. It's quite challenging but really fun!!

Current work ⚡

🔥 Modelling the thermal-mechanical structure of Great McArthur Basin.

Previous work 🌱

🤽 Subglacial Sedimentary Basin distribution and groundwater in Antarctica (paper link)

🌐 Crustal Heterogeneity model in Antarctica (poster link) Check my thesis ..

🔥 Modelling the thermal-mechanical structure of Antarctic lithopshere. (My EGU 2023 talk) Check my thesis ..

Research interest ✨

🧊 Ice-sheet dynamics, subglacial hydrology and subglacial geology

🧐 Large scale inversion method and application

🔐 Open source software - check the amazing 🌎 Fatiando a Terra

Random Stuff 🙈

🗃️ PF toolbox: Some tutorial and pacakge about potentiel field

🗺️ ANTBGM: Personal compliation of (topography,) gravity and magnetic data in Antarctica

antbgm's People

Contributors

Stargazers

Watchers

antbgm's Issues

1. Large Scale Up/Downward continuation

Pointing out a gravity reference

Hey Lu, awesome idea to start compiling this stuff 😄. I just thought I'd point out a reference for a preliminary update to AntGG which I found recently: https://ftp.space.dtu.dk/pub/RF/4D-ANTARCTICA/.

They have some limited documentation in a pdf in that link.

Here's a download function I've been using for it:

def fetch_gravity(
    type: str, plot: bool = False, info: bool = False, region=None, spacing=10e3
) -> xr.DataArray:
    """
    Loads an Antarctic gravity grid
    Preliminary compilation of Antarctica gravity and gravity gradient data.
    Updates on 2016 AntGG compilation.
    Accessed from https://ftp.space.dtu.dk/pub/RF/4D-ANTARCTICA/.
    Parameters
    ----------
    type : str
        either 'FA' or 'BA', for free-air and bouguer anomalies, respectively.
    plot : bool, optional
        choose to plot grid, by default False
    info : bool, optional
        choose to print info on grid, by default False
    region : str or np.ndarray, optional
        GMT-format region to clip the loaded grid to, by default doesn't clip
    spacing : str or int, optional
        grid spacing to resample the loaded grid to, by default 10e3
    Returns
    -------
    xr.DataArray
        Returns a loaded, and optional clip/resampled grid of either free-air or
        Bouguer gravity anomalies.
    """

    if region is None:
        region = (-3330000, 3330000, -3330000, 3330000)
    path = pooch.retrieve(
        url="https://ftp.space.dtu.dk/pub/RF/4D-ANTARCTICA/ant4d_gravity.zip",
        known_hash=None,
        processor=pooch.Unzip(),
        progressbar=True,
    )
    file = [p for p in path if p.endswith(".dat")][0]
    df = pd.read_csv(
        file,
        delim_whitespace=True,
        skiprows=3,
        names=["id", "lat", "lon", "FA", "Err", "DG", "BA"],
    )
    transformer = Transformer.from_crs("epsg:4326", "epsg:3031")
    df["x"], df["y"] = transformer.transform(df.lat.tolist(), df.lon.tolist())
    df = pygmt.blockmedian(
        df[["x", "y", type]], spacing=spacing, region=region, verbose="q"
    )
    grd = pygmt.surface(
        data=df[["x", "y", type]], spacing=spacing, region=region, M="2c", verbose="q"
    )
    if plot is True:
        grd.plot(robust=True)
    if info is True:
        print(pygmt.grdinfo(grd))
    return grd