SCECpedia - User contributions [en]

https://strike.scec.org/scecwiki/api.php?action=feedcontributions&feedformat=atom&user=Mei SCECpedia - User contributions [en] 2026-10-03T22:02:22Z User contributions MediaWiki 1.34.2 https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=31004 MUSCAL 2026-09-28T21:25:42Z

<p>Mei: /* TileDB version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> ncdump -h model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> <br /> <pre><br /> netcdf model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd {<br /> dimensions:<br /> depth = 210 ;<br /> latitude = 1251 ;<br /> longitude = 1301 ;<br /> variables:<br /> float vp(depth, latitude, longitude) ;<br /> vp:_FillValue = NaNf ;<br /> vp:units = "m/s" ;<br /> float vs(depth, latitude, longitude) ;<br /> vs:_FillValue = NaNf ;<br /> vs:units = "m/s" ;<br /> float rho(depth, latitude, longitude) ;<br /> rho:_FillValue = NaNf ;<br /> rho:units = "kg/m^3" ;<br /> float longitude(longitude) ;<br /> longitude:_FillValue = NaNf ;<br /> longitude:units = "degrees_east" ;<br /> float latitude(latitude) ;<br /> latitude:_FillValue = NaNf ;<br /> latitude:units = "degrees_north" ;<br /> float depth(depth) ;<br /> depth:_FillValue = NaNf ;<br /> depth:units = "meters" ;<br /> <br /> // global attributes:<br /> :_NCProperties = "version=2,netcdf=4.9.2,hdf5=1.14.3" ;<br /> }<br /> </pre><br /> <br /> longitudes, latitudes and depth are array of floats while vp, vp, and rho are 3 dimensional data array<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> <pre><br /> tiledb_dump -a model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.tiledb <br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"<br /> </pre><br /> <br /> code is in SCECcode/muscaltdb.git.<br /> <br /> tools/convert_nc_to_tiledb.c : converting NetCDF data to TileDB dataDimensions: (depth: 210, latitude: 1251, longitude: 1301)<br /> <br /> <pre><br /> Coordinates:<br /> * longitude (longitude) float32 -126.0 -126.0 -126.0 ... -113.0 -113.0 -113.0<br /> * latitude (latitude) float32 31.0 31.01 31.02 31.03 ... 43.48 43.49 43.5<br /> * depth (depth) float32 0.0 50.0 100.0 150.0 ... 9.7e+04 9.8e+04 9.9e+04<br /> Data variables:<br /> vp (depth, latitude, longitude) float32 dask.array<chunksize=(20, 200, 200), meta=np.ndarray><br /> vs (depth, latitude, longitude) float32 dask.array<chunksize=(20, 200, 200), meta=np.ndarray><br /> rho (depth, latitude, longitude) float32 dask.array<chunksize=(20, 200, 200), meta=np.ndarray><br /> <br /> </pre><br /> <br /> <br /> Note:<br /> <br /> <pre><br /> <br /> import tiledb<br /> <br /> with tiledb.open("model_array", "r") as A:<br /> res = A[100, 200, 50] # (i, j, k)<br /> print(res)<br /> <br /> and vp=res["vp"][l0]<br /> grab the data point using a offset<br /> <br /> or,<br /> <br /> with tiledb.open("model_array", "r", attr="vp") as A:<br /> vp = A[100, 200, 50][0]<br /> <br /> </pre></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=31003 MUSCAL 2026-09-28T21:23:43Z

<p>Mei: /* TileDB version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> ncdump -h model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> <br /> <pre><br /> netcdf model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd {<br /> dimensions:<br /> depth = 210 ;<br /> latitude = 1251 ;<br /> longitude = 1301 ;<br /> variables:<br /> float vp(depth, latitude, longitude) ;<br /> vp:_FillValue = NaNf ;<br /> vp:units = "m/s" ;<br /> float vs(depth, latitude, longitude) ;<br /> vs:_FillValue = NaNf ;<br /> vs:units = "m/s" ;<br /> float rho(depth, latitude, longitude) ;<br /> rho:_FillValue = NaNf ;<br /> rho:units = "kg/m^3" ;<br /> float longitude(longitude) ;<br /> longitude:_FillValue = NaNf ;<br /> longitude:units = "degrees_east" ;<br /> float latitude(latitude) ;<br /> latitude:_FillValue = NaNf ;<br /> latitude:units = "degrees_north" ;<br /> float depth(depth) ;<br /> depth:_FillValue = NaNf ;<br /> depth:units = "meters" ;<br /> <br /> // global attributes:<br /> :_NCProperties = "version=2,netcdf=4.9.2,hdf5=1.14.3" ;<br /> }<br /> </pre><br /> <br /> longitudes, latitudes and depth are array of floats while vp, vp, and rho are 3 dimensional data array<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> <pre><br /> tiledb_dump -a model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.tiledb <br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"<br /> </pre><br /> <br /> code is in SCECcode/muscaltdb.git.<br /> <br /> tools/convert_nc_to_tiledb.c : converting NetCDF data to TileDB dataDimensions: (depth: 210, latitude: 1251, longitude: 1301)<br /> <br /> <pre><br /> Coordinates:<br /> * longitude (longitude) float32 -126.0 -126.0 -126.0 ... -113.0 -113.0 -113.0<br /> * latitude (latitude) float32 31.0 31.01 31.02 31.03 ... 43.48 43.49 43.5<br /> * depth (depth) float32 0.0 50.0 100.0 150.0 ... 9.7e+04 9.8e+04 9.9e+04<br /> Data variables:<br /> vp (depth, latitude, longitude) float32 dask.array<chunksize=(20, 200, 200), meta=np.ndarray><br /> vs (depth, latitude, longitude) float32 dask.array<chunksize=(20, 200, 200), meta=np.ndarray><br /> rho (depth, latitude, longitude) float32 dask.array<chunksize=(20, 200, 200), meta=np.ndarray><br /> <br /> </pre></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=31002 MUSCAL 2026-09-28T20:50:08Z

<p>Mei: /* TileDB version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> ncdump -h model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> <br /> <pre><br /> netcdf model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd {<br /> dimensions:<br /> depth = 210 ;<br /> latitude = 1251 ;<br /> longitude = 1301 ;<br /> variables:<br /> float vp(depth, latitude, longitude) ;<br /> vp:_FillValue = NaNf ;<br /> vp:units = "m/s" ;<br /> float vs(depth, latitude, longitude) ;<br /> vs:_FillValue = NaNf ;<br /> vs:units = "m/s" ;<br /> float rho(depth, latitude, longitude) ;<br /> rho:_FillValue = NaNf ;<br /> rho:units = "kg/m^3" ;<br /> float longitude(longitude) ;<br /> longitude:_FillValue = NaNf ;<br /> longitude:units = "degrees_east" ;<br /> float latitude(latitude) ;<br /> latitude:_FillValue = NaNf ;<br /> latitude:units = "degrees_north" ;<br /> float depth(depth) ;<br /> depth:_FillValue = NaNf ;<br /> depth:units = "meters" ;<br /> <br /> // global attributes:<br /> :_NCProperties = "version=2,netcdf=4.9.2,hdf5=1.14.3" ;<br /> }<br /> </pre><br /> <br /> longitudes, latitudes and depth are array of floats while vp, vp, and rho are 3 dimensional data array<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB<br /> <br /> <pre><br /> tiledb_dump -a xx_array<br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"<br /> </pre></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=31001 MUSCAL 2026-09-28T20:49:30Z

<p>Mei: /* NetCDF version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> ncdump -h model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> <br /> <pre><br /> netcdf model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd {<br /> dimensions:<br /> depth = 210 ;<br /> latitude = 1251 ;<br /> longitude = 1301 ;<br /> variables:<br /> float vp(depth, latitude, longitude) ;<br /> vp:_FillValue = NaNf ;<br /> vp:units = "m/s" ;<br /> float vs(depth, latitude, longitude) ;<br /> vs:_FillValue = NaNf ;<br /> vs:units = "m/s" ;<br /> float rho(depth, latitude, longitude) ;<br /> rho:_FillValue = NaNf ;<br /> rho:units = "kg/m^3" ;<br /> float longitude(longitude) ;<br /> longitude:_FillValue = NaNf ;<br /> longitude:units = "degrees_east" ;<br /> float latitude(latitude) ;<br /> latitude:_FillValue = NaNf ;<br /> latitude:units = "degrees_north" ;<br /> float depth(depth) ;<br /> depth:_FillValue = NaNf ;<br /> depth:units = "meters" ;<br /> <br /> // global attributes:<br /> :_NCProperties = "version=2,netcdf=4.9.2,hdf5=1.14.3" ;<br /> }<br /> </pre><br /> <br /> longitudes, latitudes and depth are array of floats while vp, vp, and rho are 3 dimensional data array<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB<br /> <br /> <br /> tiledb_dump -a xx_array<br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=31000 MUSCAL 2026-09-28T20:45:34Z

<p>Mei: /* NetCDF version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> ncdump -h model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> <br /> netcdf model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd {<br /> dimensions:<br /> depth = 210 ;<br /> latitude = 1251 ;<br /> longitude = 1301 ;<br /> variables:<br /> float vp(depth, latitude, longitude) ;<br /> vp:_FillValue = NaNf ;<br /> vp:units = "m/s" ;<br /> float vs(depth, latitude, longitude) ;<br /> vs:_FillValue = NaNf ;<br /> vs:units = "m/s" ;<br /> float rho(depth, latitude, longitude) ;<br /> rho:_FillValue = NaNf ;<br /> rho:units = "kg/m^3" ;<br /> float longitude(longitude) ;<br /> longitude:_FillValue = NaNf ;<br /> longitude:units = "degrees_east" ;<br /> float latitude(latitude) ;<br /> latitude:_FillValue = NaNf ;<br /> latitude:units = "degrees_north" ;<br /> float depth(depth) ;<br /> depth:_FillValue = NaNf ;<br /> depth:units = "meters" ;<br /> <br /> // global attributes:<br /> :_NCProperties = "version=2,netcdf=4.9.2,hdf5=1.14.3" ;<br /> }<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB<br /> <br /> <br /> tiledb_dump -a xx_array<br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30999 MUSCAL 2026-09-28T20:34:35Z

<p>Mei: /* NetCDF version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> ncdump -h /var/www/html/CVM_DATASET_DIRECTORY/model/muscal/model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.<br /> nc<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB<br /> <br /> <br /> tiledb_dump -a xx_array<br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30998 MUSCAL 2026-09-28T20:33:13Z

<p>Mei: /* TileDB version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB<br /> <br /> <br /> tiledb_dump -a xx_array<br /> python -c "import tiledb; print(tiledb.open('muscal_canvas_tiledb_array').schema)"</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30997 MUSCAL 2026-09-28T20:29:00Z

<p>Mei: /* TileDB version of MUSCAL */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == NetCDF version of MUSCAL ==<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30996 MUSCAL 2026-09-28T20:27:46Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> The surface and kdtree code are in SCECcode/muscalnc.git, src directory: muscalnc_surface.c, kdtree_util.c<br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30995 MUSCAL 2026-09-28T19:54:11Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> <br /> == TileDB version of MUSCAL ==<br /> <br /> Prototype TileDB as the base data storage solution, translation and conversion code were implemented and verified with code generate using various AI tools, ie. Gemini Pro, copilot and chapGPT.<br /> <br /> MUSCALTDB</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30994 MUSCAL 2026-09-28T19:51:07Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree), which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model); Query the Kdtree for the nearest neighbors.</li><br /> </ol></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30993 MUSCAL 2026-09-28T19:50:18Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>Compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree) which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model): Query the Kdtree for the nearest neighbors.</li><br /> </ol></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30992 MUSCAL 2026-09-28T19:49:57Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree) which is good for binary tree search</li><br /> <li>If the model query returns NaNs (indicating the point is beyond the bounds of the 3D model): Query the Kdtree for the nearest neighbors.</li><br /> </ol></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30991 MUSCAL 2026-09-28T19:49:21Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree) which is good for binary tree search</li><br /> <ol?<br /> <li>If the model returns NaNs (indicating the point is beyond the bounds of the 3D model): Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> </ol></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30990 MUSCAL 2026-09-28T19:47:19Z

