PROGRAM SHIFT_GSRM2p2x Program to create several versions (x = a, b, c, d, e, f,) of a long-term, indefinite-window global shallow seismicity map based on the Global Strain Rate Map, version 2.2 (GSRM2.2) created in 2016 by Corne Kreemer. The methods of scaling strain-rate to seismicity progress from simple to complex in the various versions, with the more complex versions using the set of methods known as Seismic Hazard Inferred From Tectonics (SHIFT), as presented by [Bird & Liu, 2007, Seismol. Res. Lett.] and adapted to global strain-rate maps by Bird et al. [2010, Seismol. Res. Lett.] and by Bird & Kreemer [2015, Bull. Seismol. Soc. Amer.]. By Peter Bird, UCLA, August 2019, for the SHIFT_GSRM2.2x project. Output will be in Peter Bird's .GRD format, in units of epicenters/m^2/s. Now, opening log file SHIFT_GSRM2.2x_log.txt. Creating global grids of 899 rows * 1440 columns = 1294560 points for input strain rate tensors, output epicentral rate densities, & mean cz. Total grid area = 5.10064E+14 square meters, compared to ideal sphere area of 5.10064E+14 square meters. Difference is due to 2 polar caps of 3.88E+08 square meters each. According to the GSRM2.2 data: Plate boundary area = 7.09330E+13 m^2 (13.907%). Intraplate area = 4.39131E+14 m^2 (86.093%). and we can assume the 2 polar caps are intraplate, as well. Enter name of (global, shallow) seismic catalog .EQC file for calibration: GCMT_shallow_m5p767_1977-2013.eqc Here are the first and last earthquakes in the catalog: GlobalCMT 1977.01.01 19:01:40.7 126.150 -2.880 43 6.35 3 46 83 157 7 316 GlobalCMT 2013.12.29 15:37:40.5 126.330 1.410 29 5.78 2 122 17 31 73 217 What is the duration of this catalog, in years? 37.0000000000000 What is the completeness-threshold magnitude of this catalog? 5.76700000000000 Catalog threshold (minimum) scalar moment = 5.018E+17 N m This threshold will be used for calibration of the forecast, regardless of the requested threshold for forecast reporting. PLEASE SELECT DESIRED MODEL VERSION FROM THE FOLLOWING LIST: 1(a): Uniform (calibrated) coupled thickness in all boundaries [2 DOF]. 2(b): Add isotropic smoothing of offshore non-SUB boundaries [2 DOF]. 3(c): Add seaward smoothing of SUB boundaries to outer rises [2 DOF]. 4(d): Add separate calibrations in tectonic zones #0~4 [5 DOF]. 5(e): Add classification & partitioning of strain rate tensors [5 DOF]. 6(f): Add velocity-dependence of SUB and CCB coupling [5 DOF]. CHOOSE INTEGER CODE: 6 Enter desired magnitude threshold for output forecast: 5.76700000000000 Forecast threshold (minimum) scalar moment = 5.018E+17 N m Create a .GRD file-name for the output: SHIFT_GSRM2.2f_m5.767_2014plus.grd File underwater.grd must be available. Please check this. File PB2002_steps.dat must be available. Please check this. PB2002_steps.dat has been converted to a list of 81114 internal points. Finding PB2002 plate boundary nearest to each active GSRM2 grid cell... Search completed. Isotropically smoothing rates adjacent to offshore non-SUB plate boundaries... Isotropic smoothing of offshore non-SUB plate boundaries completed. Smoothing SUB rates toward the outer rise regions... Smoothing of SUB rates toward outer rises completed. File PB2002_tectonic_zones.grd must be available. Please check this. Tectonic zones grid has 1801 rows of 3601 columns. and was stored as INTEGER*1 to keep memory down to ~6 MB. Finished selecting tectonic zones for GSRM2 grid from finer input grid. Inferring zone #s (> 0) for active cells in zone 0, by referring to neighbors. Iteration 1: 79993 problems; 14333 assigned (through common edges). Iteration 2: 65660 problems; 12823 assigned (through common edges). Iteration 3: 52837 problems; 11210 assigned (through common edges). Iteration 4: 41627 problems; 7527 assigned (through common edges). Iteration 5: 34100 problems; 4862 assigned (through common edges). Iteration 6: 29238 problems; 3532 assigned (through common edges). Iteration 7: 25706 problems; 2736 assigned (through common edges). Iteration 8: 22970 problems; 2297 assigned (through common edges). Iteration 9: 20673 problems; 1940 assigned (through common edges). Iteration 10: 18733 problems; 1651 assigned (through common edges). Iteration 11: 17082 problems; 1414 assigned (through common edges). Iteration 12: 15668 problems; 1180 assigned (through common edges). Iteration 13: 14488 problems; 1035 assigned (through common edges). Iteration 14: 13453 problems; 883 assigned (through common edges). Iteration 15: 12570 problems; 778 assigned (through common edges). Iteration 16: 11792 problems; 699 assigned (through common edges). Iteration 17: 11093 problems; 633 assigned (through common edges). Iteration 18: 10460 problems; 578 assigned (through common edges). Iteration 19: 9882 problems; 524 assigned (through common edges). Iteration 20: 9358 problems; 492 assigned (through common edges). Iteration 21: 8866 problems; 461 assigned (through common edges). Iteration 22: 8405 problems; 438 assigned (through common edges). Iteration 23: 7967 problems; 416 assigned (through common edges). Iteration 24: 7551 problems; 387 assigned (through common edges). Iteration 25: 7164 problems; 356 assigned (through common edges). Iteration 26: 6808 problems; 340 assigned (through common edges). Iteration 27: 6468 problems; 328 assigned (through common edges). Iteration 28: 6140 problems; 305 assigned (through common edges). Iteration 29: 5835 problems; 293 assigned (through common edges). Iteration 30: 5542 problems; 272 assigned (through common edges). Iteration 31: 5270 problems; 255 assigned (through common edges). Iteration 32: 5015 problems; 238 assigned (through common edges). Iteration 33: 4777 problems; 214 assigned (through common edges). Iteration 34: 4563 problems; 198 assigned (through common edges). Iteration 35: 4365 problems; 181 assigned (through common edges). Iteration 36: 4184 problems; 173 assigned (through common edges). Iteration 37: 4011 problems; 167 assigned (through common edges). Iteration 38: 3844 problems; 163 assigned (through common edges). Iteration 39: 3681 problems; 162 assigned (through common edges). Iteration 40: 3519 problems; 158 assigned (through common edges). Iteration 41: 3361 problems; 154 assigned (through common edges). Iteration 42: 3207 problems; 149 assigned (through common edges). Iteration 43: 3058 problems; 132 assigned (through common edges). Iteration 44: 2926 problems; 129 assigned (through common edges). Iteration 45: 2797 problems; 124 assigned (through common edges). Iteration 46: 2673 problems; 119 assigned (through common edges). Iteration 47: 2554 problems; 117 assigned (through common edges). Iteration 48: 2437 problems; 110 assigned (through common edges). Iteration 49: 2327 problems; 85 assigned (through common edges). Iteration 50: 2242 problems; 77 assigned (through common edges). Iteration 51: 2165 problems; 71 assigned (through common edges). Iteration 52: 2094 problems; 69 assigned (through common edges). Iteration 53: 2025 problems; 70 assigned (through common edges). Iteration 54: 1955 problems; 69 assigned (through common edges). Iteration 55: 1886 problems; 72 assigned (through common edges). Iteration 56: 1814 problems; 73 assigned (through common edges). Iteration 57: 1741 problems; 76 assigned (through common edges). Iteration 58: 1665 problems; 76 assigned (through common edges). Iteration 59: 1589 problems; 67 assigned (through common edges). Iteration 60: 1522 problems; 57 assigned (through common edges). Iteration 61: 1465 problems; 52 assigned (through common edges). Iteration 62: 1413 problems; 42 assigned (through common edges). Iteration 63: 1371 problems; 43 assigned (through common edges). Iteration 64: 1328 problems; 43 assigned (through common edges). Iteration 65: 1285 problems; 45 assigned (through common edges). Iteration 66: 1240 problems; 45 assigned (through common edges). Iteration 67: 1195 problems; 45 assigned (through common edges). Iteration 68: 1150 problems; 47 assigned (through common edges). Iteration 69: 1103 problems; 46 assigned (through common edges). Iteration 70: 1057 problems; 45 assigned (through common edges). Iteration 71: 1012 problems; 46 assigned (through common edges). Iteration 72: 966 problems; 45 assigned (through common edges). Iteration 73: 921 problems; 44 assigned (through common edges). Iteration 74: 877 problems; 44 assigned (through common edges). Iteration 75: 833 problems; 44 assigned (through common edges). Iteration 76: 789 problems; 44 assigned (through common edges). Iteration 77: 745 problems; 43 assigned (through common edges). Iteration 78: 702 problems; 38 assigned (through common edges). Iteration 79: 664 problems; 36 assigned (through common edges). Iteration 80: 628 problems; 36 assigned (through common edges). Iteration 81: 592 problems; 35 assigned (through common edges). Iteration 82: 557 problems; 24 assigned (through common edges). Iteration 83: 533 problems; 24 assigned (through common edges). Iteration 84: 509 problems; 26 assigned (through common edges). Iteration 85: 483 problems; 27 assigned (through common edges). Iteration 86: 456 problems; 27 assigned (through common edges). Iteration 87: 429 problems; 29 assigned (through common edges). Iteration 88: 400 problems; 28 assigned (through common edges). Iteration 89: 372 problems; 28 assigned (through common edges). Iteration 90: 344 problems; 30 assigned (through common edges). Iteration 91: 314 problems; 22 assigned (through common edges). Iteration 92: 292 problems; 21 assigned (through common edges). Iteration 93: 271 problems; 20 assigned (through common edges). Iteration 94: 251 problems; 18 assigned (through common edges). Iteration 95: 233 problems; 17 assigned (through common edges). Iteration 96: 216 problems; 16 assigned (through common edges). Iteration 97: 200 problems; 15 assigned (through common edges). Iteration 98: 185 problems; 13 assigned (through common edges). Iteration 99: 172 problems; 12 assigned (through common edges). Iteration 100: 160 problems; 11 assigned (through common edges). Iteration 101: 149 problems; 10 assigned (through common edges). Iteration 102: 139 problems; 9 assigned (through common edges). Iteration 103: 130 problems; 7 assigned (through common edges). Iteration 104: 123 problems; 6 assigned (through common edges). Iteration 105: 117 problems; 4 assigned (through common edges). Iteration 106: 113 problems; 4 assigned (through common edges). Iteration 107: 109 problems; 2 assigned (through common edges). Iteration 108: 107 problems; 0 assigned (through common edges). Now giving up and zeroing strain rates in 107 orphan cells. Largest strain rate zeroed out was 3.61E-17/s ( 1.1 nanostrain/year). Finding PB2002 plate boundary nearest to each newly-active GSRM2 grid cell... Search completed. Tectonic zone 1: area 4.7061E+13 m^2 ( 9.23% of Earth) Tectonic zone 2: area 2.0909E+13 m^2 ( 4.10% of Earth) Tectonic zone 3: area 1.5306E+13 m^2 ( 3.00% of Earth) Tectonic zone 4: area 4.2205E+13 m^2 ( 8.27% of Earth) Intraplate : area 3.8458E+14 m^2 ( 75.40% of Earth) Earthquakes of magnitude >= 5.767 in plate boundary zone 1: 880 ( 13.59%) Earthquakes of magnitude >= 5.767 in plate boundary zone 2: 427 ( 6.59%) Earthquakes of magnitude >= 5.767 in plate boundary zone 3: 644 ( 9.94%) Earthquakes of magnitude >= 5.767 in plate boundary zone 4: 4442 ( 68.59%) Earthquakes of magnitude >= 5.767 in plate interiors: 83 ( 1.28%) Intraplate seismicity rate (per m^2, per s) = 1.84834E-22 Global plate-boundary seismicity scaling factor in zone 1 = 0.91325. Global plate-boundary seismicity scaling factor in zone 2 = 2.39170. Global plate-boundary seismicity scaling factor in zone 3 = 2.33521. Global plate-boundary seismicity scaling factor in zone 4 = 3.10455. Writing seismicity (EQs/m^2/s) .grd file and computing area-integral... Area integral of seismicity = 5.5463E-06 earthquakes/s = 1.7503E+04 earthquakes/century. Writing effective_cz_meters .grd file... Writing moment_factors_in_SI.dat file... Total seismic moment rate in SI = 4.6240E+14 N m /s. Job completed.