CTM

From SCECpedia
Revision as of 23:26, 27 July 2020 by Maechlin (talk | contribs)
Jump to navigationJump to search

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.

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.

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

Google Pages

Google Colab notebook: CTM_tools.ipynb

files used by the query code:

Related Entries

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