Basin Query Result

Basin Query Algorithms in Actual Implementation

For each grid point within the user supplied boundary, a stack of query points is generated with different depths in the increment of 20m. This set of query points are submitted to the backend model for their matching Vs.

Search of Vs threshold crossings of a grid point starts by looking at the Vs of this grid point from the first depth layer to the very last depth layer.

If the Vs is less than 0.0, it is invalid, just skip it.

If the Vs is greater than 0.0, check if this is a crossing point. A crossing is found if Vs of the earlier depth layer is lower than the preset threshold and the Vs of the current depth layer is higher than the preset threshold, ie. an upward sloping crossing.

5 different types of result for all the grid points are being collected and outputted into 5 different binary data files. User may selectively collect each or all data files by supplying the matching data file names on command line.

• First Crossing -- If there is no crossing, return -1 or else return the first upward crossing's depth

• Second or First Crossing -- If there is no crossing, return -1. If there are more than one crossing, return the second to the first or else return the depth of the first crossing

• Last Crossing -- If there is no crossing, return -1 or return the last upward crossing's depth

• Second only Crossing -- If there is no crossing or less than two crossings, return -1. Return the second crossing depth even if there are more than two crossings.

• Last but not Second or First Crossing -- If there is no crossing or there are less than or equal to two crossings, return -1, or return the last crossing's depth

Basin Query Algorithm Result for CCA 06

Algorithm	Threshold = 1000m/s	Threshold = 2500m/s
First Crossing	data png	data png
First Or Second Crossing	data png	data png
Last Crossing	data png	data png|
Second Only Crossing	data png	data png
Last but not First or Second Crossing	data png	data png

The basin query command that was used to generate the result,

CCA06: [minLat = 33.35, minLon = -123.0, nx = 1551, ny = 1201, gridSpacing = 0.005 degrees]

mpirun -np 20 ./basin_query_mpi -b ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.first, ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.firstOrSecond, ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.last, ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.secondOnly, ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.threeLast -o ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.result -f ucvm.conf -m cca -i 20 -v 2500 -l 33.35,-123.0 -s 0.005 -x 1551 -y 1201

A plot command that was used to generate one of the png file, (data file were renamed to simplify the batch script)

./plot_z25_map.py -b 33.35,-123.0 -u 39.35,-115.25 -c cca -s 0.005 -f cca_z2.5.first -o cca_z2.5_first_map.png

Basin Query Result

These result were ran on HPC and Blue Waters

Model	First Upward Crossing	Last Upward Crossing
CVM-S4.26 (cvms5)	threshold = 1000m/s threshold = 2500m/s	threshold = 1000m/s threshold = 2500m/s
CVM-S4.26-M01 (cvmsi)	threshold = 1000m/s threshold = 2500m/s	threshold = 1000m/s threshold = 2500m/s
CCA06 (cca)	threshold = 1000m/s threshold = 2500m/s	threshold = 1000m/s threshold = 2500m/s

Here are comparable results from OpenSHA. First crossing results were calculated by David Gill, either in UCVMC or python, I'm not sure (but 4.26 files were provided years ago and are thus likely with the C version). Second crossing results were calculated by Scott Callaghan after modifying the python version of UCVM, Summer 2017. Note that the 3rd column is second crossing, not last.

Model	First Upward Crossing	Second Upward Crossing
CVM-S4.26 (cvms5)	None	None
CVM-S4.26-M01 (cvmsi)	First crossing in old tree of opensha-commons/master	Second crossing on opensha-commons/master, now with metadata
CCA06 (cca)	First crossing in old tree of opensha-commons/master	Second crossing on opensha-commons/master, now with metadata

Related Entries

• Basin Query Notes
• UCVM
• UCVM Review
• USR