<p>Mei: /* Extended Model */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM with nearest neighbor extrapolation to generate physically consistent properties beyond defined model boundary . <br /> <br /> <ol><br /> <li>Iterate through each layer of MUSCAL dataset and collect the material properties of all edge points</li><br /> <li>compose a surface point dataset from the collected points</li><br /> <li>During UCVM query process, load the surface dataset and create a kdtree structure (k-dimensional tree) which is good for binary tree search</li><br /> <li>If the model returns NaNs (indicating the point is beyond the bounds of the 3D model):<br /> <ol><br /> <li>Query the Kdtree for the nearest neighbors.</li><br /> </ol><br /> </li><br /> </ol></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30989 MUSCAL 2026-09-28T19:36:28Z

<p>Mei: /* Model overview */</p> <hr /> <div>This page describes the Multi-Scale CALifornia (MUSCAL) velocity model and its integration into UCVM.<br /> <br /> == Model overview ==<br /> <br /> The MUSCAL Vp and Vs model was created by Te-Yang Yeh, Kim B. Olsen and Yehuda Ben-Zion in 2025-6 by starting with the CANVAS tomography model as a base and then integrating multiple regional and local high-resolution models. <br /> <br /> === Model Components ===<br /> <br /> The model was built by starting with a base model (the CANVAS (Doody et al., 2023) tomography model). Next, regional and local models are added where they improve the fit, based on M4 simulations. These simulations were run to 1 Hz with topography using a minimum Vs of 450 m/s. Finally, a spatially varying taper following the Ely approach is added.<br /> <br /> ==== Regional Models ====<br /> <br /> Regional models evaluated for inclusion are:<br /> <ul><br /> <li>Northern California<br /> <ul><br /> <li>Lin et al. (2010)</li><br /> <li>Furlong et al. (2024)</li><br /> <li>Guo et al. (2024)</li><br /> </ul></li><br /> <li>Central California<br /> <ul><br /> <li>CCA-06</li><br /> <li>CVM-H v15.1.1</li><br /> <li>CVM-S4.26</li><br /> <li>Guo et al. (2024)</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> <li>Southern California<br /> <ul><br /> <li>Fang et al. (2022)</li><br /> <li>CVM-S4.26</li><br /> <li>CVM-H v15.1.1</li><br /> <li>Berg et al. (2021)</li><br /> </ul></li><br /> </ul><br /> <br /> ==== Local Models ====<br /> <br /> Below is a list of local high-resolution models that were evaluated for inclusion.<br /> <br /> Northern and Central California:<br /> <ul><br /> <li>Ridgecrest regional (Li and Ben-Zion (2024))</li><br /> <li>Ridgecrest fault zone (Zhou et al. (2022))</li><br /> <li>San Joaquin basin (Shaw & Plesch (2016))</li><br /> <li>SFCVM v21.1 (Hirakawa & Aagaard (2022))</li><br /> </ul><br /> <br /> Southern California:<br /> <ul><br /> <li>San Jacinto fault zone (Fang et al. (2019))</li><br /> <li>San Gabriel-San Bernardino Basin (Li et al. (2023))</li><br /> <li>Imperial Valley (Persaud et al. (2016))</li><br /> <li>Coachella Valley (Ajala et al. (2019))</li><br /> <li>Santa Maria Basin (Plesch et al. (2020))</li><br /> </ul><br /> <br /> ==== Taper ====<br /> <br /> MUSCAL includes a low-velocity taper using the Ely-Jordan approach. Vs30 for this method is produced by merging the Thompson (2018) Vs30 map with the [https://earthquake.usgs.gov/data/vs30/ USGS global Vs30 dataset].<br /> <br /> The depth of the taper varies depending on the site and ranges from 300m to 1500m.<br /> <br /> == Delivered Model ==<br /> <br /> The final version of MUSCAL is delivered in netCDF format. The dimensions of the model in grid points are 1251 (lat) x 1301 (lon) x 210 (depth). Grid points have a spacing of 0.01 degrees in the horizontal; the depth is variable with denser sampling near the surface. The model is interpolated using trilinear interpolation.<br /> <br /> Details about the model creation are provided in the paper ([https://essopenarchive.org/doi/full/10.22541/essoar.15001929/v1 preprint]) and the [https://zenodo.org/records/19243477 Zenodo page].<br /> <br /> Below is a plot of the 3D model extent.<br /> <br /> {|<br /> | [[File:MUSCAL_region_extent.png|thumb|400px]]<br /> |}<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitude:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]<br /> <br /> == Extended Model ==<br /> <br /> The 3D MUSCAL model is defined over the same volume as the CANVAS base model. However, some simulation regions may extend beyond the boundaries, like the blue corners in the example CyberShake volume for USC below: <br /> <br /> {|<br /> | [[File:USC_muscal_horiz_slice_z0.png|thumb|400px]]<br /> |}<br /> <br /> To solve this in a consistent and reproducible way, we will extend the MUSCAL model through UCVM. The initial plan is to use the following algorithm:<br /> <ol><br /> <li>Query 3D MUSCAL model at a given point.</li><br /> <li>If the model returns NaNs (indicating the point is beyond the bounds of the 3D model):<br /> <ol><br /> <li>Query the 1D MUSCAL model (Te-Yang to provide this).</li><br /> <li>Determine Vs30 from the USGS Global Vs30 model.</li><br /> <li>Apply the Ely taper (LVT) to the top 300m.</li><br /> </ol><br /> </li><br /> </ol></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:MUSCAL_test_points_shallow.final.txt&diff=30735 File:MUSCAL test points shallow.final.txt 2026-04-01T23:08:50Z

<p>Mei: </p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:MUSCAL_test_points_deep.final.txt&diff=30734 File:MUSCAL test points deep.final.txt 2026-04-01T23:08:28Z

<p>Mei: </p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30733 MUSCAL 2026-04-01T23:08:08Z

<p>Mei: /* Validation */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.final.txt]]</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30732 MUSCAL 2026-04-01T23:05:32Z

