An interface for interoperability between the Generic Mapping Tools (GMT), a tool used by geophysicists to create research-quality figures, and Matlab has recently been developed that allows GMT users to interact with Matlab and Matlab users to make use of GMT.
GMT wrappers are currently also being developed for the Python programming language, particularly to be used in the IPython/Jupyter notebook due to an initiative by Leonardo Uieda (and his professor Paul Wessel) whose Postdoc is being funded by the NSF. You can watch his talk at the SciPy 2017 conference below.
Some of the mentioned advantages to which I attest include:
- Begin and End statements are introduced to eliminate the need to pipe postscript results into a file in each line of code being written. This also eliminates the need to use the -K and -O flags which keep the file open and updates it, respectively. The -K and -O flags are a major confusion for newcomers to GMT.
- temporary files are created under the /tmp directory, in Linux, so they will automatically be cleaned once the jupyter notebook is closed or the operating system is rebooted. Moreover, every project will have its own directory so files from different projects don’t get mixed up.
- GMT documentation straight in the Jupyter notebook
- Matplotlib- & Basemap-like behaviour, particularly inline viewing of figures, using gmt.show()
- Pythonic aliases make the compact GMT flags
To contribute: github.com/GenericMappingTools
Cook, T. (2017), A powerful new tool for research, Eos, 98, https://doi.org/10.1029/2017EO077489. Published on 17 July 2017.
One of the ways scientists are attempting to reduce greenhouse gases is to inject these gases into the ground. Specifically, they are testing injecting them, into basalt which is type of igneous rocks usually forming the first (rock) layer (sandwiched between the sedimentary & gabbro layers) in the oceanic crust basalt and in volcanic regions.
The above video features Iceland and its geothermal plants. Iceland is a heaven for geothermal energy as it lies along the Mid-Atlantic Ridge (MAR) (i.e. where Mid-Atlantic ocean floor is spreading apart in opposite directions forming a ridge). Most notably Iceland lies along the V-Shaped Reykjavik ridge (figure below; Google maps) which is part of the Norther MAR.
A “real-time” animation of the seismic (vertical) velocity of Italy’s 5.5 Mw earthquake that hit between the Aquila and Rieti provinces. The second half of the animation shows the whole country and surrounding area.
Real-time indicates that the video reflects how the waves propagated in real-time. It’s neither slowed down not sped up.
Red color indicates relatively higher vertical velocities indicating the ground is moving upwards while the blue color indicates the lower (negative) velocities indicated a downward movement of the ground. Color intensity refers to the magnitude.
In a previous post we selected some resources on fluid dynamics. In that post we suggested the videos from the National Committe for Fluid Mechanics (notes on the website) that date to the 1960’s and is an excellent exposition to fluid mechanics. We share them below following a very clear demonstration of Euler’s equation and Bernoulli’s integral:
What do you think an earthquake would sound like if you could hear it?
Here’s a sample (opens in a new tab).
The audio was generate by processing seismic data.
Credit: Peng, Z., C. Aiken*, D. Kilb, D. Shelly, B. Enescu (2012), Listening to the 2011 magnitude 9.0 Tohoku-Oki, Japan earthquake, Seismol. Res. Lett., 83(2), 287-293, doi: 10.1785/gssrl.83.2.287. , and Kilb, D., Z. Peng, D. Simpson, A. Michael and M. Fisher* (2012), Listen, watch, learn: SeisSound video products, Seismol. Res. Lett., 83(2), 281-286, doi: 10.1785/gssrl.83.2.281.
As a follow-up here’s a nice video from USGS