Stack

MSNoise is capable of using a reference function defined by absolute or relative dates span. For example, an absolute range could be “from 1 January 2010 to 31 December 2011” and a relative range could be “the last 200 days”. In the latter case, the REF will need to be exported at every run, meaning the following steps (MWCS and DTT) will be executed on the whole configured period. If the REF is defined between absolute dates, excluding “today”, the MWCS and DTT will only be calculated for new data (e.g. “yesterday” and “today”). The corresponding configuration bits are ref_begin and ref_end. In the future, we plan on allowing multiple references to be defined.

Only data for new/modified dates need to be exported. If any CC-job has been marked “Done” within the last day and triggered the creation of STACK jobs, the stacks will be calculated and a new MWCS job will be inserted in the database. For dates in the period of interest, the moving-window stack will only be exported if new/modified CCF is available. The export directory are “REF/” and “DAY%03i/” where %03i will be replaced by the number of days stacked together (DAYS_005 for a 5-days stack, e.g.).

Please note that within MSNoise, stacks are always inclusive of the time/day mentioned. For example, a 5-days stack on January 10, will contain cross-correlation functions computed for January 6, 7, 8, 9 AND 10! The graphical representation centered on a “January 10” tick might then display changes in the CCF that occurred on the 10th !

Moving-window stacks are configured using the mov_stack parameter in msnoise admin.

If stack_method is ‘linear’, then a simple mean CFF of all daily is saved as the mov or ref CCF. On the other hand, if stack_method is ‘pws’, then all the Phase Weighted Stack (PWS) is computed and saved as the mov or ref CCF. The PWS is done in two steps: first the mean coherence between the instantaneous phases of all windows is calculated, and eventually serves a weighting factor on the mean. The smoothness of this weighting array is defined using the pws_timegate parameter in the configuration. The weighting array is the power of the mean coherence array. If pws_power is equal to 0, a linear stack is done (then it’s faster to do set stack_method = ‘linear’). Usual value is 2.

Warning

PWS is largely untested, not cross-validated. It looks good, but that doesn’t mean a lot, does it? Use with Caution! And if you cross-validate it, please let us know!!

Schimmel, M. and Paulssen H., “Noise reduction and detection of weak, coherent signals through phase-weighted stacks”. Geophysical Journal International 130, 2 (1997): 497-505.

Configuration Parameters

  • ref_begin : Beginning or REF stacks. Can be absolute (2012-01-01) or relative (-100) days (default=1970-01-01)

  • ref_end : End or REF stacks. Same as ref_begin (default=2100-01-01)

  • mov_stack : A list of two parameters: the time to “roll” over (default 1 day) and the granularity (step) of the resulting stacked CCFs (default 1 day) to stack for the Moving-window stacks. This can be a list of tuples, e.g. ((‘1d’,’1d’),(‘2d’,’1d’)) corresponds to the MSNoise 1.6 “1,2” before. Time deltas can be anything pandas can interpret (“d”, “min”, “sec”, etc). (default=((‘1d’,’1d’)))

  • stack_method : Stack Method: Linear Mean or Phase Weighted Stack (default=linear) | new in 1.4

  • pws_timegate : If stack_method =’pws’, width of the smoothing (in seconds) (default=10) | new in 1.4

  • pws_power : If stack_method =”pws”, Power of the Weighting (default=2) | new in 1.4

  • hpc : Is MSNoise going to run on an HPC? (default=N) | new in 1.6

Once done, each job is marked “D”one in the database and, unless hpc is Y, MWCS jobs are inserted/updated in the database.

Usage:

msnoise cc stack --help

Usage:  [OPTIONS]

  Stacks the [REF] or [MOV] windows. Computes the STACK jobs.

Options:
  -r, --ref   Compute the REF Stack
  -m, --mov   Compute the MOV Stacks
  -s, --step  Compute the STEP Stacks
  --help      Show this message and exit.

For most users, the REF stack will need to be computed only once for specific dates and then, on routine basis, only compute the MOV stacks:

$ msnoise cc stack -r
$ msnoise reset STACK
$ msnoise cc stack -m

as for all other steps, this procedure can be run in parallel:

$ msnoise -t 4 cc stack -r
$ msnoise reset STACK
$ msnoise -t 4 cc stack -m

New in version 1.4: The Phase Weighted Stack.

New in version 1.6: The hpc parameter that can prevent the automatic creation of MWCS jobs. The REF and MOV stacks have been separated and need to be run independently.