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 CTm notebook.
+
 
 +
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>
  
== KML file showing model coverage regions ==
+
== Testing GFM cells using CTM QueryPoints ==
*[http://hypocenter.usc.edu/research/CME/gfm_cvmh_coverage_regions.kml GFM/CVMH Coverage Regions KML File]
+
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.

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

Preliminary Research Page

Google Pages

Google Colab notebook: CTM_tools.ipynb

files used by the query code:

Related Entries