<p>Mei: /* Validation */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]<br /> <br /> <br /> Final result from ucvm_query, (with interpolation)<br /> <br /> [[FILE:MUSCAL_test_points_deep.final]]<br /> [[FILE:MUSCAL_test_points_shallow.final]]</div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:MUSCAL_test_points_shallow.txt&diff=30731 File:MUSCAL test points shallow.txt 2026-04-01T23:03:32Z

<p>Mei: </p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:MUSCAL_test_points_deep.txt&diff=30730 File:MUSCAL test points deep.txt 2026-04-01T23:02:43Z

<p>Mei: </p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30729 MUSCAL 2026-04-01T23:01:15Z

<p>Mei: /* Validation */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30728 MUSCAL 2026-04-01T23:00:29Z

<p>Mei: /* Validation */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32<br /> <br /> [[FILE:MUSCAL_test_points_deep.txt]]<br /> [[FILE:MUSCAL_test_points_shallow.txt]]</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30727 MUSCAL 2026-04-01T22:58:53Z

<p>Mei: /* ... */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> === Validation ===<br /> <br /> <br /> Several steps of pre/post processing are done to original MUSCAL model(netcdf format) in order to <br /> incorporate it into UCVM and CVM explorer.<br /> <br /> Because of the speed of netcdf-C code is too slow to support in-time nature of the explorer and would <br /> like to have interpolation on query,<br /> <br /> * Number of depth layer is reduced with deeper layers merged into fewer layers<br /> * preprocessing the netcdf data into binary data files<br /> * data are all in float32</div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp-1000_matprops.json.txt&diff=30726 File:Muscal small interp-1000 matprops.json.txt 2026-04-01T22:41:55Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp-1000 matprops.json.txt</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp-500_matprops.json.txt&diff=30725 File:Muscal small interp-500 matprops.json.txt 2026-04-01T22:41:36Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp-500 matprops.json.txt</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp-100_matprops.json.txt&diff=30724 File:Muscal small interp-100 matprops.json.txt 2026-04-01T22:40:13Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp-100 matprops.json.txt</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp-50_matprops.json.txt&diff=30723 File:Muscal small interp-50 matprops.json.txt 2026-04-01T22:39:50Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp-50 matprops.json.txt</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp_1000.png&diff=30722 File:Muscal small interp 1000.png 2026-04-01T22:39:29Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp 1000.png</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp_500.png&diff=30721 File:Muscal small interp 500.png 2026-04-01T22:39:04Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp 500.png</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp_100.png&diff=30720 File:Muscal small interp 100.png 2026-04-01T22:38:35Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp 100.png</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp_50.png&diff=30719 File:Muscal small interp 50.png 2026-04-01T22:37:50Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp 50.png</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30718 MUSCAL 2026-04-01T22:37:17Z

