Difference between revisions of "Registering CS173 into UCVM"

From SCECpedia
Jump to navigationJump to search
 
(41 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 
This page details the steps involved in delivering a velocity mesh to John Shaw and Andreas Plesch, so that the Harvard group's Central California basins can be integrated into a velocity model, which can be fed back into UCVM and used for CyberShake.
 
This page details the steps involved in delivering a velocity mesh to John Shaw and Andreas Plesch, so that the Harvard group's Central California basins can be integrated into a velocity model, which can be fed back into UCVM and used for CyberShake.
 +
 +
== Test Points ==
 +
To evaluate the registration of CS17.3 into UCVM we are querying a set of point using Scott's original CS17.3 and the recently updated UCVMC version. The following test points are sites from the CS17.3 study. We will query each point at 5 depths 0km, 100m, 1km, 10km, and 40km
 +
 +
<pre>
 +
<name>Osito Audit: Castaic Lake Dam (OSI)</name>
 +
<coordinates>-118.723500,34.614500</coordinates>
 +
 +
<name>Santa Barbara (SBR)</name>
 +
<coordinates>-119.698189,34.420831</coordinates>
 +
 +
<name>s1356 (s1356)</name>
 +
<coordinates>-120.582740,34.968060</coordinates>
 +
 +
<name>Calstate Bakersfield (BAK)</name>
 +
<coordinates>-119.104450,35.344440</coordinates>
 +
 +
<name>Lime Mtn. Quarry (LMQ)</name>
 +
<coordinates>-120.990669,35.668835</coordinates>
 +
 +
<name>Parkfield (PARK)</name>
 +
<coordinates>-120.432800,35.899700</coordinates>
 +
 +
<name>Huron (HUR)</name>
 +
<coordinates>-120.102917,36.202731</coordinates>
 +
 +
<name>Lemoore (LEM)</name>
 +
<coordinates>-119.782911,36.300783</coordinates>
 +
 +
</pre>
 +
 +
== Test points and depths list ==
 +
<pre>
 +
-118.723500 34.614500 0.0
 +
-118.723500 34.614500 100.0
 +
-118.723500 34.614500 1000.0
 +
-118.723500 34.614500 10000.0
 +
-118.723500 34.614500 40000.0
 +
-119.698189 34.420831 0.0
 +
-119.698189 34.420831 100.0
 +
-119.698189 34.420831 1000.0
 +
-119.698189 34.420831 10000.0
 +
-119.698189 34.420831 40000.0
 +
-120.582740 34.968060 0.0
 +
-120.582740 34.968060 100.0
 +
-120.582740 34.968060 1000.0
 +
-120.582740 34.968060 10000.0
 +
-120.582740 34.968060 40000.0
 +
-119.104450 35.344440 0.0
 +
-119.104450 35.344440 100.0
 +
-119.104450 35.344440 1000.0
 +
-119.104450 35.344440 10000.0
 +
-119.104450 35.344440 40000.0
 +
-120.990669 35.668835 0.0
 +
-120.990669 35.668835 100.0
 +
-120.990669 35.668835 1000.0
 +
-120.990669 35.668835 10000.0
 +
-120.990669 35.668835 40000.0
 +
-120.432800 35.899700 0.0
 +
-120.432800 35.899700 100.0
 +
-120.432800 35.899700 1000.0
 +
-120.432800 35.899700 10000.0
 +
-120.432800 35.899700 40000.0
 +
-120.102917 36.202731 0.0
 +
-120.102917 36.202731 100.0
 +
-120.102917 36.202731 1000.0
 +
-120.102917 36.202731 10000.0
 +
-120.102917 36.202731 40000.0
 +
-119.782911 36.300783 0.0
 +
-119.782911 36.300783 100.0
 +
-119.782911 36.300783 1000.0
 +
-119.782911 36.300783 10000.0
 +
-119.782911 36.300783 40000.0
 +
</pre>
 +
 +
 +
 +
ucvm_query result from UCVM-C 18.5 alpha
 +
 +
<pre>
 +
../bin/ucvm_query -f ../conf/ucvm.conf -m cs173 << EOF &> cs173-it.log
 +
-118.723500 34.614500 0.0
 +
...
 +
EOF
 +
</pre>
 +
 +
<pre>
 +
WARNING: Could not load model into memory. Reading the model from the
 +
hard disk may result in slow performance.
 +
Using Geo Depth coordinates as default mode.
 +
-118.7235    34.6145      0.000    713.454    390.000      cs173  2047.660  1089.404  2139.522      none      0.000      0.000      0.000      crust  2047.660  1089.404  2139.522
 +
-118.7235    34.6145    100.000    713.454    390.000      cs173  2734.936  1518.460  2232.448      none      0.000      0.000      0.000      crust  2734.936  1518.460  2232.448
 +
-118.7235    34.6145  1000.000    713.454    390.000      cs173  5460.268  3237.855  2605.311      none      0.000      0.000      0.000      crust  5460.268  3237.855  2605.311
 +
-118.7235    34.6145  10000.000    713.454    390.000      cs173  6561.615  3871.144  2890.995      none      0.000      0.000      0.000      crust  6561.615  3871.144  2890.995
 +
-118.7235    34.6145  40000.000    713.454    390.000      cs173  7545.088  4324.382  3141.390      none      0.000      0.000      0.000      crust  7545.088  4324.382  3141.390
 +
-119.6982    34.4208      0.000    15.086    387.000      cs173  4670.389  2585.928  2430.952      none      0.000      0.000      0.000      crust  4670.389  2585.928  2430.952
 +
-119.6982    34.4208    100.000    15.086    387.000      cs173  4668.188  2583.502  2430.503      none      0.000      0.000      0.000      crust  4668.188  2583.502  2430.503
 +
-119.6982    34.4208  1000.000    15.086    387.000      cs173  4625.225  2526.679  2420.277      none      0.000      0.000      0.000      crust  4625.225  2526.679  2420.277
 +
-119.6982    34.4208  10000.000    15.086    387.000      cs173  5915.015  3283.019  2621.902      none      0.000      0.000      0.000      crust  5915.015  3283.019  2621.902
 +
-119.6982    34.4208  40000.000    15.086    387.000      cs173  7619.421  4327.448  3143.100      none      0.000      0.000      0.000      crust  7619.421  4327.448  3143.100
 +
-120.5827    34.9681      0.000    18.443    280.000      cs173  4767.372  2397.728  2399.038      none      0.000      0.000      0.000      crust  4767.372  2397.728  2399.038
 +
-120.5827    34.9681    100.000    18.443    280.000      cs173  4771.665  2403.488  2399.941      none      0.000      0.000      0.000      crust  4771.665  2403.488  2399.941
 +
-120.5827    34.9681  1000.000    18.443    280.000      cs173  4833.515  2478.506  2412.065      none      0.000      0.000      0.000      crust  4833.515  2478.506  2412.065
 +
-120.5827    34.9681  10000.000    18.443    280.000      cs173  5818.455  3107.401  2560.882      none      0.000      0.000      0.000      crust  5818.455  3107.401  2560.882
 +
-120.5827    34.9681  40000.000    18.443    280.000      cs173  7883.487  4559.606  3267.956      none      0.000      0.000      0.000      crust  7883.487  4559.606  3267.956
 +
-119.1044    35.3444      0.000    114.259    280.000      cs173  3372.060  1940.204  2335.459      none      0.000      0.000      0.000      crust  3372.060  1940.204  2335.459
 +
-119.1044    35.3444    100.000    114.259    280.000      cs173  3378.782  1943.201  2335.872      none      0.000      0.000      0.000      crust  3378.782  1943.201  2335.872
 +
-119.1044    35.3444  1000.000    114.259    280.000      cs173  3507.057  1995.574  2343.018      none      0.000      0.000      0.000      crust  3507.057  1995.574  2343.018
 +
-119.1044    35.3444  10000.000    114.259    280.000      cs173  6903.631  4181.306  3061.306      none      0.000      0.000      0.000      crust  6903.631  4181.306  3061.306
 +
-119.1044    35.3444  40000.000    114.259    280.000      cs173  7822.489  4575.821  3276.215      none      0.000      0.000      0.000      crust  7822.489  4575.821  3276.215
 +
-120.9907    35.6688      0.000    602.291    782.000      cs173  4317.442  2337.356  2389.908      none      0.000      0.000      0.000      crust  4317.442  2337.356  2389.908
 +
-120.9907    35.6688    100.000    602.291    782.000      cs173  4324.663  2342.934  2390.741      none      0.000      0.000      0.000      crust  4324.663  2342.934  2390.741
 +
-120.9907    35.6688  1000.000    602.291    782.000      cs173  4501.126  2466.457  2410.151      none      0.000      0.000      0.000      crust  4501.126  2466.457  2410.151
 +
-120.9907    35.6688  10000.000    602.291    782.000      cs173  5709.418  3078.566  2551.708      none      0.000      0.000      0.000      crust  5709.418  3078.566  2551.708
 +
-120.9907    35.6688  40000.000    602.291    782.000      cs173  7870.905  4547.270  3261.620      none      0.000      0.000      0.000      crust  7870.905  4547.270  3261.620
 +
-120.4328    35.8997      0.000    467.078    349.000      cs173  4388.868  2370.534  2394.871      none      0.000      0.000      0.000      crust  4388.868  2370.534  2394.871
 +
-120.4328    35.8997    100.000    467.078    349.000      cs173  4390.334  2371.916  2395.081      none      0.000      0.000      0.000      crust  4390.334  2371.916  2395.081
 +
-120.4328    35.8997  1000.000    467.078    349.000      cs173  4371.586  2367.767  2394.453      none      0.000      0.000      0.000      crust  4371.586  2367.767  2394.453
 +
-120.4328    35.8997  10000.000    467.078    349.000      cs173  5795.831  3301.412  2628.905      none      0.000      0.000      0.000      crust  5795.831  3301.412  2628.905
 +
-120.4328    35.8997  40000.000    467.078    349.000      cs173  7586.191  4296.916  3126.092      none      0.000      0.000      0.000      crust  7586.191  4296.916  3126.092
 +
-120.1029    36.2027      0.000    113.530    280.000      cs173  2851.014  1495.107  2259.639      none      0.000      0.000      0.000      crust  2851.014  1495.107  2259.639
 +
-120.1029    36.2027    100.000    113.530    280.000      cs173  2859.429  1498.931  2260.468      none      0.000      0.000      0.000      crust  2859.429  1498.931  2260.468
 +
-120.1029    36.2027  1000.000    113.530    280.000      cs173  3050.430  1580.739  2277.315      none      0.000      0.000      0.000      crust  3050.430  1580.739  2277.315
 +
-120.1029    36.2027  10000.000    113.530    280.000      cs173  5612.277  3107.839  2561.259      none      0.000      0.000      0.000      crust  5612.277  3107.839  2561.259
 +
-120.1029    36.2027  40000.000    113.530    280.000      cs173  7696.131  4393.316  3179.477      none      0.000      0.000      0.000      crust  7696.131  4393.316  3179.477
 +
-119.7829    36.3008      0.000    69.541    280.000      cs173  3031.683  1489.972  2258.508      none      0.000      0.000      0.000      crust  3031.683  1489.972  2258.508
 +
-119.7829    36.3008    100.000    69.541    280.000      cs173  3040.815  1494.782  2259.555      none      0.000      0.000      0.000      crust  3040.815  1494.782  2259.555
 +
-119.7829    36.3008  1000.000    69.541    280.000      cs173  3253.721  1604.038  2281.744      none      0.000      0.000      0.000      crust  3253.721  1604.038  2281.744
 +
-119.7829    36.3008  10000.000    69.541    280.000      cs173  5852.581  3303.596  2629.943      none      0.000      0.000      0.000      crust  5852.581  3303.596  2629.943
 +
-119.7829    36.3008  40000.000    69.541    280.000      cs173  7657.497  4359.726  3160.989      none      0.000      0.000      0.000      crust  7657.497  4359.726  3160.989
 +
</pre>
 +
 +
== Another set of Test Points ==
 +
<pre>
 +
-118.72254 34.61444 0
 +
-118.72254 34.61444 175
 +
-118.72254 34.61444 1050
 +
-118.72254 34.61444 9975
 +
-118.72254 34.61444 40075
 +
-119.69846 34.42187 0
 +
-119.69846 34.42187 175
 +
-119.69846 34.42187 1050
 +
-119.69846 34.42187 9975
 +
-119.69846 34.42187 40075
 +
-120.58229 34.96759 0
 +
-120.58229 34.96759 175
 +
-120.58229 34.96759 1050
 +
-120.58229 34.96759 9975
 +
-120.58229 34.96759 40075
 +
-119.10517 35.34494 0
 +
-119.10517 35.34494 175
 +
-119.10517 35.34494 1050
 +
-119.10517 35.34494 9975
 +
-119.10517 35.34494 40075
 +
-120.99114 35.66917 0
 +
-120.99114 35.66917 175
 +
-120.99114 35.66917 1050
 +
-120.99114 35.66917 9975
 +
-120.99114 35.66917 40075
 +
-120.43380 35.89934 0
 +
-120.43380 35.89934 175
 +
-120.43380 35.89934 1050
 +
-120.43380 35.89934 9975
 +
-120.43380 35.89934 40075
 +
-120.50437 35.81076 0
 +
-120.50437 35.81076 175
 +
-120.50437 35.81076 1050
 +
-120.50437 35.81076 9975
 +
-120.50437 35.81076 40075
 +
-119.78270 36.30157 0
 +
-119.78270 36.30157 175
 +
-119.78270 36.30157 1050
 +
-119.78270 36.30157 9975
 +
-119.78270 36.30157 40075
 +
</pre>
 +
 +
 +
ucvm_query result
 +
 +
<pre>
 +
-118.7225    34.6144      0.000    720.976    390.000      cs173  2043.192  1085.995  2138.201      none      0.000      0.000      0.000      crust  2043.192  1085.995  2138.201
 +
-118.7225    34.6144    175.000    720.976    390.000      cs173  3248.173  1837.993  2301.267      none      0.000      0.000      0.000      crust  3248.173  1837.993  2301.267
 +
-118.7225    34.6144  1050.000    720.976    390.000      cs173  5459.535  3237.474  2605.180      none      0.000      0.000      0.000      crust  5459.535  3237.474  2605.180
 +
-118.7225    34.6144  9975.000    720.976    390.000      cs173  6556.426  3867.625  2889.135      none      0.000      0.000      0.000      crust  6556.426  3867.625  2889.135
 +
-118.7225    34.6144  40075.000    720.976    390.000      cs173  7551.173  4330.115  3144.578      none      0.000      0.000      0.000      crust  7551.173  4330.115  3144.578
 +
-119.6985    34.4219      0.000    16.241    387.000      cs173  4656.098  2577.254  2429.342      none      0.000      0.000      0.000      crust  4656.098  2577.254  2429.342
 +
-119.6985    34.4219    175.000    16.241    387.000      cs173  4652.338  2573.052  2428.572      none      0.000      0.000      0.000      crust  4652.338  2573.052  2428.572
 +
-119.6985    34.4219  1050.000    16.241    387.000      cs173  4609.664  2514.169  2418.086      none      0.000      0.000      0.000      crust  4609.664  2514.169  2418.086
 +
-119.6985    34.4219  9975.000    16.241    387.000      cs173  5913.173  3281.206  2621.216      none      0.000      0.000      0.000      crust  5913.173  3281.206  2621.216
 +
-119.6985    34.4219  40075.000    16.241    387.000      cs173  7622.714  4330.343  3144.709      none      0.000      0.000      0.000      crust  7622.714  4330.343  3144.709
 +
-120.5823    34.9676      0.000    18.570    280.000      cs173  4765.085  2395.066  2398.626      none      0.000      0.000      0.000      crust  4765.085  2395.066  2398.626
 +
-120.5823    34.9676    175.000    18.570    280.000      cs173  4772.822  2405.407  2400.245      none      0.000      0.000      0.000      crust  4772.822  2405.407  2400.245
 +
-120.5823    34.9676  1050.000    18.570    280.000      cs173  4837.920  2483.522  2412.900      none      0.000      0.000      0.000      crust  4837.920  2483.522  2412.900
 +
-120.5823    34.9676  9975.000    18.570    280.000      cs173  5814.321  3105.312  2560.201      none      0.000      0.000      0.000      crust  5814.321  3105.312  2560.201
 +
-120.5823    34.9676  40075.000    18.570    280.000      cs173  7883.104  4559.139  3267.717      none      0.000      0.000      0.000      crust  7883.104  4559.139  3267.717
 +
-119.1052    35.3449      0.000    114.266    280.000      cs173  3371.219  1940.716  2335.530      none      0.000      0.000      0.000      crust  3371.219  1940.716  2335.530
 +
-119.1052    35.3449    175.000    114.266    280.000      cs173  3382.928  1945.926  2336.248      none      0.000      0.000      0.000      crust  3382.928  1945.926  2336.248
 +
-119.1052    35.3449  1050.000    114.266    280.000      cs173  3516.983  2000.467  2343.681      none      0.000      0.000      0.000      crust  3516.983  2000.467  2343.681
 +
-119.1052    35.3449  9975.000    114.266    280.000      cs173  6906.483  4185.805  3063.830      none      0.000      0.000      0.000      crust  6906.483  4185.805  3063.830
 +
-119.1052    35.3449  40075.000    114.266    280.000      cs173  7823.070  4576.244  3276.429      none      0.000      0.000      0.000      crust  7823.070  4576.244  3276.429
 +
-120.9911    35.6692      0.000    612.945    782.000      cs173  4317.940  2337.993  2390.002      none      0.000      0.000      0.000      crust  4317.940  2337.993  2390.002
 +
-120.9911    35.6692    175.000    612.945    782.000      cs173  4330.583  2347.756  2391.460      none      0.000      0.000      0.000      crust  4330.583  2347.756  2391.460
 +
-120.9911    35.6692  1050.000    612.945    782.000      cs173  4516.919  2477.496  2411.950      none      0.000      0.000      0.000      crust  4516.919  2477.496  2411.950
 +
-120.9911    35.6692  9975.000    612.945    782.000      cs173  5711.902  3079.826  2552.101      none      0.000      0.000      0.000      crust  5711.902  3079.826  2552.101
 +
-120.9911    35.6692  40075.000    612.945    782.000      cs173  7870.165  4546.529  3261.238      none      0.000      0.000      0.000      crust  7870.165  4546.529  3261.238
 +
-120.4338    35.8993      0.000    468.744    349.000      cs173  4401.550  2376.688  2395.809      none      0.000      0.000      0.000      crust  4401.550  2376.688  2395.809
 +
-120.4338    35.8993    175.000    468.744    349.000      cs173  4403.934  2378.959  2396.156      none      0.000      0.000      0.000      crust  4403.934  2378.959  2396.156
 +
-120.4338    35.8993  1050.000    468.744    349.000      cs173  4379.906  2371.956  2395.088      none      0.000      0.000      0.000      crust  4379.906  2371.956  2395.088
 +
-120.4338    35.8993  9975.000    468.744    349.000      cs173  5789.239  3302.620  2629.347      none      0.000      0.000      0.000      crust  5789.239  3302.620  2629.347
 +
-120.4338    35.8993  40075.000    468.744    349.000      cs173  7586.570  4297.113  3126.202      none      0.000      0.000      0.000      crust  7586.570  4297.113  3126.202
 +
-120.5044    35.8108      0.000    553.646    455.000      cs173  5124.013  2859.918  2490.301      none      0.000      0.000      0.000      crust  5124.013  2859.918  2490.301
 +
-120.5044    35.8108    175.000    553.646    455.000      cs173  5115.474  2852.077  2488.361      none      0.000      0.000      0.000      crust  5115.474  2852.077  2488.361
 +
-120.5044    35.8108  1050.000    553.646    455.000      cs173  5039.287  2789.464  2473.323      none      0.000      0.000      0.000      crust  5039.287  2789.464  2473.323
 +
-120.5044    35.8108  9975.000    553.646    455.000      cs173  5692.855  3347.414  2646.662      none      0.000      0.000      0.000      crust  5692.855  3347.414  2646.662
 +
-120.5044    35.8108  40075.000    553.646    455.000      cs173  7612.118  4317.427  3137.526      none      0.000      0.000      0.000      crust  7612.118  4317.427  3137.526
 +
-119.7827    36.3016      0.000    69.588    280.000      cs173  3032.590  1490.514  2258.628      none      0.000      0.000      0.000      crust  3032.590  1490.514  2258.628
 +
-119.7827    36.3016    175.000    69.588    280.000      cs173  3048.569  1498.927  2260.458      none      0.000      0.000      0.000      crust  3048.569  1498.927  2260.458
 +
-119.7827    36.3016  1050.000    69.588    280.000      cs173  3273.004  1613.803  2283.640      none      0.000      0.000      0.000      crust  3273.004  1613.803  2283.640
 +
-119.7827    36.3016  9975.000    69.588    280.000      cs173  5843.598  3296.038  2626.953      none      0.000      0.000      0.000      crust  5843.598  3296.038  2626.953
 +
-119.7827    36.3016  40075.000    69.588    280.000      cs173  7657.531  4359.616  3160.929      none      0.000      0.000      0.000      crust  7657.531  4359.616  3160.929
 +
</pre>
 +
 +
 +
 +
== Comparison Plots CS173 ==
 +
 +
=== Raw mesh plots ===
 +
 +
Below are cross-section plots produced by using the CyberShake mesh plotting tools on the raw cs173 velocity mesh files.  The color scale is basemap.cm.GMT_seis, which was originally defined in GMT and then ported to basemap.
 +
 +
==== Horizontal cross-sections ====
 +
 +
The surface cross-section includes lines indicating where the vertical cross-section plots are.
 +
 +
{|
 +
| [[File:cs173_horiz_slice_0z.png|thumb|400px|z=0, depth=0 km, with vertical cross-section locations indicated]]
 +
| [[File:cs173_horiz_slice_6z.png|thumb|400px|z=6, depth=1.05 km]]
 +
| [[File:cs173_horiz_slice_57z.png|thumb|400px|z=57, depth=9.975 km]]
 +
| [[File:cs173_horiz_slice_114z.png|thumb|400px|z=114, depth=19.95 km]]
 +
|}
 +
 +
==== Vertical cross-sections ====
 +
 +
Parallel to Y-axis
 +
 +
X=1000: (32.550700, -115.409380) to (38.446980, -126.695380)
 +
 +
X=2000: (33.824270, -114.301790) to (39.813370, -125.716460)
 +
 +
{|
 +
| [[File:cs173_vert_slice_1000x.png|thumb|400px|X=1000 (blue line on horizontal plot)]]
 +
| [[File:cs173_vert_slice_2000x.png|thumb|400px|X=2000 (orange line on horizontal plot)]]
 +
|}
 +
 +
Parallel to X-axis
 +
 +
Y=1740: (32.824810, -119.141480) to (36.959180, -115.655030)
 +
 +
Y=3480: (34.320700, -121.888010) to (38.531550, -118.487980)
 +
 +
Y=5220: (35.744050, -124.726900) to (40.029510, -121.442660)
 +
 +
{|
 +
| [[File:cs173_vert_slice_1740y.png|thumb|400px|Y=1740 (green line on horizontal plot)]]
 +
| [[File:cs173_vert_slice_3480y.png|thumb|400px|Y=3480 (red line on horizontal plot)]]
 +
| [[File:cs173_vert_slice_5220y.png|thumb|400px|Y=5220 (purple line on horizontal plot)]]
 +
|}
 +
 +
=== UCVM plots ===
 +
 +
{|
 +
| [[File:scott_cs173_hort_0.png|thumb|400px|z=0, depth=0 km]]
 +
| [[File:scott_cs173_hort_1050.png|thumb|400px|z=6, depth=1.05 km]]
 +
|}
 +
{|
 +
| [[File:scott_cs173_hort_9975.png|thumb|400px|z=57, depth=9.975 km]]
 +
| [[File:scott_cs173_hort_19950.png|thumb|400px|z=114, depth=19.95 km]]
 +
|}
 +
{|
 +
| [[File:scott_cs173_hort_50400.png|thumb|400px|z=287, depth=50.4 km]]
 +
|}
 +
 +
Parallel to Y-axis
 +
 +
{|
 +
| [[File:scott_cs173_cross_x1.png|thumb|400px|X=1000]]
 +
| [[File:scott_cs173_cross_x2.png|thumb|400px|X=2000]]
 +
|}
 +
 +
Parallel to X-axis
 +
 +
{|
 +
| [[File:scott_cs173_cross_y1.png|thumb|400px|Y=1740]]
 +
| [[File:scott_cs173_cross_y2.png|thumb|400px|Y=3480]]
 +
|}
 +
{|
 +
| [[File:scott_cs173_cross_y3.png|thumb|400px|Y=5220]]
 +
|}
 +
 +
== Example Comparison Data Points ==
 +
* CS173
 +
* UCVM CS173
 +
 +
Vs at top 6 layers(counting down at steps of 175m) at lon/lats,(36.02718,-119.446), that translates to grid point (1600, 3000) on the mesh
 +
 +
    x 1600 y 3000 z 287  >>  vs 1776.336060
 +
    x 1600 y 3000 z 286  >>  vs 1789.887817
 +
    x 1600 y 3000 z 285  >>  vs 1803.439453
 +
    x 1600 y 3000 z 284  >>  vs 1818.368774
 +
    x 1600 y 3000 z 283  >>  vs 1841.563721
 +
    x 1600 y 3000 z 282  >>  vs 1864.758545
  
 
== Velocity mesh region ==
 
== Velocity mesh region ==
  
−
The region we will extract is 560 km x 1218 km x 50.4 km, with 55 degree rotation counter-clockwise.  The mesh will have 175m grid point spacing. The corner points are approximately
+
The region we will extract is 560 km x 1218 km x 50.4 km (as measured in the WGS84 coordinate system), with 55 degree rotation counter-clockwise.  The mesh will have 175m grid point spacing.
−
 
 
−
W: (37.30489, -127.62831)
 
−
N: (41.59023, -124.28641)
 
−
E: (35.17113, -112.97233)
 
−
S: (31.19740, -116.67400)
 
  
 
== Mesh extraction process ==
 
== Mesh extraction process ==
Line 15: Line 321:
  
 
#CyberShake processing tools will be used to construct a mesh, with the location of each point defined as km offsets from the center of the mesh at the surface.
 
#CyberShake processing tools will be used to construct a mesh, with the location of each point defined as km offsets from the center of the mesh at the surface.
−
## The mesh will evenly discretize the region in x,y,z, but do to the rotation, this is not North,South. The resolution will be 175km as measured in the UTM coordinate system.
+
## The mesh will evenly discretize the region in x,y,z, but due to the rotation, this is not North,South. The resolution will be 175m as measured in the UTM coordinate system.
−
##A number of CyberShake codes, originally from RWG, are involved in this step. These codes use his great circle projection.  These codes figure out the outer bounds of the mesh given the location of a site and fault surfaces. They use great circle will be a problem, and return a lat/lon of the mesh center.  
+
##A number of CyberShake codes, originally from RWG, are involved in this step. These codes use his great circle projection.  These codes figure out the outer bounds of the mesh given the location of a site and fault surfaces. They return a lat/lon of the mesh center.  
−
## The codes used are in version control here: [CyberShake Projection Codes]
+
## The codes used are in version control here: [http://source.usc.edu/svn/cybershake/import/trunk/PreCVM].
−
##Then we use CyberShake codes convert this center point to UTM Zone 11, which is the southern California UTM zone. While some of the mesh region extends north into UTM zone 10, we will use Zone 11 since the majority of our region of interest is in that zone.  We'll do the remaining next steps in UTM.   
+
##Then we use the Proj 4 library to convert this center point to UTM Zone 11, which is the southern California UTM zone. While some of the mesh region extends north into UTM zone 10, we will use Zone 11 since the majority of our region of interest is in that zone.  We'll do the remaining next steps in UTM.   
 +
#The location of each mesh point will be redefined as a UTM easting and northing, by translating the center of the mesh from lat/lon to UTM coordinates in Zone 11 and then applying the offsets.
 +
##The UTM mesh points will be written to a file, in fast y, x, z format, starting with all the mesh points are the surface, first changing y values, then x values, then moving down one the Z axis, and listing the next layer of points.  These points will be in the same order as the points in the velocity mesh.
 +
##This file in internal, and will not be part of the data sent to the USR group.
 +
#The pyproj package will be used to convert the UTM coordinates at each point into a latitude and longitude value, to 5 decimal places.
 +
##The proj4 coordinate system to coordinate system (cs2cs) conversion string used will be proj='utm', zone=11, ellps='WGS84'.
 +
#This list of points will be used to query UCVM v15.10 through the C API, using CyberShake's [[CyberShake_Code_Base#UCVM|UCVM wrapper]], same as used in CyberShake Study 17.3, and create a velocity mesh.  For a given point, we will first query the CCA-06 model, then CVM-S4.26, then USGS Bay Area, and finally a background 1D model.  Additionally, we will query the USGS topo model in UCVM 17.3 to obtain the surface elevation at each point and write it to a file.
 +
#Smoothing will be applied to the velocity mesh for all points within 10 km of an interface, by averaging the neighbors for 10km in the N, S, E, and W directions.  We will use the [[CyberShake_Code_Base#Smoothing|CyberShake smoothing code]].  Sample plots illustrating the effect of this smoothing are available [[Study_16.9_Velocity_Plots|here]].
 +
#Scott will deliver three meshes (one each with Vp, Vs, and rho) in fast y, x, z binary format with 4-byte floats, a list of the elevation of each surface point, a list of the coordinates of each surface point in lat/lon and UTM zone 11, and MD5 sums for all products.
 +
 
 +
== Mesh region details ==
 +
 
 +
The actual mesh we extracted has its center at (36.44287, -120.22745), or 210707 easting, 4037914 northing in UTM Zone 11.
 +
 
 +
The actual corner points are:
 +
S: (31.26643, -116.48562)
 +
W: (37.08416, -127.65648)
 +
E: (35.33518, -112.92896)
 +
N: (41.45284, -124.51032)
 +
 
 +
Using proj4 project, we project these points into UTM Zone 10 (since most of the mesh is in northern California:
  
 +
Lat/Lon WGS84 to UTM Zone 10
 +
S: (31.26643, -116.48562) : 1120786.790041  3477495.533344
 +
W: (37.08416, -127.65648) : 86016.577605    4114365.092550
 +
E: (35.33518, -112.92896)  : 1416797.228564  3957099.088864
 +
N: (41.45284, -124.51032) : 373851.268121 4590129.392433
  
  
−
#The location of each mesh point will be redefined as a UTM easting and northing, by translating the center of the mesh from lat/lon to UTM coordinates in Zone 11 and then applying the offsets.
+
Lat/Lon WGS84 to UTM Zone 11
 +
S: (31.26643, -116.48562) : 548969.079292 3459243.769232
 +
W: (37.08416, -127.65648) : -448610.671359 4157755.975881
 +
E: (35.33518, -112.92896) : 870071.761707 3917825.705794
 +
N: (41.45284, -124.51032) :  -127508.133290 4616338.305196
 +
 
 +
== Content of the config file ==
  
 +
  # UTM Zone
 +
  utm_zone = 11
 +
  #
 +
  # Model directory
 +
  model_dir = cs173
 +
  #
 +
  # Number of cells in x, y, and z.
 +
  nx = 3200
 +
  ny = 6960
 +
  nz = 288
 +
  #
 +
  # Depth, in meters
 +
  depth = 50400
 +
  depth_interval = 175
 +
  #
 +
  # Four corners
 +
  bottom_left_corner_e = 548969.079292
 +
  bottom_left_corner_n = 3459243.769232
 +
  top_left_corner_e    = -448610.671359
 +
  top_left_corner_n    = 4157755.975881
 +
  bottom_right_corner_e = 870071.761707
 +
  bottom_right_corner_n = 3917825.705794
 +
  top_right_corner_e  = -127508.133290
 +
  top_right_corner_n  = 4616338.305196
 +
  #
 +
  # Seek method, fast-X or fast-Y
 +
  seek_axis = fast-Y
 +
  seek_direction = top-down
  
−
#The Proj4 package will be used to convert the UTM coordinates at each point into a latitude and longitude value, to 5 decimal places.
+
== Mesh integration process ==
−
#This list of points will be used to query UCVM and create a velocity mesh.  Additionally, we will query the USGS topo model in UCVM to obtain the surface elevation at each point.
+
# Returned mesh will be registered into UCVM
−
#Smoothing will be applied to the velocity mesh for all points within 10 km of an interface, by averaging the neighbors for 10km in the N, S, E, and W directions.
+
#We will develop evaluation plan to determine whether returned points are registered properly.
−
#Scott will deliver three meshes (one each with Vp, Vs, and rho) in fast y, x, z binary format with 4-byte floats, a list of the elevation of each surface point, a list of the coordinates of each surface point in lat/lon and UTM zone 11, and MD5 sums for all products.
 
  
 
== Model Coverage Regions ==
 
== Model Coverage Regions ==
Line 34: Line 398:
  
 
== Related Entries ==
 
== Related Entries ==
−
*[USR]
+
 
−
*[CVM-S]
+
*[[USR]]
−
*[CVM-H]
+
*[[CVM-S]]
−
*[UCVM]
+
*[[CVM-H]]
 +
*[[UCVM]]
 +
*[[Proj4]]

Latest revision as of 21:48, 15 May 2018

This page details the steps involved in delivering a velocity mesh to John Shaw and Andreas Plesch, so that the Harvard group's Central California basins can be integrated into a velocity model, which can be fed back into UCVM and used for CyberShake.

Test Points

To evaluate the registration of CS17.3 into UCVM we are querying a set of point using Scott's original CS17.3 and the recently updated UCVMC version. The following test points are sites from the CS17.3 study. We will query each point at 5 depths 0km, 100m, 1km, 10km, and 40km 

<name>Osito Audit: Castaic Lake Dam (OSI)</name>
<coordinates>-118.723500,34.614500</coordinates>

<name>Santa Barbara (SBR)</name>
<coordinates>-119.698189,34.420831</coordinates>

<name>s1356 (s1356)</name>
<coordinates>-120.582740,34.968060</coordinates>

<name>Calstate Bakersfield (BAK)</name>	
<coordinates>-119.104450,35.344440</coordinates>

<name>Lime Mtn. Quarry (LMQ)</name>
<coordinates>-120.990669,35.668835</coordinates>

<name>Parkfield (PARK)</name>
<coordinates>-120.432800,35.899700</coordinates>

<name>Huron (HUR)</name>
<coordinates>-120.102917,36.202731</coordinates>

<name>Lemoore (LEM)</name>
<coordinates>-119.782911,36.300783</coordinates>

Test points and depths list

-118.723500 34.614500 0.0
-118.723500 34.614500 100.0
-118.723500 34.614500 1000.0
-118.723500 34.614500 10000.0
-118.723500 34.614500 40000.0
-119.698189 34.420831 0.0
-119.698189 34.420831 100.0
-119.698189 34.420831 1000.0
-119.698189 34.420831 10000.0
-119.698189 34.420831 40000.0
-120.582740 34.968060 0.0
-120.582740 34.968060 100.0
-120.582740 34.968060 1000.0
-120.582740 34.968060 10000.0
-120.582740 34.968060 40000.0
-119.104450 35.344440 0.0
-119.104450 35.344440 100.0
-119.104450 35.344440 1000.0
-119.104450 35.344440 10000.0
-119.104450 35.344440 40000.0
-120.990669 35.668835 0.0
-120.990669 35.668835 100.0
-120.990669 35.668835 1000.0
-120.990669 35.668835 10000.0
-120.990669 35.668835 40000.0
-120.432800 35.899700 0.0
-120.432800 35.899700 100.0
-120.432800 35.899700 1000.0
-120.432800 35.899700 10000.0
-120.432800 35.899700 40000.0
-120.102917 36.202731 0.0
-120.102917 36.202731 100.0
-120.102917 36.202731 1000.0
-120.102917 36.202731 10000.0
-120.102917 36.202731 40000.0
-119.782911 36.300783 0.0
-119.782911 36.300783 100.0
-119.782911 36.300783 1000.0
-119.782911 36.300783 10000.0
-119.782911 36.300783 40000.0


ucvm_query result from UCVM-C 18.5 alpha 

../bin/ucvm_query -f ../conf/ucvm.conf -m cs173 << EOF &> cs173-it.log
-118.723500 34.614500 0.0
...
EOF

WARNING: Could not load model into memory. Reading the model from the
hard disk may result in slow performance.
Using Geo Depth coordinates as default mode.
 -118.7235    34.6145      0.000    713.454    390.000      cs173   2047.660   1089.404   2139.522       none      0.000      0.000      0.000      crust   2047.660   1089.404   2139.522
 -118.7235    34.6145    100.000    713.454    390.000      cs173   2734.936   1518.460   2232.448       none      0.000      0.000      0.000      crust   2734.936   1518.460   2232.448
 -118.7235    34.6145   1000.000    713.454    390.000      cs173   5460.268   3237.855   2605.311       none      0.000      0.000      0.000      crust   5460.268   3237.855   2605.311
 -118.7235    34.6145  10000.000    713.454    390.000      cs173   6561.615   3871.144   2890.995       none      0.000      0.000      0.000      crust   6561.615   3871.144   2890.995
 -118.7235    34.6145  40000.000    713.454    390.000      cs173   7545.088   4324.382   3141.390       none      0.000      0.000      0.000      crust   7545.088   4324.382   3141.390
 -119.6982    34.4208      0.000     15.086    387.000      cs173   4670.389   2585.928   2430.952       none      0.000      0.000      0.000      crust   4670.389   2585.928   2430.952
 -119.6982    34.4208    100.000     15.086    387.000      cs173   4668.188   2583.502   2430.503       none      0.000      0.000      0.000      crust   4668.188   2583.502   2430.503
 -119.6982    34.4208   1000.000     15.086    387.000      cs173   4625.225   2526.679   2420.277       none      0.000      0.000      0.000      crust   4625.225   2526.679   2420.277
 -119.6982    34.4208  10000.000     15.086    387.000      cs173   5915.015   3283.019   2621.902       none      0.000      0.000      0.000      crust   5915.015   3283.019   2621.902
 -119.6982    34.4208  40000.000     15.086    387.000      cs173   7619.421   4327.448   3143.100       none      0.000      0.000      0.000      crust   7619.421   4327.448   3143.100
 -120.5827    34.9681      0.000     18.443    280.000      cs173   4767.372   2397.728   2399.038       none      0.000      0.000      0.000      crust   4767.372   2397.728   2399.038
 -120.5827    34.9681    100.000     18.443    280.000      cs173   4771.665   2403.488   2399.941       none      0.000      0.000      0.000      crust   4771.665   2403.488   2399.941
 -120.5827    34.9681   1000.000     18.443    280.000      cs173   4833.515   2478.506   2412.065       none      0.000      0.000      0.000      crust   4833.515   2478.506   2412.065
 -120.5827    34.9681  10000.000     18.443    280.000      cs173   5818.455   3107.401   2560.882       none      0.000      0.000      0.000      crust   5818.455   3107.401   2560.882
 -120.5827    34.9681  40000.000     18.443    280.000      cs173   7883.487   4559.606   3267.956       none      0.000      0.000      0.000      crust   7883.487   4559.606   3267.956
 -119.1044    35.3444      0.000    114.259    280.000      cs173   3372.060   1940.204   2335.459       none      0.000      0.000      0.000      crust   3372.060   1940.204   2335.459
 -119.1044    35.3444    100.000    114.259    280.000      cs173   3378.782   1943.201   2335.872       none      0.000      0.000      0.000      crust   3378.782   1943.201   2335.872
 -119.1044    35.3444   1000.000    114.259    280.000      cs173   3507.057   1995.574   2343.018       none      0.000      0.000      0.000      crust   3507.057   1995.574   2343.018
 -119.1044    35.3444  10000.000    114.259    280.000      cs173   6903.631   4181.306   3061.306       none      0.000      0.000      0.000      crust   6903.631   4181.306   3061.306
 -119.1044    35.3444  40000.000    114.259    280.000      cs173   7822.489   4575.821   3276.215       none      0.000      0.000      0.000      crust   7822.489   4575.821   3276.215
 -120.9907    35.6688      0.000    602.291    782.000      cs173   4317.442   2337.356   2389.908       none      0.000      0.000      0.000      crust   4317.442   2337.356   2389.908
 -120.9907    35.6688    100.000    602.291    782.000      cs173   4324.663   2342.934   2390.741       none      0.000      0.000      0.000      crust   4324.663   2342.934   2390.741
 -120.9907    35.6688   1000.000    602.291    782.000      cs173   4501.126   2466.457   2410.151       none      0.000      0.000      0.000      crust   4501.126   2466.457   2410.151
 -120.9907    35.6688  10000.000    602.291    782.000      cs173   5709.418   3078.566   2551.708       none      0.000      0.000      0.000      crust   5709.418   3078.566   2551.708
 -120.9907    35.6688  40000.000    602.291    782.000      cs173   7870.905   4547.270   3261.620       none      0.000      0.000      0.000      crust   7870.905   4547.270   3261.620
 -120.4328    35.8997      0.000    467.078    349.000      cs173   4388.868   2370.534   2394.871       none      0.000      0.000      0.000      crust   4388.868   2370.534   2394.871
 -120.4328    35.8997    100.000    467.078    349.000      cs173   4390.334   2371.916   2395.081       none      0.000      0.000      0.000      crust   4390.334   2371.916   2395.081
 -120.4328    35.8997   1000.000    467.078    349.000      cs173   4371.586   2367.767   2394.453       none      0.000      0.000      0.000      crust   4371.586   2367.767   2394.453
 -120.4328    35.8997  10000.000    467.078    349.000      cs173   5795.831   3301.412   2628.905       none      0.000      0.000      0.000      crust   5795.831   3301.412   2628.905
 -120.4328    35.8997  40000.000    467.078    349.000      cs173   7586.191   4296.916   3126.092       none      0.000      0.000      0.000      crust   7586.191   4296.916   3126.092
 -120.1029    36.2027      0.000    113.530    280.000      cs173   2851.014   1495.107   2259.639       none      0.000      0.000      0.000      crust   2851.014   1495.107   2259.639
 -120.1029    36.2027    100.000    113.530    280.000      cs173   2859.429   1498.931   2260.468       none      0.000      0.000      0.000      crust   2859.429   1498.931   2260.468
 -120.1029    36.2027   1000.000    113.530    280.000      cs173   3050.430   1580.739   2277.315       none      0.000      0.000      0.000      crust   3050.430   1580.739   2277.315
 -120.1029    36.2027  10000.000    113.530    280.000      cs173   5612.277   3107.839   2561.259       none      0.000      0.000      0.000      crust   5612.277   3107.839   2561.259
 -120.1029    36.2027  40000.000    113.530    280.000      cs173   7696.131   4393.316   3179.477       none      0.000      0.000      0.000      crust   7696.131   4393.316   3179.477
 -119.7829    36.3008      0.000     69.541    280.000      cs173   3031.683   1489.972   2258.508       none      0.000      0.000      0.000      crust   3031.683   1489.972   2258.508
 -119.7829    36.3008    100.000     69.541    280.000      cs173   3040.815   1494.782   2259.555       none      0.000      0.000      0.000      crust   3040.815   1494.782   2259.555
 -119.7829    36.3008   1000.000     69.541    280.000      cs173   3253.721   1604.038   2281.744       none      0.000      0.000      0.000      crust   3253.721   1604.038   2281.744
 -119.7829    36.3008  10000.000     69.541    280.000      cs173   5852.581   3303.596   2629.943       none      0.000      0.000      0.000      crust   5852.581   3303.596   2629.943
 -119.7829    36.3008  40000.000     69.541    280.000      cs173   7657.497   4359.726   3160.989       none      0.000      0.000      0.000      crust   7657.497   4359.726   3160.989

Another set of Test Points

-118.72254 34.61444 0
-118.72254 34.61444 175
-118.72254 34.61444 1050
-118.72254 34.61444 9975
-118.72254 34.61444 40075
-119.69846 34.42187 0
-119.69846 34.42187 175
-119.69846 34.42187 1050
-119.69846 34.42187 9975
-119.69846 34.42187 40075
-120.58229 34.96759 0
-120.58229 34.96759 175
-120.58229 34.96759 1050
-120.58229 34.96759 9975
-120.58229 34.96759 40075
-119.10517 35.34494 0
-119.10517 35.34494 175
-119.10517 35.34494 1050
-119.10517 35.34494 9975
-119.10517 35.34494 40075
-120.99114 35.66917 0
-120.99114 35.66917 175
-120.99114 35.66917 1050
-120.99114 35.66917 9975
-120.99114 35.66917 40075
-120.43380 35.89934 0
-120.43380 35.89934 175
-120.43380 35.89934 1050
-120.43380 35.89934 9975
-120.43380 35.89934 40075
-120.50437 35.81076 0
-120.50437 35.81076 175
-120.50437 35.81076 1050
-120.50437 35.81076 9975
-120.50437 35.81076 40075
-119.78270 36.30157 0
-119.78270 36.30157 175
-119.78270 36.30157 1050
-119.78270 36.30157 9975
-119.78270 36.30157 40075


ucvm_query result 

-118.7225    34.6144      0.000    720.976    390.000      cs173   2043.192   1085.995   2138.201       none      0.000      0.000      0.000      crust   2043.192   1085.995   2138.201
 -118.7225    34.6144    175.000    720.976    390.000      cs173   3248.173   1837.993   2301.267       none      0.000      0.000      0.000      crust   3248.173   1837.993   2301.267
 -118.7225    34.6144   1050.000    720.976    390.000      cs173   5459.535   3237.474   2605.180       none      0.000      0.000      0.000      crust   5459.535   3237.474   2605.180
 -118.7225    34.6144   9975.000    720.976    390.000      cs173   6556.426   3867.625   2889.135       none      0.000      0.000      0.000      crust   6556.426   3867.625   2889.135
 -118.7225    34.6144  40075.000    720.976    390.000      cs173   7551.173   4330.115   3144.578       none      0.000      0.000      0.000      crust   7551.173   4330.115   3144.578
 -119.6985    34.4219      0.000     16.241    387.000      cs173   4656.098   2577.254   2429.342       none      0.000      0.000      0.000      crust   4656.098   2577.254   2429.342
 -119.6985    34.4219    175.000     16.241    387.000      cs173   4652.338   2573.052   2428.572       none      0.000      0.000      0.000      crust   4652.338   2573.052   2428.572
 -119.6985    34.4219   1050.000     16.241    387.000      cs173   4609.664   2514.169   2418.086       none      0.000      0.000      0.000      crust   4609.664   2514.169   2418.086
 -119.6985    34.4219   9975.000     16.241    387.000      cs173   5913.173   3281.206   2621.216       none      0.000      0.000      0.000      crust   5913.173   3281.206   2621.216
 -119.6985    34.4219  40075.000     16.241    387.000      cs173   7622.714   4330.343   3144.709       none      0.000      0.000      0.000      crust   7622.714   4330.343   3144.709
 -120.5823    34.9676      0.000     18.570    280.000      cs173   4765.085   2395.066   2398.626       none      0.000      0.000      0.000      crust   4765.085   2395.066   2398.626
 -120.5823    34.9676    175.000     18.570    280.000      cs173   4772.822   2405.407   2400.245       none      0.000      0.000      0.000      crust   4772.822   2405.407   2400.245
 -120.5823    34.9676   1050.000     18.570    280.000      cs173   4837.920   2483.522   2412.900       none      0.000      0.000      0.000      crust   4837.920   2483.522   2412.900
 -120.5823    34.9676   9975.000     18.570    280.000      cs173   5814.321   3105.312   2560.201       none      0.000      0.000      0.000      crust   5814.321   3105.312   2560.201
 -120.5823    34.9676  40075.000     18.570    280.000      cs173   7883.104   4559.139   3267.717       none      0.000      0.000      0.000      crust   7883.104   4559.139   3267.717
 -119.1052    35.3449      0.000    114.266    280.000      cs173   3371.219   1940.716   2335.530       none      0.000      0.000      0.000      crust   3371.219   1940.716   2335.530
 -119.1052    35.3449    175.000    114.266    280.000      cs173   3382.928   1945.926   2336.248       none      0.000      0.000      0.000      crust   3382.928   1945.926   2336.248
-119.1052    35.3449   1050.000    114.266    280.000      cs173   3516.983   2000.467   2343.681       none      0.000      0.000      0.000      crust   3516.983   2000.467   2343.681
 -119.1052    35.3449   9975.000    114.266    280.000      cs173   6906.483   4185.805   3063.830       none      0.000      0.000      0.000      crust   6906.483   4185.805   3063.830
 -119.1052    35.3449  40075.000    114.266    280.000      cs173   7823.070   4576.244   3276.429       none      0.000      0.000      0.000      crust   7823.070   4576.244   3276.429
 -120.9911    35.6692      0.000    612.945    782.000      cs173   4317.940   2337.993   2390.002       none      0.000      0.000      0.000      crust   4317.940   2337.993   2390.002
 -120.9911    35.6692    175.000    612.945    782.000      cs173   4330.583   2347.756   2391.460       none      0.000      0.000      0.000      crust   4330.583   2347.756   2391.460
 -120.9911    35.6692   1050.000    612.945    782.000      cs173   4516.919   2477.496   2411.950       none      0.000      0.000      0.000      crust   4516.919   2477.496   2411.950
 -120.9911    35.6692   9975.000    612.945    782.000      cs173   5711.902   3079.826   2552.101       none      0.000      0.000      0.000      crust   5711.902   3079.826   2552.101
 -120.9911    35.6692  40075.000    612.945    782.000      cs173   7870.165   4546.529   3261.238       none      0.000      0.000      0.000      crust   7870.165   4546.529   3261.238
 -120.4338    35.8993      0.000    468.744    349.000      cs173   4401.550   2376.688   2395.809       none      0.000      0.000      0.000      crust   4401.550   2376.688   2395.809
 -120.4338    35.8993    175.000    468.744    349.000      cs173   4403.934   2378.959   2396.156       none      0.000      0.000      0.000      crust   4403.934   2378.959   2396.156
 -120.4338    35.8993   1050.000    468.744    349.000      cs173   4379.906   2371.956   2395.088       none      0.000      0.000      0.000      crust   4379.906   2371.956   2395.088
 -120.4338    35.8993   9975.000    468.744    349.000      cs173   5789.239   3302.620   2629.347       none      0.000      0.000      0.000      crust   5789.239   3302.620   2629.347
 -120.4338    35.8993  40075.000    468.744    349.000      cs173   7586.570   4297.113   3126.202       none      0.000      0.000      0.000      crust   7586.570   4297.113   3126.202
 -120.5044    35.8108      0.000    553.646    455.000      cs173   5124.013   2859.918   2490.301       none      0.000      0.000      0.000      crust   5124.013   2859.918   2490.301
 -120.5044    35.8108    175.000    553.646    455.000      cs173   5115.474   2852.077   2488.361       none      0.000      0.000      0.000      crust   5115.474   2852.077   2488.361
 -120.5044    35.8108   1050.000    553.646    455.000      cs173   5039.287   2789.464   2473.323       none      0.000      0.000      0.000      crust   5039.287   2789.464   2473.323
 -120.5044    35.8108   9975.000    553.646    455.000      cs173   5692.855   3347.414   2646.662       none      0.000      0.000      0.000      crust   5692.855   3347.414   2646.662
-120.5044    35.8108  40075.000    553.646    455.000      cs173   7612.118   4317.427   3137.526       none      0.000      0.000      0.000      crust   7612.118   4317.427   3137.526
 -119.7827    36.3016      0.000     69.588    280.000      cs173   3032.590   1490.514   2258.628       none      0.000      0.000      0.000      crust   3032.590   1490.514   2258.628
 -119.7827    36.3016    175.000     69.588    280.000      cs173   3048.569   1498.927   2260.458       none      0.000      0.000      0.000      crust   3048.569   1498.927   2260.458
 -119.7827    36.3016   1050.000     69.588    280.000      cs173   3273.004   1613.803   2283.640       none      0.000      0.000      0.000      crust   3273.004   1613.803   2283.640
 -119.7827    36.3016   9975.000     69.588    280.000      cs173   5843.598   3296.038   2626.953       none      0.000      0.000      0.000      crust   5843.598   3296.038   2626.953
 -119.7827    36.3016  40075.000     69.588    280.000      cs173   7657.531   4359.616   3160.929       none      0.000      0.000      0.000      crust   7657.531   4359.616   3160.929


Comparison Plots CS173

Raw mesh plots

Below are cross-section plots produced by using the CyberShake mesh plotting tools on the raw cs173 velocity mesh files. The color scale is basemap.cm.GMT_seis, which was originally defined in GMT and then ported to basemap.

Horizontal cross-sections

The surface cross-section includes lines indicating where the vertical cross-section plots are.

z=0, depth=0 km, with vertical cross-section locations indicated
z=6, depth=1.05 km
z=57, depth=9.975 km
z=114, depth=19.95 km

Vertical cross-sections

Parallel to Y-axis

X=1000: (32.550700, -115.409380) to (38.446980, -126.695380)

X=2000: (33.824270, -114.301790) to (39.813370, -125.716460)

X=1000 (blue line on horizontal plot)
X=2000 (orange line on horizontal plot)

Parallel to X-axis

Y=1740: (32.824810, -119.141480) to (36.959180, -115.655030)

Y=3480: (34.320700, -121.888010) to (38.531550, -118.487980)

Y=5220: (35.744050, -124.726900) to (40.029510, -121.442660)

Y=1740 (green line on horizontal plot)
Y=3480 (red line on horizontal plot)
Y=5220 (purple line on horizontal plot)

UCVM plots

z=0, depth=0 km
z=6, depth=1.05 km
z=57, depth=9.975 km
z=114, depth=19.95 km
z=287, depth=50.4 km

Parallel to Y-axis

X=1000
X=2000

Parallel to X-axis

Y=1740
Y=3480
Y=5220

Example Comparison Data Points

  • CS173
  • UCVM CS173

Vs at top 6 layers(counting down at steps of 175m) at lon/lats,(36.02718,-119.446), that translates to grid point (1600, 3000) on the mesh 

   x 1600 y 3000 z 287  >>  vs 1776.336060
   x 1600 y 3000 z 286  >>  vs 1789.887817
   x 1600 y 3000 z 285  >>  vs 1803.439453
   x 1600 y 3000 z 284  >>  vs 1818.368774
   x 1600 y 3000 z 283  >>  vs 1841.563721
   x 1600 y 3000 z 282  >>  vs 1864.758545

Velocity mesh region

The region we will extract is 560 km x 1218 km x 50.4 km (as measured in the WGS84 coordinate system), with 55 degree rotation counter-clockwise. The mesh will have 175m grid point spacing.

Mesh extraction process

To deliver the mesh, we will follow the following procedure:

  1. CyberShake processing tools will be used to construct a mesh, with the location of each point defined as km offsets from the center of the mesh at the surface.
    1. The mesh will evenly discretize the region in x,y,z, but due to the rotation, this is not North,South. The resolution will be 175m as measured in the UTM coordinate system.
    2. A number of CyberShake codes, originally from RWG, are involved in this step. These codes use his great circle projection. These codes figure out the outer bounds of the mesh given the location of a site and fault surfaces. They return a lat/lon of the mesh center.
    3. The codes used are in version control here: [1].
    4. Then we use the Proj 4 library to convert this center point to UTM Zone 11, which is the southern California UTM zone. While some of the mesh region extends north into UTM zone 10, we will use Zone 11 since the majority of our region of interest is in that zone. We'll do the remaining next steps in UTM.
  2. The location of each mesh point will be redefined as a UTM easting and northing, by translating the center of the mesh from lat/lon to UTM coordinates in Zone 11 and then applying the offsets.
    1. The UTM mesh points will be written to a file, in fast y, x, z format, starting with all the mesh points are the surface, first changing y values, then x values, then moving down one the Z axis, and listing the next layer of points. These points will be in the same order as the points in the velocity mesh.
    2. This file in internal, and will not be part of the data sent to the USR group.
  3. The pyproj package will be used to convert the UTM coordinates at each point into a latitude and longitude value, to 5 decimal places.
    1. The proj4 coordinate system to coordinate system (cs2cs) conversion string used will be proj='utm', zone=11, ellps='WGS84'.
  4. This list of points will be used to query UCVM v15.10 through the C API, using CyberShake's UCVM wrapper, same as used in CyberShake Study 17.3, and create a velocity mesh. For a given point, we will first query the CCA-06 model, then CVM-S4.26, then USGS Bay Area, and finally a background 1D model. Additionally, we will query the USGS topo model in UCVM 17.3 to obtain the surface elevation at each point and write it to a file.
  5. Smoothing will be applied to the velocity mesh for all points within 10 km of an interface, by averaging the neighbors for 10km in the N, S, E, and W directions. We will use the CyberShake smoothing code. Sample plots illustrating the effect of this smoothing are available here.
  6. Scott will deliver three meshes (one each with Vp, Vs, and rho) in fast y, x, z binary format with 4-byte floats, a list of the elevation of each surface point, a list of the coordinates of each surface point in lat/lon and UTM zone 11, and MD5 sums for all products.

Mesh region details

The actual mesh we extracted has its center at (36.44287, -120.22745), or 210707 easting, 4037914 northing in UTM Zone 11.

The actual corner points are: 

S: (31.26643, -116.48562)
W: (37.08416, -127.65648)
E: (35.33518, -112.92896)
N: (41.45284, -124.51032)

Using proj4 project, we project these points into UTM Zone 10 (since most of the mesh is in northern California:

Lat/Lon WGS84 to UTM Zone 10 

S: (31.26643, -116.48562) : 1120786.790041  3477495.533344
W: (37.08416, -127.65648) : 86016.577605    4114365.092550
E: (35.33518, -112.92896)  : 1416797.228564  3957099.088864
N: (41.45284, -124.51032) : 373851.268121 4590129.392433


Lat/Lon WGS84 to UTM Zone 11 

S: (31.26643, -116.48562) : 548969.079292 3459243.769232
W: (37.08416, -127.65648) : -448610.671359 4157755.975881
E: (35.33518, -112.92896) : 870071.761707 3917825.705794
N: (41.45284, -124.51032) :  -127508.133290 4616338.305196

Content of the config file

 # UTM Zone
 utm_zone = 11
 #
 # Model directory
 model_dir = cs173
 #
 # Number of cells in x, y, and z.
 nx = 3200
 ny = 6960
 nz = 288
 #
 # Depth, in meters
 depth = 50400
 depth_interval = 175
 # 
 # Four corners
 bottom_left_corner_e = 548969.079292
 bottom_left_corner_n = 3459243.769232
 top_left_corner_e    = -448610.671359
 top_left_corner_n    = 4157755.975881
 bottom_right_corner_e = 870071.761707
 bottom_right_corner_n = 3917825.705794
 top_right_corner_e   = -127508.133290
 top_right_corner_n   = 4616338.305196
 #
 # Seek method, fast-X or fast-Y
 seek_axis = fast-Y
 seek_direction = top-down

Mesh integration process

  1. Returned mesh will be registered into UCVM
  2. We will develop evaluation plan to determine whether returned points are registered properly.

Model Coverage Regions

The mesh is depicted in white; the region covered by CCA-06 is in blue, by CVM-S4.26 in red, and by USGS Bay Area in green.

Related Entries

Retrieved from "https://strike.scec.org/scecwiki/index.php?title=Registering_CS173_into_UCVM&oldid=20626"