Difference between revisions of "CyberShake interpolation"

From SCECpedia
Jump to navigationJump to search
Line 39: Line 39:
  
 
I chose two metrics by which to compare the accuracy of my interpolations - max percent difference and average percent difference.
 
I chose two metrics by which to compare the accuracy of my interpolations - max percent difference and average percent difference.
 +
 +
=== Results ===
  
 
== Intensity Measures ==
 
== Intensity Measures ==

Revision as of 18:35, 21 June 2024

This page documents research into techniques for interpolating various CyberShake data products, performed by 2024 SOURCES intern Amelia Kratzer.

Source Code

Hazard Curves

Procedure

Selection of Sites

  • To start, I selected 4 sites (input sites), all in the basin, that formed a 10km by 10km square inside which there was at least 1 site to interpolate (interp site). This way I could compare the simulated values for that site to the interpolated to check the accuracy of my interpolation.
    • 1st chosen area: USC as there are three sites on the exterior for 1D linear interpolation and two sites on the interior for 2D linear interpolation
    • Square formed by sites S345, S387, S389, S347
    • 2nd chosen area: CSDH
    • Square formed by sites S385, S429, S431, S387
  • Next, I selected two locations where some of the sites were not in the basin, to test how different velocity structures may affect the accuracy of the interpolation.
    • 1st location: STG
    • Square formed by sites S491, S632, S634, S593
    • 2nd location: SGRTT
    • Square formed by sites S437, S478, S480, S439

Linear Interpolation of Hazard Curves

  • Download and plot the input sites' hazard curves.
  • Query CyberShake database to get the X shaking values and Y probability values for input and interp sites.
  • Use UTM to convert locations of sites to meters from lat/lon.
  • Calculate the X’ and Y’ distances to my interpolated site from the left and top edge of my square since the 10km by 10km square is on a rotated axis.
  • For each shaking value, use the probabilities of shaking at the input sites and these X’ and Y’ distances to perform bilinear interpolation.
  • Plot the (shaking value, interpolated probability) points to form the interpolated hazard curve.
  • Overlay the interpolated hazard curve and the simulated hazard curve on single plot for comparison.

Error Calculations
I compared the simulated probabilities to the interpolated probabilities for the given site, and calculated the absolute percent difference at each acceleration value. Percent difference accurately models the accuracy of the hazard curve since it reflects the larger differences in the simulated and interpolated results at the higher acceleration values.

I chose two metrics by which to compare the accuracy of my interpolations - max percent difference and average percent difference.

Results

Intensity Measures

Procedure

Results