<p>Mei: /* another small dataset in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> ==== ... ====</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30717 MUSCAL 2026-04-01T22:36:48Z

<p>Mei: /* Another small dataset in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> ==== ... ====</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30716 MUSCAL 2026-04-01T22:36:29Z

<p>Mei: /* another small dataset in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== another small dataset in netcdf ====<br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> m</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30715 MUSCAL 2026-04-01T22:35:11Z

<p>Mei: /* Another small dataset in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== another small dataset in netcdf ====<br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> m</div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30714 MUSCAL 2026-04-01T22:34:48Z

<p>Mei: /* small data in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> ==== another small dataset in netcdf ====<br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30713 MUSCAL 2026-04-01T19:23:05Z

<p>Mei: /* Datasets */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command use:<br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth_1000.png -C 'Multi-Scale Statewide California Velocity Model'<br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30712 MUSCAL 2026-04-01T19:20:56Z

<p>Mei: Undo revision 30707 by Mei (talk)</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30711 MUSCAL 2026-04-01T19:18:33Z

<p>Mei: Undo revision 30710 by Mei (talk)</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command used: <br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth1000.png -C 'Multi-Scale Statewide California Velocity Model' <br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_depth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_depth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_depth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> * the interp used has incorrect axis -- ignore<br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30710 MUSCAL 2026-04-01T19:18:07Z

