WG02 Hazard Calcs - Instructions for figures


This page explains how to reproduce the figures in the paper:

Hazard Calculations for the WGCEP-2002 Forecast Using OpenSHA and Distributed Object Technologies

By Edward H. Field, Nitin Gupta, Vipin Gupta, Michael Blanpied, Philip Maechling, and Thomas H Jordan

Seismological Research Letters 76, 161-167.

Download PDF (WG02_HazCalcsPaper_V15.pdf)


You will need Internet access to reproduce these figures.

This does not describe how to reproduce Figure 1 (which shows the conceptual framework of OpenSHA) because that figure was made using PowerPoint.

First, download the application (Java jar file) at:

Note to Mac/Safari users: you must right-click the Jar-file link below (or hold down the option key if your mouse only has one button), otherwise the Jar-file contents will be shown in a browser window (not a pretty site)

Hazard Curve Calculators:

Now launch the application by double clicking (on Mac or Windows) or by launching from the command line as follows:

java -jar HazardCurveServerModeApp.jar

Read the section “How to Run Java Applications” if you have problems.

Note: you can change the size of sub-windows in the application by dragging the boundaries.


FIGURE 2:

1) Choose the IMR (or attenuation relationship) by choosing Campbell and Bozorgnia (2003) from the first list in the “Set IMR” window.

2) Set the IMT (intensity-measure type) by choosing “PGA” from the first list in the “Set IMT” window,.

3) Set the latitude and longitude in the “Set Site Params” window to be 37.8 and -122.417, respectively. This is a site near downtown San Francisco that happens to be a class “B” site according to the site-conditions map of Wills et al. (2000). An alternative way to set this latitude and longitude is to select “Sites of Interest” from the “Control Panels” list at the bottom of the application, and then choose “San Francisco Class B” from the window that pops up.

4) We now want to set the “Campbell-2002 Site Type” parameter below the latitude and longitude. We can set this to any of the options listed. However, we’ll set this to a value consistent with the Wills et al. (2000) classification map by doing the following: from the “Control Panels” list at the bottom of the application, choose “Set Site Params from Web Services”. In the window that pops up, click “Set Params For All IMRs”. Doing this requests the Wills et al. (2000) classification for our location from a Web Service, and the value that is returned (“B”) is translated into the appropriate value for the “Campbell-2002 Site Type” (the parameter should now have been automatically set to “Firm Rock”). In fact, because we clicked the “Set Params For All IMRs” button in the control panel, the site parameters for all IMRs (attenuation relationships) were automatically set; these settings are needed for creating Figure 3 below.

5) Choose an ERF (Earthquake Rupture Forecast) by doing the following: click the “ERF & Time Span” tab at the top of the application, and under the “Eqk Rup Forecast” list at the top of the sub-window, choose “WG02 Fortran Wrapped ERF List” (you might need to scroll down to see this option).

6) In the “Epistemic List Control” that automatically pops up, click the box next to “Plot Average” (the box next to “Plot all curves” should already be selected). Also choose “Plot Fractiles” from the list below, and click “Update Fractile List” at the bottom of the popup window. You can close this popup window when done.

7) In the main application, set the various ERF’s adjustable parameters (listed below the chosen ERF) as follows:

Rupture Offset = 10.0
Fault Discretization = 2.0
Delta Mag = 0.2
Backgound Seismicity = “Exclude”
GR Tail Seismicity = “Include”
Num Realizations = 100

8) Initiate the computation by clicking “Compute” at the bottom of the application. It will take several minutes to generate all the curves for the 100 realizations. The progress bar should keep you posted.

9) Convert to a log-log plot by clicking both the "YLog" and "XLog" check-boxes at the bottom of the application.

10) change the color scheme by clicking “Plot Prefs” at the bottom of the application and in the popup window set the following for each Dataset:

DATASET #1 Curves: Color as Gray and Size as 2
DATASET #1 Fractiles: Color as Blue and Size as 2
DATASET #1 Mean: Color as Red and Size as 2

now click “Done”.

11) change the axis limits by clicking “Set Axis” at the bottom of the application; then choose “Custom Scale” from the top of the popup window and set “Min X” = 1.0E-2, “Max X” = 3.0, and “Min Y” = 1.0E-6. Click “OK” to apply the changes and close the window.

The window you now see represents the screenshot shown in Figure 2.


FIGURE 3:

(this assumes you've gone through steps 1-11 above)

12) Put the present plot aside by clicking “Peel Off” at the bottom of the application (or click “Clear Plot”).

13) In the main application, make it so we’ll only plot the mean curves by doing the following: in the “Control Panels” list at the bottom of the application, choose “Epistemic List Control” (you might need to scroll down to see it). Deselect “Plot all curves”, and choose “No Fractiles”. Only the “Plot Average” option should remain selected.

14) Re-compute the average shown in Figure 2 (for Campbell and Bazorgnia (2003) by again clicking “Compute” at the bottom of the application.

15) Now redo this calculation using a different IMR (attenuation relationship) by doing the following: click the “IMR, IMT & Site” tab at the top of the application (if you haven’t already done this), under the “Set IMR” window choose “Boore et al. (1997)” from the list at the top, set the IMT to “PGA in the “Set IMT” window below, and click “Compute” at the bottom. Note that the site-type parameter under the lat and lon was already set in step (4) above.

16) Repeat step (15) for the “Abrahamson and Silva (1997)” and “Sadigh et al. (1997)” attenuation relationships in that order (and be sure to set the IMT to PGA each time).

17) After increasing the size of the lines to 2.0 (by clicking “Plot Prefs”), and changing the axis range (by clicking “Set Axis”), you have now reproduced Figure 3. Note that you can save a png version of the image by selecting the "Save" option in the "File" menu at the upper left corner of the application or by clicking on the save icon symbol at the top left of the application.


FIGURE 4:

(this assumes you gone through steps 1-17 above)

18) Click either “Peel Off” or “Clear Plot” at the bottom of the application.

19) In the main application, re-select the Campbell & Bozorgnia (2003) IMR (attenuation relationship), set the IMT to PGA, and click “Compute” at the bottom of the application. This re-plots the mean curve shown in Figure 2.

20) Click the “ERF & Time Span” tab at the top of the application, and on the right side change the time-span “Duration” from 30 to 20 years. Click “Compute” at the bottom of the application.

21) Repeat step (20) for durations of 10, 5, and 1 (in that order).

22) After increasing the size of the lines to 2.0 (by clicking “Plot Prefs”), and changing the axis range (by clicking “Set Axis”), you have now reproduced Figure 4.


FIGURE 5:

(this assumes you gone through steps 1-22 above)

23) Click either “Peel Off” or “Clear Plot” at the bottom of the application.

24) In the main application, re-select the duration of 30 years and click “Compute” (to re-plot the mean curve shown in Figure 2).

25) For each of the 7 ERF adjustable parameters labeled “[fault name] Prob. Model Wts”, click the button and set all values in the popup window to zero, except set that for the “Poisson” model to 1.0. After dong this for the 7 parameters (one for each fault), click “Compute” at the bottom of the application.

26) Repeat step (25), but rather than giving the “Poisson” model a weight of 1.0, give the “BPT”, model a weight of 1.0 in each of the 7 parameters.

27) After increasing the size of the lines to 2.0 (by clicking “Plot Prefs”), and changing the axis range (by clicking “Set Axis”), you have now reproduced Figure 5.