The QLG program can be used to determine an average or to perform a tomographic inversion. The method is described in Ottemöller et al. [2002]. Here, we use the same names for the damping parameters, and many of the other parameters should be self-explanatory. The program can also produce the input for distance trace plots. Note that using the program is no trivial task. The data set needs to be carefully selected and the instrument calibration has to be known. The input to the program is a Nordic file, which includes several events. The parameter file needs to be carefully set up.

The program can be used in the following way:

- Determine average
- Perform checker-board test to chose damping parameters
- Tomographic inversion

Note: The main purpose of including the program is to give an example source code so that the user can make use of it when implementing similar programs. The program uses a linear grid...

Example of the parameter file `qlg.par` :

KEYWORD............Comments.............Par 1.....Par 2..... FILTER for distance plot 0.01 15. DISTANCES min and max 200. 3000. GROUP VEL LG lg group vel window 3.0 3.7 GROUP VEL P p group vel window 5.0 8.0 INVERSION TYPE 1. for tomography 1. 0. for average ORIENTATION 0.=vert, 1.=horiz. 0. PHASE ONLY 1.=phase pick requ. 1. in s-file FREQUENCY frequency and 1/q 1. 5.2392E-03 FREQUENCY frequency and 1/q 1.25 4.5246E-03 FREQUENCY frequency and 1/q 1.60 4.1239E-03 FREQUENCY frequency and 1/q 2. 3.5312E-03 FREQUENCY frequency and 1/q 2.5 2.9081E-03 FREQUENCY frequency and 1/q 3.15 2.2568E-03 FREQUENCY frequency and 1/q 4. 1.7029E-03 FREQUENCY frequency and 1/q 5. 1.1228E-03 STATION MIN min # of stations 4. VELOCITY LG 3350. DAMPING ALPHA damping parameters 500. DAMPING SIGMA ------------------ 100. DAMPING BETA ------------------ 500. DAMPING LAMBDA ------------------ 0.001 NSMOOTH smooth spec # times 0. CHECKERBOARD 1. for cb-test 0. CHECKERBOARD DELTA 0.0004 FIX SITE 0. FIX SOURCE 0. SOURCE PERTURBATION 7. 0.2 GAUSSIAN NOISE 0.1 VERBOSE 0. for quite mode 1. # # Grid # X START x start of grid -92.5 Y START y start of grid 6.50 X DELTA x delta grid 1. Y DELTA y delta grid 1. X NPOINTS x # points 17. Y NPOINTS y # points 13.

Menke et al. [2006] pointed out the non-uniqueness in attenuation tomography between the source term and Q. They suggest to investigate the non-uniqueness by synthetic tests in which a perturbation is applied to the source term and the inversion for Q is done without inverting for differences in the source term. The solutions obtained are null-solutions and one needs to be careful not to mistake them for real patterns. These tests are possible within QLG by setting the parameter `SOURCE PERTURBATION, where the first parameter refers to the source that is perturbed and the second parameter gives the amount of perturbation in units of moment magnitude.

It is possible to invert real data without inverting for the site term by setting `FIX SITE'. This can be a useful test as there is also a trade-off with the site term. Fixing the site term is more problematic, as this is done based on the local magnitude, which may not be the same as the moment magnitude.

Another useful stability test is to add Gaussian noise to the spectra and check the inversion result. This can be done for both real data and the checkerboard test by setting the parameter `GAUSSIAN NOISE', units are equivalent to change in moment magnitude.