Creating Cross Sections

In scientific publications that deal with seismic tomographies, the results are often shown in 2D horizontal and vertical cross sections (slices and profiles). The Geophysical Model Generator therefore also provides tools do exactly this.

We will use the data that we have previously imported to show how to create cross sections. If you haven’t closed julia, these datasets should still be in memory, otherwise you have to load the corresponding jld2 file:

julia> load("Paffrath_Europe.jld2")

1. Creating horizontal cross sections (slices)

Once you have loaded the required data, creating a horizontal slice is quite easy using the GMG tool CrossSection. To extract a horizontal slice at a given depth, simply call the function as in the example below:

julia> Data_Paffrath_d100  =   CrossSection(Data_Paffrath, Depth_level=-100km)  

where the first argument to the function is the GMG data structure and the second argument denotes the kind of cross section to be extracted. The key word Depth_level denotes that a horizontal slice is to be extracted, with its depth value being -100 km (note that depth is negative).

2. Creating vertical cross sections (profiles)

Creating a vertical cross section is similarly straightforward, the only difference being that in this case, the starting and ending point of the cross section has to be provided:

julia> Data_Paffrath_A  =  CrossSection(Data_Paffrath, Start=(4.65,45.73), End=(17.23,43.80))

Where the start and end points have to be given as (longitude,latitude) pairs.


  1. Create cross-sections for all three profiles given below.
  2. If you have time, have a look for horizontal slices in Lippitsch et al. (2003).
  3. If you have additional time, , do the same for the P-wave tomography by Koulakov et al.
Profile ID lon start lat start lon end lat end Comment
A 4.65 45.73 17.23 43.80 Lippitsch 2002 & 2003 (location corrected)
B 5.51 51.53 12.04 43.68 Lippitsch 2002 & 2003 (location corrected)
C 17.78 50.95 11.66 43.68 Lippitsch 2002 & 2003

Marcel Thielmann
Marcel Thielmann
Principal Investigator/Group Leader

Working on localization processes in Earth materials, in particular deep earthquakes, flow of complex fluids in porous media and Earth’s lower mantle rheology.