<p>Mei: Undo revision 30709 by Mei (talk)</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command used: <br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth1000.png -C 'Multi-Scale Statewide California Velocity Model' <br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_depth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_deth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> * the interp used has incorrect axis -- ignore<br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30709 MUSCAL 2026-04-01T19:17:50Z

<p>Mei: Undo revision 30708 by Mei (talk)</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command used: <br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth1000.png -C 'Multi-Scale Statewide California Velocity Model' <br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_depth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_depth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_depth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> * the interp used has incorrect axis -- ignore<br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30708 MUSCAL 2026-04-01T19:14:12Z

<p>Mei: /* big data in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command used: <br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth1000.png -C 'Multi-Scale Statewide California Velocity Model' <br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_depth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_deth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> * the interp used has incorrect axis -- ignore<br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30707 MUSCAL 2026-04-01T19:13:08Z

<p>Mei: /* Datasets */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> <br /> command used: <br /> plot_depth_profile.py -n $UCVM_INSTALL_PATH/conf/ucvm.conf -i $UCVM_INSTALL_PATH -d vs -c muscal <br /> -o muscal_small_depth1000.png -C 'Multi-Scale Statewide California Velocity Model' <br /> -v 1000 -b 0 -s 36.5054,-119.0587 -e 30000<br /> <br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_depth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_depth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_depth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> * the interp used has incorrect axis -- ignore<br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=MUSCAL&diff=30706 MUSCAL 2026-03-30T20:36:05Z

