SAGA-Task2.7 Modelling Software Table MDabas - Modelling Software
SAGA-Task2.7 Modelling Software Table MDabas - Modelling Software
geophysics.
The subject is related to the diversity of software available for forward and inverse
modelling in geophysics with a specific focus for archaeological applications. It is
intended to make a list of them in a first step (some work is already available by the
author of this call and some wiki pages) and a comparison between them with a focus
on their use for detectability of archaeological structures (prediction of magnitude and
shape of anomalies) and their usability (environment, ergonomy, etc.). The most used
methods (magnetics, resistivity, low frequency EM and GPR) are the main methods
looked for. This study does not limit itself to open-source software, but also includes
freeware software. Visualization software is not part of the scope of this call.
ELECTROMAGNETICS
Despite instruments used for a long time (Geonics EM31 for example), this method is probably
the one which has the least software available. Among the different responses obtained from
our questionnaire, only two software packages were quoted (Inv2DVLF and pyGIMLi).
- Inv2DVLF is an inversion code for single frequency VLF data (we will not discuss in
detail VLF data, because the method is very rarely used in geophysics applied to
archaeology). It was developed by Monteiro Santos et al. in 2016 following the
resolution of the direct problem using finite element modelling by Sasaki (ref:
http://emtomo.com/news-and-events/42-a-new-version-of-the-vlf2d-em-
multifrequency).
- A new version for multifrequency data is now available: VLF2Dmf
(http://emtomo.com/products/vlf2dmf). It includes a 2D forward modelling and 2D
Occam inversion of VLF-EM data including topography. This software is not free, but
we mention it because of the paucity of available LFEM software.
- EM4 soil from the same company (EMTOMO) is a program for 1D Laterally
Constrained Inversion (or quasi-2D inversion) and 1D Spatially Constrained Inversion
(quasi-3D inversion) of data collected with multisensor and multifrequency instruments
like those from DUALEM, GEONICS, GF-CMD or GEM and GSSI-Profiler. This
software is not free (http://emtomo.com/products/em4soil).
- We also want to mention the code from Hanssens (Ghent University). He has made it
available on GitHub (https://github.com/dhanssens/FDEM1D) in November 2018. It is
a 1D forward and sensitivity modelling code for a 1D multi-layered-space.
Conductivity, susceptibility and permittivity of the different layers can be changed.
Code is written both in Python 3 and MatlabR2016 and necessitates one of these
software packages.
General comments:
Most of the researchers we know use their own routines. It must be pointed that the full
resolution of Maxwell equations is difficult and some simplifications has been used for
a long time following the work of Geonics in Canada. Using these simple equations (for
example the cumulative response for Horizontal coplanar or vertical coplanar
configurations) is straightforward but is an approximation (Born). To our knowledge
only the team of A. Tabbagh/ J. Thiesson in Paris-Sorbonne University and J.
Guillemoteau / Tronicke in Potsdam University have resolved the full problem for
archaeological applications (Guillemoteau et al., 2016). The Paris team has made pieces
of software available (written in Fortran77) on request. There is no graphical interface
and input and output of data is not straightforward. We have no information about
Guillemoteau’s software. Most of the software quoted is only modelling 1D cases with
lateral constraints making them quasi-2D case. The Aarhus software seems to us one of
the best available software packages even if we have not tested it.
N.B. To keep updated about new open-source software browse for example:
https://github.com/search?q=geophysics+electromagnetics .
SEISMICS
- ReflexW is a program for processing and interpretation of GPR data, seismic and
ultrasound data with different geometry configurations (reflection, refraction, borehole
cross-hole and tomography - 2D and 3D cases). The program is very stable and has been
refined since more than 20 years by Sandmeier (REFLEXW is written with Delphi
5/DelphiXE2).
(https://www.sandmeier-geo.de/reflexw.html). The software is licenced and includes a
discount for academic licences and also free updates for 5 years after purchase. It runs
only under Windows. There exist many pdf guides and videos.
General comments:
We did not spend much time to browse most of the available software in seismics
because of the rare use of this method for archaeological investigations. ReflexW is a
long-time stable software perfectly suited for 2D and 3D analysis that can also be used
for GPR processing.
- MG2DATA: This open-source code was developed by Stocco et al., 2009 in Matlab. It
allows the forward model and the inversion of 2D and 2.5D magnetic data of both Total,
vertical gradient and vertical component of the magnetic field. Remanent magnetization
can be included. This code can only be used for selected magnetic profiles extracted
from maps for example to estimate the size and depth of elongated source bodies. This
code applies to 2 D and 2.5D structures (2.5D means that the extension of the bodies
perpendicular to the profile is not ‘infinite’ – practically more than 10 times the
dimension on-line – but finite).
- Potensoft developed by Arisoy and Dikmen, 2011 is a Matlab code (Version 7.6 Release
2008b) not only for modelling but also processing and mapping gravity and magnetic
data. It runs under Windows but has been tested with Linux. It is easy to use (a well-
designed GUI) and is open-source. (Unfortunately the link:
http://www.eng.ankara.edu.tr/~arisoy/potensoft.htm is no longer valid).
- POTENT is a software developed by Geophysical Software Solutions Pty. Ltd. and runs
under Windows (http://www.geoss.com.au/potent.html). It provides a highly interactive
framework for 3-D modelling of magnetic and gravity data. PotentQ, is a simplified and
streamlined version of Potent that allows rapid semi-automatic modelling of a single
magnetic and/or gravity anomaly (PotentQ has been integrated with Geosoft's Oasis
montaj interface, https://www.geosoft.com/products ). Considering the price for Oasis,
we will not discuss this very powerful software. The price of POTENT is also high for
a perpetual licence but 30 days or 6 months licences are available.
GPR
- GRED HD 3D: It is also not a modelling software. The GRED HD software is able to
process data collected from all IDS GPR systems (https://idsgeoradar.com/). Data are
acquired by OneVision software and processed by GRED HD. Objects of any generic
shape are able to be directly exported to CAD or GIS maps.
General comments:
There exists much free code for the modelling of GPR anomalies. A Python/Matlab
environment is generally needed. gprMax seems to be the only software able to deal
with complex EM properties and volume geometries like those that we encounter in
archaeological prospecting. But its use is complex and there is no graphical interface at
this time. MatGPR is certainly simpler to use – but restricted to 2D cases- with a
standalone version (without Matlab). We have no information if this project is still alive.
- DC-2DPro / 3DPro: this software for 2D and 3D inversion is developed by Prof. Jung-
Ho Kim from KIGAM, the Korean Institute of Geoscience and Mineral Resources (Kim,
2009). It is based on a Finite-Element Method to solve the forward resistivity problem
and an active constraint balancing approach incorporating L1 and L2 regularization for
fine inverse solutions. It used to be commercially available. The status right now is not
known since Dr. Jung-Ho Kim has recently passed away (Papadopoulos, pers. com.).
General comments:
Despite the fact that resistivity imaging is not the most frequently used method in
geophysics applied to archaeology (possibly due to considerable work load), there exists
much free code for the modelling of resistive anomalies. This is due to the wide use of
ERT in other fields of geophysics and due to the work of Loke and Papadopoulos for
example: ready-to-use executable codes from 1D to 3D exists. Otherwise, a Python
environment is generally needed for more complex and scalable applications
P.S.
We have recently found a promising freeware: SimPEG.XYZ (Simulation and
Parameter Estimation in Geophysics, https://simpeg.xyz/). It is an open source Python
package for simulation (finite volume) and gradient based parameter estimation in
geophysical applications (Cocket et al., 2015). This software is based on a MIT Licence
(this license does not force packages that use SIMPEG to be open source nor does it
restrict commercial use). It provides simulation and inversion tools for all the methods
we have discussed in this paper. SimPEG is a very active community with a forum, lots
of documentation, examples of codes and weekly online meetings.
We must thank now Armin Schmidt for his careful reading. He has mentioned software
packages that we missed to mention:
-1 The broad range of software available from the UBC Geophysical Inversion Facility,
the team of Doug Oldenburg: https://gif.eos.ubc.ca/software/main_programs. These
software (DCIP2D and 3D/MAG3D/GRAV3D/EM1DFM/EM1DTM) are available for
academic, educational and commercial licensing. SimPEG seems to be derived from
their previous expertise. A very interesting Internet site with on-line courses, toolkit,
books, etc. including SimPEG is available at http://geosci.xyz.
-2 Seismic Un*x (often called ‘SU’) for the processing of GPR and seismic reflection
data (https://wiki.seismic-unix.org/doku.php ). The CWP/SU is an open source seismic
utilities package originally created by the Center for Wave Phenomena (CWP) at the
Colorado School of Mines (CSM) and is more than 50 years old. One must have a Unix
or Unix-like environment to run SU. This package offers more than 450 processing and
auxiliary tools….
I am sure I miss lots of packages here, but this list is open now to new contributors…
__________________________________________________________________________________
REFERENCES
Arisoi, Özgü, Dikmen, Ünal, 2011. Computer and Geosciences Potensoft. Computers & Geosciences 37,
7, 935-942.
Auken, E., Christiansen, A. V., Kirkegaard, C., Fiandaca, G., Schamper, C., Behroozmand, A. A., Binley,
A., Nielsen, E., Effersø, F., Christensen, N. B., Sørensen, K. I., Foged, N., and Vignoli, G., 2014. An
overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and
borehole electromagnetic and electric data, Exploration Geophysics.
Binley A, 2015. Tools and Techniques: DC Electrical Methods, In: Treatise on Geophysics, 2nd Edition,
G Schubert (Ed.), Elsevier, 11, 233-259.
Blanchy G., Saneiyan S., Boyd J., McLachlan P. and Binley A., 2020. “ResIPy, an Intuitive Open Source
Software for Complex Geoelectrical Inversion/Modelling.” Computers & Geosciences, February,
104423.
Cockett, Rowan, Seogi Kang, Lindsey J. Heagy, Adam Pidlisecky, and Douglas W. Oldenburg, 2015.
SimPEG: An Open Source Framework for Simulation and Gradient Based Parameter Estimation in
Geophysical Applications. Computers & Geosciences, September 2015.
doi:10.1016/j.cageo.2015.09.015.
Giannakis, I., Giannopoulos, A., Warren, C. & Davidson N., 2015. An Accurate FDTD Numerical
Modelling Framework of Ground Penetrating Radar for Landmine DetecKon, IEEE Journal Of Selected
Topics In Applied Earth ObservaKons And Remote Sensing.
Goodman D, 1994. Ground-penetrating radar simulation in engineering and archaeology. Geophysics,
59, 2, 224-232.
Götze, H.-J. and Lahmeyer, B., 1988. Application of three-dimensional interactive modelling in gravity
and magnetics: Geophysics, 53, 8, 1096-1108.
Guillemoteau, J., Simon, F.-X., Lück, E., Tronicke, J., 2016. 1D sequential inversion of portable multi-
configuration electromagnetic induction data. Near Surface Geophysics. 14, 411-420.
Günther, T., Rücker, C., and Spitzer, K., 2006. 3-d modelling and inversion of DC resistivity data
incorporating topography - Part II: Inversion. Geophys. J. Int., 166, 2, 506-517.
Hanssens, D., Delefortrie, S., De Pue, J., Van Meirvenne, M., and P. De Smedt, 2019. Frequency-Domain
Electromagnetic Forward and Sensitivity Modelling: Practical Aspects of modelling a Magnetic Dipole
in a Multilayered Half-Space. IEEE Geoscience and Remote Sensing Magazine, 7, 1, 74-85
Huber E. and Hans G., 2018. RGPR — An open-source package to process and visualize GPR data. 17th
International Conference on Ground Penetrating Radar (GPR), Switzerland, Rapperswil, 18-21 June
2018, pp. 1-4. doi: 10.1109/ICGPR.2018.8441658.
Kim, J., 2009. DC2DPro : 2D interpretation system of DC resistivity tomography. User's Manual and
Theory. Publisher: KIGAM, S. Korea.
Monterio Santos FA, Mateus A, Figueiras J, Goncalves MA, 2006. Mapping groundwater contamination
around a landfill facility using the VLF-EM method- a case study. Journal of Applied Geophysics, 60,
115-125.
Papadopoulos, N.G., Tsourlos P., Papazachos C., Tsokas G.N., Sarris A. and Kim J.H., 2011. An algorithm
for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried
antiquities. Geophysical Prospecting, 59, 557-575.
Rücker, C., Günther, T., Wagner, F.M., 2017. pyGIMLi: An open-source library for modelling and
inversion in geophysics, Computers and Geosciences, 109, 106-123, doi: 10.1016/j.cageo.2017.07.011.
Sasaki Y., 1989. Two-dimensional joint inversion of magnetotelluric and dipole–dipole resistivity data,
Geophysics, 54, 254-262
Schmidt S., Götze Hans-Jürgen, Fichler C., Alvers M., 2010. 3D Gravity, FTG and magnetic modelling:
the new IGMAS + software. Geoinformatik 2010.
Schmidt S., Götze H.-J., Fichler Ch., Ebbing J., Alvers M. R., 2010. IGMAS+ a new 3D Gravity, FTG and
Magnetic Modelling Software, In: Zipf A., Behncke, K.,Hillen F., Schfermeyer J. (Hsrg.): Die Welt im Netz
(Geoinformatik 2010 Kiel 17.3 – 19.3.2010).
Stocco S., Godio A., Sambuelli L., 2009. Modelling and compact inversion of magnetic data: A Matlab
code. Computers and Geosciences, 35, 2111–2118.
Tszanis A., MATGPR Release 2: A freeware MATLAB® package for the analysis & interpretation of
common & single offset GPR data, fastTIMES, 15, 1, 2010, 15-43.
Uieda, L, Oliveira Jr, V C, Ferreira, A, Santos, H B; Caparica Jr, J F, 2014. Fatiando a Terra: a Python
package for modelling and inversion in geophysics. figshare. doi:10.6084/m9.figshare.1115194.
Warren, C., Giannopoulos, A., & Giannakis, I., 2016. gprMax: Open source software to simulate
electromagnetic wave propagation for Ground Penetrating Radar, Computer Physics Communications,
209, 163-170.
Geophysical Modelling Software Informations
Geophysical Method
Software Name Modelling Forward
Owner Web Site User Mode Source Reference
dimension solution
Electrical Resistivity - Induced Polarization
InvVERIS Emtomo Commercial Executable
https://www.veristech.com/the-sensors 1
pyGIMLi-BERT Günther T, Rücker C, Wagner Free Open source
https://www.pygimli.org/index.html
F 2 and 3 Finite-differenceGünther,
and finite-element
T., Rücker, C., and Spitzer, K. (2006). 3-d modeling and inversion
RES1D Loke MH Free
http://www.geotomosoft.com/ Executable 1 Linear filter method
RES2Dmod /3Dmod Loke MH Free
http://www.geotomosoft.com/ Executable 2 Finite-difference and finite-element
RES2DINV 3DINV http://www.geotomosoft.com
Loke MH Commercial Executable 2 and 3 Finite-difference and finite-element
DC_2DPro /3DPro KIGAM (Korean Commercial
http://www.kigam.re.kr
Institute of Geoscience andExecutable 2 and 3 Finite-element
Mineral Resources) Kim, J., 2009. DC2DPro-2D interpretation system of DC resistivity tomograp
SensInv2D/3D Crosswellinstruments.de Commercial Executable
http://geotomographie.de/ 2 and 3 SIRT/MSIRT
2DIP/TS2DIP Zonge Commercial Executable
http://www.zonge.com 2 Finite-element
EarthImager2D / 3D Advanced Geosciences Commercial Executable
http://www.agiusa.com/
Inc. 2 and 3
R2/cR2/R3t/cR3t /ResIPy Binley A Free Executable 2 and 3 Finite-difference
http://www.es.lancs.ac.uk/people/amb/Freeware/Freeware.htm Binley
andA,finite-element
2015. Tools and Techniques: DC Electrical Methods, In: Treatise o
RESINVM3D SEG Free Open source
https://github.com/xiaoyingpu/RESINVM3D-mod 3
3DInv http://www.ims.forth.gr/index_main.php?c=90&l=e&d=7
Free Open source
GPR
ReflexW Commercial Executable
Sandmeier KJhttps://www.sandmeier-geo.de/reflexw.html 2 and 3 Finite-Difference Time-Domain
GPRSIM Commercial Executable
Goodman D https://www.gpr-survey.com/gprsim.html 2 Ray tracing Goodman D, 1994. Ground-penetrating radar simulation in engineering an
GPRSlice Commercial Executable
Goodman D https://www.gpr-survey.com 2 and 3
gprMax Free
Giannopouloshttps://www.gprmax.com
A Open source 2 and 3 Finite-DifferenceWarren,
Time-Domain
C., Giannopoulos, A., & Giannakis, I. (2016). gprMax: Open source
VIY3 Viy https://viy.ua Free Executable 2 and 3
GRED HD 3D IDS corp. https://idsgeoradar.com/
Commercial Executable
MATGPR Tzanis http://users.uoa.gr/~atzanis/matgpr
Open source
Free and licenced 2 Finite-Difference
Tszanis A., MATGPR Release 2: A freeware MATLAB® package for the analy
RGPR Huber http://emanuelhuber.github.io/RGPR
Free Open source 2 Huber E. and Hans G. (2018) RGPR — An open-source package to process
Geophysical Method
Software Name Modelling Forward
Owner Web Site User Mode Source Reference
dimension solution
POTENT Software Geophysical software solutions
Commercial
(Australia)
Executable 2 and 3D
GM-SYS Commercial
Northwest Geophysical Associates Inc (US)
Executable 2D
GMM G-Trend Ltd. (Slovakia) ? ? 2 and 2.5 ?
MOD3D Igor Cerovskyhttps://github.com/igorcerovsky/Mod3D
Free Open source 3D formulae for polyhedral bodies
IGMAS+ Univ. Kiel http://www.gravity.uni-kiel.de/igmas
? Executable 3D IGMAS+
Integral eq. Over a new 3D gravity, FTG and magnetic modeling software
polyhedra
GMG Uieda Free
https://btozer.github.io/gmgOpen source 2D Talwani eq. Uieda, L, Oliveira Jr, V C, Ferreira, A, Santos, H B; Caparica Jr, J F (2014), Fa
Seismics
ReflexW Commercial Executable 2 and 3
Sandmeier KJhttps://www.sandmeier-geo.de/reflexw.html Finite-Difference Time-Domain
Seismic Un*x Colorado School
https://wiki.seismic-unix.org/doku.php
of Mines/CWP Free Open source 1 to 3
Electromagnetics
Inv2DVLF Monteiro Santos Licenced Executable 2
http://emtomo.com/news-and-events/42-a-new-version-of-the-vlf2d-em-multifrequency
FA Finite-ElementMonteiro Santos FA, António Mateus, Jorge Figueiras, Mário A. Gonçalves,
pyGIMLi Günther T, Rücker
https://www.pygimli.org/index.html
C, Wagner Free
F Open source 1 to 3 Finite-difference
Rücker
andC,finite-element
Günther T, Wagner FM, 2017. pyGIMLi: An open-source library fo
trie, S., De Pue, J., Van Meirvenne, M., and
FDEM1D P. Ghent
Univ. De Smedt,
https://github.com/dhanssens/FDEM1D
2019. Frequency-Domain
(Belgium) Free Open
Electromagnetic
source 1ForwardFinite-difference
and Sensitivity Modeling: Practical Aspects of modeling a Magnetic Dipole in a Multila
and finite-element
General software
UBC Edu. Licence Open source
https://gif.eos.ubc.ca/software/main_programs
DCIP2D and 3D/MAG3D/GRAV3D/EM1DFM/EM1DTM 1 to 3
SIMPEG.XYZ Community https://simpeg.xyz/Free Open source 1 to 3 Finite-ElementCockett, Rowan, Seogi Kang, Lindsey J. Heagy, Adam Pidlisecky, and Dougla
PyGIMLl Günther T, Rücker
https://www.pygimli.org/index.html
C, Wagner Free
F Open source 1 to 3 Finite-difference
Rücker
andC,finite-element
Günther T, Wagner FM, 2017. pyGIMLi: An open-source library fo