Auxiliary material for Paper 2009JB006317 Long-term fault slip rates, distributed deformation rates, and forecast of seismicity in the western United States from joint fitting of community geologic, geodetic, and stress-direction datasets Peter Bird Department of Earth and Space Sciences, University of California, Los Angeles, California, USA Bird, P. (2009), Long-term fault slip rates, distributed deformation rates, and forecast of seismicity in the western United States from joint fitting of community geologic, geodetic, and stress-direction datasets, J. Geophys. Res., 114, XXXXXX, doi:10.1029/2009JB006317. Introduction This auxiliary material is a set of 15 ASCII computer files, including this one. It provides the essential files, and instructions, to permit reproduction of the preferred kinematic finite-element model GCN2008088 which was presented in the text, tables, and figures of the paper. It also includes the gridded long-term seismicity matrix predicted by that model, a table of ranges of fault offset-rate predicted by all 24 acceptable models, and 2 files permitting the use of GCN2008088 as a "deformation model" in alternative unified rupture forecasts. The necessary steps to reproduce GCN2008088 would include: (1) Strip the second filename extension ".txt" from all files provided. (2) Compile NeoKinema_v2p2_Guadalupe.f90 with a Fortran90 compiler. -Note that the source statement: USE Sphere refers to another Fortran90 module by Peter Bird, named Sphere.f90. This is available at: http://peterbird.name/oldFTP/MapTools/ and must also be compiled and linked into the executable. -Note that the source statement: USE Dislocation refers to another Fortran90 module by Peter Bird, named Dislocation.f90. This is available at: http://peterbird.name/oldFTP/NeoKinema/ and must also be compiled and linked into the executable. -Note that the source statement: USE Numerical_Libraries is used to link to subprograms LSLPB, LSLSF, EVCRG, & LINDS, which are parts of I.M.S.L. (International Mathematics Subroutine Library); my code is customized to work with IMSL's storage conventions for matrices representing the coefficients of linear systems. If different solvers are substituted then NeoKinema will need to be modified. -Note that the source statement: USE DFLIB, ARCQQ => ARC ! provided with Digital Visual Fortran allows reference to subprogram GETFILEINFOQQ, which provides names of files matching specifications like "p*.nki". This helps user to select input file. If no substitute is available on your system when you compile, just omit SUBROUTINE File_List (and any CALLs to it). [N.B. If you wish to run NeoKinema under 32-bit Windows, you may not need to compile it yourself. Obtain NeoKinema_v2p2_Guadalupe.exe from: http://peterbird.name/oldFTP/NeoKinema/ ] (3) Assemble the necessary input data files in a folder/directory which will be the current/active directory for NeoKinema_v2p2_Guadalupe.exe: p_GCN2008088.nki fGCN_merged_WGCEPFM2p2_200810.nki fGCN_merged_WGCEPFM2p2_200810.dig s_Gorda-Cal-Nev.nki GPS2006_NA_over2kmMergedFM2p2.gps GPS2006_NA_over2kmMergedFM2p2.gp2 GCN8p9.feg b_GCN8p9.nki (4) Start NeoKinema_v2p2_Guadalupe.exe and answer the (only) prompt by specifying parameter file name p_GCN2008088.nki. (At this point execution may be left unattended. It may take ~3 days.) (5) Use program NeoKineMap (from http://peterbird.name/oldFTP/NeoKinema/ ) to make any desired plots of input and/or output fields. (Note that compiling NeoKineMap requires a Fortran90 compiler, plus modules Adobe_Illustrator.f90, Map_Projections.f90, Map_Tools.f90, and Icosahedron.f90; all are available at: http://peterbird.name/oldFTP/Map_Tools/ .) 1. 2007jb005460-SHELLS.f90.txt This is the Fortran90 source code for thin-shell dynamic neotectonic finite-element program SHELLS. Version of 29 August 2006, including important new features added for this paper. 2. 2007jb005460-Earth5R.feg.txt This is a finite-element grid representing the global lithosphere, including its faults which divide it into plates (and subdivide some plates into blocks within orogens). For each node, the values given in SI units are: longitude (degrees East), latitude (degrees North), elevation, heat-flow, crustal thickness, mantle lithosphere thickness, density anomaly of compositional origin, & geothermal curvature due to transient cooling. To edit this grid (or simply to visualize it) you will want to use our program OrbWin, from: http://peterbird.name/oldFTP/OrbWin/ 3. iEarth5-046.in.txt This is an input-parameter file, corresponding to Table 1 in the paper. (Also note that Table 2 of the paper shows how we determined optimal parameters.) This is the first of a group of 4 input files needed to reproduce Earth5-049. Most lines in this group of files do not change; however, note changes in: -line 1: title line (changed to reflect number of this run: ..."-046") -line 14: value of parameter ICONVE (= 3, but irrelevant due to following...) -line 15: value of parameter TRHMAX (=0, preventing any distributed basal traction) 4. iEarth5-047.in.txt This is an input-parameter file, corresponding to Table 1 in the paper. (Also note that Table 2 of the paper shows how we determined optimal parameters.) This is the second of a group of 4 input files needed to reproduce Earth5-049. Most lines in this group of files do not change; however, note changes in: -line 1: title line (changed to reflect number of this run: ..."-047") -line 14: value of parameter ICONVE (= 6; use inferred basal tractions) -line 15: value of parameter TRHMAX (= 20.E6; limit raised to 20 MPa) 5. iEarth5-048.in.txt This is an input-parameter file, corresponding to Table 1 in the paper. (Also note that Table 2 of the paper shows how we determined optimal parameters.) This is the third of a group of 4 input files needed to reproduce Earth5-049. Most lines in this group of files do not change; however, note changes in: -line 1: title line (changed to reflect number of this run: ..."-048") -line 14: value of parameter ICONVE (= 6; use inferred basal tractions) -line 15: value of parameter TRHMAX (= 20.E6; limit raised to 20 MPa) 6. iEarth5-049.in.txt This is an input-parameter file, corresponding to Table 1 in the paper. (Also note that Table 2 of the paper shows how we determined optimal parameters.) This is the last of a group of 4 input files needed to reproduce Earth5-049. Most lines in this group of files do not change; however, note changes in: -line 1: title line (changed to reflect number of this run: ..."-049") -line 14: value of parameter ICONVE (= 6; use inferred basal tractions) -line 15: value of parameter TRHMAX (= 20.E6; limit raised to 20 MPa) 7. Earth5R-type4AplusA.bcs.txt This is the first of two versions of the boundary-conditions file appropriate for use with finite element grid Earth5R.feg. This version includes "extra" velocity boundary conditions in plate interiors, and is used for steps (4), (5), (6). 8. Earth5R-type4A.bcs.txt This is the second of two versions of the boundary-conditions file appropriate for use with finite element grid Earth5R.feg. This version does not include any "extra" velocity boundary conditions, and is used only in the last step (7). 9. 2007jb005460-qEarth5-048.out.txt This is a "torque report" file with content similar to Table 4 (and Figures 15-16) in the paper. It will be needed if you jump directly to step (7), because normally it would be produced by running step (6).