<p>Mei: /* small data in netcdf */</p> <hr /> <div>== about MUSCAL ==<br /> <br /> From Te-Yang Yeh and Yehuda Ben-Zion<br /> <br /> The Multi-Scale CALifornia (MUSCAL) statewide Vp and Vs velocity models provide high-quality integrated description of seismic structures across the state. Starting with the CANVAS base model (Doody et al., 2023), MUSCAL incorporates multiple regional and local velocity datasets into a unified structure, capturing features ranging from broad crustal-mantle structures to fine-scale local anomalies such as sedimentary basins.<br /> <br /> To ensure quality, the merged multi-scale models underwent a data-informed refinement process guided by simulations of small validation events. A key feature of MUSCAL is the inclusion of a locally optimized near-surface low-velocity taper (LVT), specifically designed to better represent under-resolved shallow structures and improve the accuracy of ground-motion predictions.<br /> <br /> <br /> === Datasets ===<br /> <br /> zone: 11<br /> <br /> ==== big data in netcdf ====<br /> <br /> model_MSCAL_CANVAS_dll0.01_dz50_cmpd.nc<br /> 4.5G<br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:671 from 0 to 100,000 <br /> 50 increments until 30,000<br /> 1000 increments til 100000<br /> <br /> Plot Depth profile at 36.5054,-119.0587 in different step increments. No interpretation and access data using nc api<br /> <br /> {|<br /> | [[FILE:muscal_big_deth_50.png|thumb|300px|muscal big 50 vs]]<br /> | [[FILE:muscal_big_deth_100.png|thumb|300px|muscal big 100 vs]]<br /> | [[FILE:muscal_big_depth_500.png|thumb|300px|muscal big 500 vs]] <br /> | [[FILE:muscal_big_deth_1000.png|thumb|300px|muscal big 1000 vs]]<br /> |}<br /> <br /> <br /> [[FILE:muscal_big_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_big_depth-1000_matprops.json.txt]]<br /> <br /> ==== small data in netcdf ====<br /> model_MUSCAL_CANVAS_dll0.01_vardz_cmpd.nc <br /> 2.1G <br /> <br /> longitdue:1301 from -126 to -113 <br /> latitude:1251 from 31 to 43.5<br /> depth:210 from 0 to 99,000 <br /> 50 increments upto 3000 <br /> 100 increments upto 5000<br /> 250 increments upto 10000<br /> 500 increments upto 30000<br /> 1000 increments upto 99000<br /> <br /> <br /> Depth profiles in different steps : 50m,100m,500m,1000m<br /> <br /> <br /> Direct from netcdf as external file<br /> <br /> {|<br /> |[[FILE:muscal_small_depth_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_deth_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_depth_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_deth_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_depth-50_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-100_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-500_matprops.json.txt]]<br /> [[FILE:muscal_small_depth-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data no interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_no_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_no_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_no_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_no_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_no_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_no_interp-1000_matprops.json.txt]]<br /> <br /> <br /> Loaded in-memory as binary data with interpolation<br /> <br /> {|<br /> |[[FILE:muscal_small_interp_50.png|thumb|300px|muscal small 50 vs]] <br /> |[[FILE:muscal_small_interp_100.png|thumb|300px|muscal small 100 vs]]<br /> |[[FILE:muscal_small_interp_500.png|thumb|300px|muscal small 500 vs]] <br /> |[[FILE:muscal_small_interp_1000.png|thumb|300px|muscal small 1000 vs]]<br /> |}<br /> <br /> [[FILE:muscal_small_interp-50_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-100_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-500_matprops.json.txt]]<br /> [[FILE:muscal_small_interp-1000_matprops.json.txt]]<br /> <br /> <br /> ==== Another small dataset in netcdf ====<br /> <br /> <br /> model_MUSCAL_CANVAS_dll0.01_vardz_float32_cmpd.nc<br /> 1.4G<br /> <br /> All longitude, latitude, and depth points are now saved as float32, as well as the vp, vs, and rho. </div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp-100_matprops.json.txt&diff=30705 File:Muscal small interp-100 matprops.json.txt 2026-03-30T19:58:20Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp-100 matprops.json.txt</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp-50_matprops.json.txt&diff=30704 File:Muscal small interp-50 matprops.json.txt 2026-03-30T19:57:55Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp-50 matprops.json.txt</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp_100.png&diff=30703 File:Muscal small interp 100.png 2026-03-30T19:55:35Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp 100.png</p> <hr /> <div></div>

Mei https://strike.scec.org/scecwiki/index.php?title=File:Muscal_small_interp_50.png&diff=30702 File:Muscal small interp 50.png 2026-03-30T19:54:16Z

<p>Mei: Mei uploaded a new version of File:Muscal small interp 50.png</p> <hr /> <div></div>

Mei