diff --git a/paper/figure.png b/paper/figure.png new file mode 100644 index 0000000..d8b4cc1 Binary files /dev/null and b/paper/figure.png differ diff --git a/paper/paper.bib b/paper/paper.bib index 2548d40..ce54825 100644 --- a/paper/paper.bib +++ b/paper/paper.bib @@ -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} } @@ -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} +} diff --git a/paper/paper.md b/paper/paper.md index 48a553d..7ce064d 100644 --- a/paper/paper.md +++ b/paper/paper.md @@ -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. - - +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: @@ -39,9 +56,11 @@ 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. @@ -49,17 +68,16 @@ Our Python package, `pycurious`, ingests maps of the Earth's magnetic anomaly an - 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