Difference between revisions of "CTM"
From SCECpedia
Jump to navigationJump to search| (7 intermediate revisions by the same user not shown) | |||
| Line 1: | Line 1: | ||
We are using jupyter notebook developed by Liz Hearn to calculate CTM values at points defined for the GFM. | We are using jupyter notebook developed by Liz Hearn to calculate CTM values at points defined for the GFM. | ||
| + | |||
| + | == GFM/CTM Data Files == | ||
| + | The GFM cells with CTM temperature and heat flow region ID added for each cell. There are to version of this data set. Unsmoothed is CTM temp from the Thatcher CTM model. Smoothed version uses kridging between points to avoid sharp temperature discontinuities. | ||
| + | |||
| + | *Unsmoothed CTM data. | ||
| + | **[http://hypocenter.usc.edu/research/CME/GFM_box_10x10x1km_elev_geo_results_unsmoothed.dat.zip GFM cell data updated with Thatcher CTM temp and heatflow region. Unsmoothed] 16.7Mb Zip files 65Mb unzipped txt file | ||
| + | |||
| + | *Smoothed CTM data | ||
| + | **[http://hypocenter.usc.edu/research/CME/GFM_box_10x10x1km_elev_geo_results_smoothed.dat.zip GFM cell data updated with Thatcher CTM temp and heatflow region. Smoothed] 16.7Mb Zip files 65Mb unzipped txt file | ||
== GFM Points == | == GFM Points == | ||
| Line 7: | Line 16: | ||
The resolution of the GFM 3D voxel data set is 10km x 10km x 1km (depth). | The resolution of the GFM 3D voxel data set is 10km x 10km x 1km (depth). | ||
| − | Each level of the GFM has 8256 points, and the depth goes from -100km to +4km (8256 x 105 = 866880) points. | + | Each level of the GFM has 8256 points, with 960km lon and 860km in lat (96 x 86 = 8256) and the depth goes from -100km to +4km (8256 x 105 = 866880) points. |
Corner Points from GFM file: | Corner Points from GFM file: | ||
| Line 16: | Line 25: | ||
The cell location in this file use the CVM-H convention of elevation relative to sea level in meter, with positive values above sea level, and negative values depth below sea level. The cells extend from -100km to 4km. | The cell location in this file use the CVM-H convention of elevation relative to sea level in meter, with positive values above sea level, and negative values depth below sea level. The cells extend from -100km to 4km. | ||
| + | |||
| + | == KML file showing model coverage regions == | ||
| + | *[http://hypocenter.usc.edu/research/CME/gfm_cvmh_coverage_regions.kml GFM/CVMH Coverage Regions KML File] | ||
== Example GFM File Data Format == | == Example GFM File Data Format == | ||
| Line 40: | Line 52: | ||
== CTM Test Points in GFM format == | == CTM Test Points in GFM format == | ||
| − | The CTM notebook defines a series of test points. The query points are given as depth below surface. | + | The CTM notebook defines a series of test points. The query points are given as depth below surface as a positive number in KM. |
| − | We've converted these test points into the GFM file format, that expresses them as depth below earth surface. We used the elevation information to express the test points in GFM format. We'll use this GFM test file, to confirm our CTM query script is returning the same CTM results, for both smoothed and unsmoothed modes, as the | + | |
| + | We've converted these test points into the GFM file format, that expresses them as depth below earth surface as a negative number. We used the elevation information to express the test points in GFM format. We'll use this GFM test file, to confirm our CTM query script is returning the same CTM results, for both smoothed and unsmoothed modes, as the CTM notebook. | ||
CTM QueryPoint.txt format | CTM QueryPoint.txt format | ||
| Line 70: | Line 83: | ||
</pre> | </pre> | ||
| − | == | + | == Testing GFM cells using CTM QueryPoints == |
| − | + | Unsmoothed results from CTM notebook QueryPoints.txt file | |
| + | <pre> | ||
| + | # Lon Lat Depth Province Temp | ||
| + | -120.000 35.000 5.000 3.000 156.000 | ||
| + | -119.000 35.000 5.000 10.000 131.000 | ||
| + | -118.000 35.000 5.000 9.000 137.000 | ||
| + | -120.000 34.000 6.000 15.000 203.690 | ||
| + | -119.000 34.000 6.000 13.000 184.000 | ||
| + | -118.000 34.000 6.000 7.000 154.000 | ||
| + | -120.000 33.000 8.000 15.000 264.210 | ||
| + | -119.000 33.000 8.000 15.000 264.210 | ||
| + | -118.000 33.000 8.000 15.000 264.210 | ||
| + | </pre> | ||
| + | |||
| + | Unsmoothed results using GFM input files | ||
| + | <pre> | ||
| + | long. lat. Z regionID temp_median heatRegionID elevation | ||
| + | -120.00 35.00 -4532.0 1.000000 156.011 3.0 468.40 | ||
| + | -119.00 35.00 -4556.0 1.000000 131.0 10.0 444.00 | ||
| + | -118.00 35.00 -4290.0 1.000000 137.0 9.0 710.00 | ||
| + | -120.00 34.00 -5820.0 1.000000 203.69 15.0 180.00 | ||
| + | -119.00 34.00 -5820.0 1.000000 184.0 13.0 180.00 | ||
| + | -118.00 34.00 -5532.0 1.000000 154.0 7.0 468.00 | ||
| + | -120.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 | ||
| + | -119.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 | ||
| + | -118.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 | ||
| + | </pre> | ||
| + | |||
| + | Smoothed results using CTM QueryPoints file | ||
| + | <pre> | ||
| + | # Lon Lat Depth Province Temp | ||
| + | -120.000 35.000 5.000 3.000 156.000 | ||
| + | -119.000 35.000 5.000 10.000 131.028 | ||
| + | -118.000 35.000 5.000 9.000 137.000 | ||
| + | -120.000 34.000 6.000 15.000 187.490 | ||
| + | -119.000 34.000 6.000 13.000 176.906 | ||
| + | -118.000 34.000 6.000 7.000 154.000 | ||
| + | -120.000 33.000 8.000 15.000 264.210 | ||
| + | -119.000 33.000 8.000 15.000 264.210 | ||
| + | -118.000 33.000 8.000 15.000 264.210 | ||
| + | </pre> | ||
| + | |||
| + | Smoothed results using CFM format Input file: | ||
| + | <pre> | ||
| + | long. lat. Z regionID temp_median heatRegionID elevation | ||
| + | -120.00 35.00 -4532.0 1.000000 156.011 3.0 468.40 | ||
| + | -119.00 35.00 -4556.0 1.000000 131.028 10.0 444.00 | ||
| + | -118.00 35.00 -4290.0 1.000000 137.0 9.0 710.00 | ||
| + | -120.00 34.00 -5820.0 1.000000 187.49 15.0 180.00 | ||
| + | -119.00 34.00 -5820.0 1.000000 176.906 13.0 180.00 | ||
| + | -118.00 34.00 -5532.0 1.000000 154.0 7.0 468.00 | ||
| + | -120.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 | ||
| + | -119.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 | ||
| + | -118.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 | ||
| + | </pre> | ||
== Heatflow Region IDs and Names == | == Heatflow Region IDs and Names == | ||
| Line 102: | Line 169: | ||
</pre> | </pre> | ||
| + | |||
| + | == Preliminary Research Page == | ||
| + | *[https://www.scec.org/research/ctm CTM Research Page] | ||
== Google Pages == | == Google Pages == | ||
Latest revision as of 17:14, 28 July 2020
We are using jupyter notebook developed by Liz Hearn to calculate CTM values at points defined for the GFM.
Contents
GFM/CTM Data Files
The GFM cells with CTM temperature and heat flow region ID added for each cell. There are to version of this data set. Unsmoothed is CTM temp from the Thatcher CTM model. Smoothed version uses kridging between points to avoid sharp temperature discontinuities.
- Unsmoothed CTM data.
- GFM cell data updated with Thatcher CTM temp and heatflow region. Unsmoothed 16.7Mb Zip files 65Mb unzipped txt file
- Smoothed CTM data
- GFM cell data updated with Thatcher CTM temp and heatflow region. Smoothed 16.7Mb Zip files 65Mb unzipped txt file
GFM Points
AP has created a discretized version of the GFM. The coverage region for the GFM is similar to the CVM-H GTL coverage region. Corner points from the GFM grid are listed below.
The resolution of the GFM 3D voxel data set is 10km x 10km x 1km (depth).
Each level of the GFM has 8256 points, with 960km lon and 860km in lat (96 x 86 = 8256) and the depth goes from -100km to +4km (8256 x 105 = 866880) points.
Corner Points from GFM file:
- -122.681200, 29.709900 -123.205500, 37.335600 -112.481200, 37.412300 -112.864200, 29.767300
A white-space delimited file that defines the GFM region types for each of these points is posted here. This is the updated file that include elevation for each point.
- GFM Cell Points and Data (Warning 13MB zip file download. 65Mb expanded as text file)
The cell location in this file use the CVM-H convention of elevation relative to sea level in meter, with positive values above sea level, and negative values depth below sea level. The cells extend from -100km to 4km.
KML file showing model coverage regions
Example GFM File Data Format
This shows the format of the GFM cell data, before it is updated with CTM heatRegionID and CTM temps. The elevation data is from the GFM topography surface. The last 9 lines are CTM test points, with known results, used to verify that the script used to populate the heatRegionID and Temp are producing the correct results for the GFM points.
long. lat. Z regionID temp_median heatRegionID elevation -122.681222 29.709863 -100000.000000 1.000000 1316.400024 1.000000 -4253.000000 -122.578339 29.714251 -100000.000000 1.000000 1316.400024 1.000000 -4312.000000 -122.475433 29.718560 -100000.000000 1.000000 1316.400024 1.000000 -4310.586426 -122.372505 29.722790 -100000.000000 1.000000 1316.400024 1.000000 -4233.651367 -122.269555 29.726940 -100000.000000 1.000000 1316.400024 1.000000 -4059.687012 -122.166584 29.731010 -100000.000000 1.000000 1316.400024 1.000000 -4157.000000 -122.063592 29.735001 -100000.000000 1.000000 1316.400024 1.000000 -4212.185059 -120.00 35.00 -4532.0 1.000000 1316.400024 1.000000 468.40 -119.00 35.00 -4556.0 1.000000 1316.400024 1.000000 444.00 -118.00 35.00 -4290.0 1.000000 1316.400024 1.000000 710.00 -120.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -119.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -118.00 34.00 -5532.0 1.000000 1316.400024 1.000000 468.00 -120.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -119.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -118.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00
CTM Test Points in GFM format
The CTM notebook defines a series of test points. The query points are given as depth below surface as a positive number in KM.
We've converted these test points into the GFM file format, that expresses them as depth below earth surface as a negative number. We used the elevation information to express the test points in GFM format. We'll use this GFM test file, to confirm our CTM query script is returning the same CTM results, for both smoothed and unsmoothed modes, as the CTM notebook.
CTM QueryPoint.txt format
-120.00 35.00 5.0 -119.00 35.00 5.0 -118.00 35.00 5.0 -120.00 34.00 6.0 -119.00 34.00 6.0 -118.00 34.00 6.0 -120.00 33.00 8.0 -119.00 33.00 8.0 -118.00 33.00 8.0
Same points in GFM cell file format
long. lat. Z regionID temp_median heatRegionID elevation -120.00 35.00 -4532.0 1.000000 1316.400024 1.000000 468.40 -119.00 35.00 -4556.0 1.000000 1316.400024 1.000000 444.00 -118.00 35.00 -4290.0 1.000000 1316.400024 1.000000 710.00 -120.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -119.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -118.00 34.00 -5532.0 1.000000 1316.400024 1.000000 468.00 -120.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -119.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -118.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00
Testing GFM cells using CTM QueryPoints
Unsmoothed results from CTM notebook QueryPoints.txt file
# Lon Lat Depth Province Temp -120.000 35.000 5.000 3.000 156.000 -119.000 35.000 5.000 10.000 131.000 -118.000 35.000 5.000 9.000 137.000 -120.000 34.000 6.000 15.000 203.690 -119.000 34.000 6.000 13.000 184.000 -118.000 34.000 6.000 7.000 154.000 -120.000 33.000 8.000 15.000 264.210 -119.000 33.000 8.000 15.000 264.210 -118.000 33.000 8.000 15.000 264.210
Unsmoothed results using GFM input files
long. lat. Z regionID temp_median heatRegionID elevation -120.00 35.00 -4532.0 1.000000 156.011 3.0 468.40 -119.00 35.00 -4556.0 1.000000 131.0 10.0 444.00 -118.00 35.00 -4290.0 1.000000 137.0 9.0 710.00 -120.00 34.00 -5820.0 1.000000 203.69 15.0 180.00 -119.00 34.00 -5820.0 1.000000 184.0 13.0 180.00 -118.00 34.00 -5532.0 1.000000 154.0 7.0 468.00 -120.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 -119.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 -118.00 33.00 -7820.0 1.000000 264.21 15.0 180.00
Smoothed results using CTM QueryPoints file
# Lon Lat Depth Province Temp -120.000 35.000 5.000 3.000 156.000 -119.000 35.000 5.000 10.000 131.028 -118.000 35.000 5.000 9.000 137.000 -120.000 34.000 6.000 15.000 187.490 -119.000 34.000 6.000 13.000 176.906 -118.000 34.000 6.000 7.000 154.000 -120.000 33.000 8.000 15.000 264.210 -119.000 33.000 8.000 15.000 264.210 -118.000 33.000 8.000 15.000 264.210
Smoothed results using CFM format Input file:
long. lat. Z regionID temp_median heatRegionID elevation -120.00 35.00 -4532.0 1.000000 156.011 3.0 468.40 -119.00 35.00 -4556.0 1.000000 131.028 10.0 444.00 -118.00 35.00 -4290.0 1.000000 137.0 9.0 710.00 -120.00 34.00 -5820.0 1.000000 187.49 15.0 180.00 -119.00 34.00 -5820.0 1.000000 176.906 13.0 180.00 -118.00 34.00 -5532.0 1.000000 154.0 7.0 468.00 -120.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 -119.00 33.00 -7820.0 1.000000 264.21 15.0 180.00 -118.00 33.00 -7820.0 1.000000 264.21 15.0 180.00
Heatflow Region IDs and Names
ID's for the 14 original CTM heat flow regions, from Thatcher and Chapman (2020): 1 WBR (Western Basin and Range) 2 SN (Sierra Nevada) 3 CCR (Central Coast Ranges) 4 PR (Peninsular Ranges) 5 ST (Outer Salton Trough, aka oST) 6 WTR (Western Transverse Ranges) 7 LA (Los Angeles) 8 SG (San Gabriels) 9 MD (Mojave Desert) 10 GV (Great Valley) 11 ETR (Eastern Transverse Ranges) 12 iST (Inner Salton Trough, also called STi5 in some files) 13 iCB (Inner Continental Borderland) 14 VB (Ventura Basin) (the following were added to expand CTM region to cover the Geologic Framework region) 15 Ocean 16 CCR extension 17 GV extension 18 SN extension 19 WBR extension 20 ST extension 21 PR extension
Preliminary Research Page
Google Pages
Google Colab notebook: CTM_tools.ipynb
files used by the query code:
