Skip to content

Commit

Permalink
Update paper.md (#31)
Browse files Browse the repository at this point in the history
* Update paper.md

A few grammatical updates to the paper.

* Better description of the motivation for using `pycurious`

- Added a figure.png
- Updated some references

* Fixing grammar suggestions

* Fixing citation styles

* Whedon doesn't like the colons in each reference

* Richer image caption
  • Loading branch information
lheagy authored and Ben Mather committed Jul 24, 2019
1 parent a257379 commit 385cf6e
Show file tree
Hide file tree
Showing 3 changed files with 45 additions and 15 deletions.
Binary file added paper/figure.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
14 changes: 13 additions & 1 deletion paper/paper.bib
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ @article{Tanaka1999
number = {3-4},
pages = {461--470},
title = {{Curie point depth based on spectrum analysis of the magnetic anomaly data in East and Southeast Asia}},
url = {http://linkinghub.elsevier.com/retrieve/pii/S0040195199000724},
url = {https://dx.doi.org/10.1016/S0040-1951(99)00072-4},
volume = {306},
year = {1999}
}
Expand All @@ -35,3 +35,15 @@ @article{Meyer2017
urldate = {2018-09-01},
year = {2017}
}

@article{Mather2019,
author = {Mather, Ben and Fullea, Javier},
doi = {10.5194/se-10-839-2019},
journal = {Solid Earth},
number = {3},
pages = {839--850},
title = {{Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data}},
url = {https://doi.org/10.5194/se-10-839-2019},
volume = {10},
year = {2019}
}
46 changes: 32 additions & 14 deletions paper/paper.md
Original file line number Diff line number Diff line change
Expand Up @@ -24,11 +24,28 @@ bibliography: paper.bib

# Summary

Magnetic data is one of the most prevelant geophysical datasets available on the surface of the Earth. Curie depth is the depth at which rocks lose their magnetism. The most prevalent magnetic mineral is magnetite, which has a Curie point of 580°C, thus the Curie depth is often interpreted as the 580°C isotherm. Current methods to derive Curie depth first compute the (fast) Fourier transform over a square window of a magnetic anomaly that has been reduced to the pole. The depth and thickness of magnetic sources is estimated from the slope of the radial power spectrum.

Our Python package, `pycurious`, ingests maps of the Earth's magnetic anomaly and distributes the computation of Curie depth across multiple CPUs. `pycurious` implements the [@Tanaka1999] and [@Bouligand2009] methods for computing the thickness of a buried magnetic source. The former selects portions of the radial power spectrum in the low and high frequency domain to compute the depth of magnetic sources, while the latter fits an analytical solution to the entire power spectrum. We cast the [@Bouligand2009] method within a Bayesian framwork to estimate the uncertainty of Curie depth calculations. Common computational workflows and geospatial manipulation of magnetic data are covered in the Jupyter notebooks bundled with this package. The `mapping` module includes a set of functions that help to wrangle maps of the magnetic anomaly into a useable form for `pycurious`. Such an approach is commonly encountered for transforming global compilations of the magnetic anomaly, e.g. EMAG2 [@Meyer2017], from latitudinal/longitudinal coordinates to a local projection in eastings/northings.

<!-- `pycurious` is an object-oriented Python package that accepts grids of the Earth's magnetic anomaly and computes the Curie depth across a specified window size. The computation is distributed across all available processors to improve efficiency over large study areas. -->
Multiple geophysical methods have been proposed to resolve the thermal structure of the Earth's lithosphere with varying degrees of precision.
Geotherms may be constructed from heat flow or xenolith data, but the spatial coverage of these are often limited.
Seismic velocity has proven effective to infer upper-mantle temperature, but its application relies on building a compositional model suitable for the geological context and estimating attenuation from grainsize and water content.
In contrast, magnetic data are among the most widespread geophysical datasets available on the surface of the Earth and can offer useful insight into its thermal structure.

The magnetic anomaly is the observed magnetic field less the contribution from the core and external fields.
Geologic features enhance or depress the local magnetic field; this response is primarily controlled by the concentration of magnetite in the Earth's crust.
Magnetite is the most prevalent magnetic mineral, in terms of both its magnetic susceptibility and concentration in the crust, and has a Curie point of 580°C.
The depth at which rocks lose their permanent magnetic properties is referred to as the Curie depth.
Due to the prevalence of magnetite in the crust, the Curie depth is often interpreted as the 580°C isotherm.
Current methods to derive the Curie depth first compute the (fast) Fourier transform over a square window of a magnetic anomaly that has been reduced to the pole.
The depth and thickness of magnetic sources is estimated from the slope of the radial power spectrum.

Our Python package, `pycurious`, ingests a map of the Earth's magnetic anomaly and distributes the computation of Curie depth across multiple CPUs.
`pycurious` implements the @Tanaka1999 and @Bouligand2009 methods for computing the thickness of a buried magnetic source.
The former selects portions of the radial power spectrum in the low and high frequency domain to compute the depth of magnetic sources, while the latter fits an analytical solution to the entire power spectrum.
We cast the @Bouligand2009 method within a Bayesian framework to estimate the uncertainty of Curie depth calculations [@Mather2019].
Common computational workflows and geospatial manipulation of magnetic data are covered in the Jupyter notebooks bundled with this package.
The `mapping` module includes a set of functions that help to wrangle maps of the magnetic anomaly into a useful form for `pycurious`.
Such an approach is commonly encountered for transforming global compilations of the magnetic anomaly, e.g. EMAG2 [@Meyer2017], from latitudinal/longitudinal coordinates to a local projection in eastings/northings.

![Radial power spectra (right) computed from different sized windows of a synthetic magnetic anomaly (left) using `pycurious`. Larger window sizes increase the number of points in the low wavenumber range of the radial power spectrum, which is important for estimating the thickness of magnetic sources in the crust.](figure.png)

`pycurious` includes the following functionality:

Expand All @@ -39,27 +56,28 @@ Our Python package, `pycurious`, ingests maps of the Earth's magnetic anomaly an
- Bayesian inversion framework and sensitivity analysis.


## Documentation
## Documentation and examples

`pycurious` is bundled with a linked collection of jupyter notebooks that can act as a user guide and an introduction to the package. The notebooks are split into matching sets for frequentist and Bayesian estimates of Curie depth. The notebooks cover:
`pycurious` is bundled with a linked collection of Jupyter notebooks that can act as a user guide and an introduction to the package.
The notebooks are split into matching sets for frequentist and Bayesian estimates of Curie depth.
The notebooks cover:

- Plotting the radial and azimuthal power spectrum.
- Computing Curie depth from a synthetic magnetic anomaly.
- Exploring parameter sensitivity.
- Posing the inverse problem and objective function.
- Mapping Curie depth using the EMAG2 magnetic anomaly dataset [@Meyer2017].

All documentation can be accessed from within the module via a Python function that installs the notebooks at a filesystem location specified by the user at run time. The API documentation is kept up-to-date on [GitHub pages](https://brmather.github.io/pycurious/).


## Installation, Dependencies and Usage

`pycurious` requires `numpy` and `Cython` to compile C routines that are included with the distribution. The documentation is supplied in the form of Jupyter notebooks (the Jupyter environment is a dependency) which also have optional dependencies for the `mapping` module including `matplotlib`, `pyproj`, and `cartopy`. `pycurious` and all Python dependencies can be installed through the pypi.org `pip` package manager, however, several of the dependencies for `cartopy` may cause problems for inexperienced users. We therefore provide a Docker image built with all required and optional dependencies and a deployment of the documentation / examples on [mybinder.org](https://mybinder.org/v2/gh/brmather/pycurious/binder?filepath=Notebooks%2F0-StartHere.ipynb).
All documentation can be accessed from within the module via a Python function that installs the notebooks at a filesystem location specified by the user at run time.
The API documentation is kept up-to-date on [GitHub pages](https://brmather.github.io/pycurious/).
`pycurious` and all Python dependencies can be installed through the pypi.org `pip` package manager, however, several of the dependencies for `cartopy` may cause problems for inexperienced users.
We therefore provide a Docker image built with all required and optional dependencies and a deployment of the documentation / examples on [mybinder.org](https://mybinder.org/v2/gh/brmather/pycurious/binder?filepath=Notebooks%2F0-StartHere.ipynb).


# Acknowledgements

This work was made possible by the G.O.THERM.3D project, supported by an Irish Research Council Research for Policy & Society grant (RfPS/2016/50) co-funded by Geological Survey Ireland and by Sustainable Energy Authority Of Ireland. Development of `pycurious` was inspired from [`pycpd`](https://github.com/groupeLIAMG/pycpd) and [`fatiando a terra`](https://github.com/fatiando/fatiando).
This work was made possible by the G.O.THERM.3D project, supported by an Irish Research Council Research for Policy & Society grant (RfPS/2016/50) co-funded by Geological Survey Ireland and by Sustainable Energy Authority Of Ireland.
Development of `pycurious` was inspired from [`pycpd`](https://github.com/groupeLIAMG/pycpd) and [`fatiando a terra`](https://github.com/fatiando/fatiando).


# References

0 comments on commit 385cf6e

Please sign in to comment.