Auxiliary material for Paper 2009TC002565
Exploratory models of long-term crustal flow and resulting seismicity across the Alpine-Aegean orogen.
Tracy M. Howe and Peter Bird
(Department of Earth and Space Sciences, University of California, Los Angeles, CA 90095-1567)
Tectonics, doi: 10.1029/2009TC002565.

Introduction
To encourage experimentation and improvement of the model, 
we attach a table with data sources for faults, 
and also all files (not previously published) that are needed to reproduce 
the preferred "creeping subduction" model #200921.
All files have filename extension ".txt" and are plain ASCII text.  
For source codes of NeoKinema and Long_Term_Seismicity, 
see the electronic supplements attached to Bird (2009, 
J. Geophys. Res., 114, B11403, doi:10.1029/2009JB006317).  
For grid-editing and map-plotting utility software
(e.g., programs OrbWin and NeoKineMap, respectively), 
see http://peterbird.name/oldFTP/.

1. f_Alps_2009_01_21.dig.txt
Traces of active and potentially-active faults, in (longitude, latitude) coordinates,
with East longitude and North latitude positive.
The letter(s) "D, L, N, P, R, S, T" following the fault index number (e.g., F0001)
indicate the nominal sense(s) of slip for plotting purposes (as in Figure 1),
but these senses of slip are occasionally reversed in NeoKinema solutions.
If dip is not vertical, then dip is to the left when traveling along the trace
in the direction of digitization.
Note that most traces were hand-digitized from small-scale regional maps,
and these locations do not have sufficient precision for local mapping purposes.
The format of .dig files is further documented at:
http://peterbird.name/guide/dig_format.htm

2. f_Alps_nki_2009_04_16.txt
Table with data sources for faults, saved as tab-delimited ASCII text.
Column headers and contents are as follows:
2.01 FAULT = letter "F" followed by a 4-digit NeoKinema serial number (arbitrary, and not always sequential).
2.02 OFFSET = sense of slip, represented by "D, L, N, P, R, S, T" as explained in text.
            (Note that faults with oblique slip will have two entries, on adjacent lines; e.g., F0009N & F0009R.)
2.03 Descriptive text = name of fault (if known to us), and general location.
2.04 Offset-rate = long-term-average offset rate (if known) in mm/a.  Otherwise, zero.
2.05 Sigma = uncertainty (standard error) of offset rate (if known) in mm/a.  Otherwise, a large number, as discussed in text.
2.06 Creeping? = logical indicator, set "T" if fault creeps steadily without producing earthquakes.
2.07 CULxKm = upper locking depth in km, if known; otherwise "-1.0".
2.08 LLxKm = lower locking depth in km, if known; otherwise "-1.0".
2.09 Data? = TRUE if offset rate is based on offset geologic feature(s) we found in literature; usually FALSE.
2.10 Value? = offset rate (range), as found in literature, in free-format, in mm/a.  Only entered if Data? = TRUE.
2.11 Source(s) = short citation of literature source; full citation is found in main paper.
2.12 Note(s) = any additional comments, in free-format text.

3. f_Alps_2009_04_16.nki.txt
Fault offset-rates and uncertainties, in the input format required by NeoKinema.
Obtained by deleting unneeded columns from the table above, saving it as 
space-delimited text, and adding a FORTRAN FORMAT line as a header.
For meaning of columns, see the paragraph above in this file.

4. s_Alps.nki.txt
Azimuths (in degrees clockwise from local North) of most-compressive horizontal principal stress axes;
obtained from the World Stress Map and converted to the input format required by NeoKinema.
Note that line 1 is a FORTRAN FORMAT needed to read the rest of the file, and
that line 2 contains text identifying the columns:
4.1 text1 = site identifiers, copied from WSM data file.
4.2 text2 = further subsite identifier (if any), copied from WSM data file.
4.3 E_long. = East longitude of stress measurement site, in decimal degrees.
4.4 N_lat. = North latitude of stress measurement site, in decimal degrees.
4.5 azi = azimuth (in degrees clockwise from local North) of most-compressive horizontal principal stress axis.
4.6 quality = A, B, C, D, E (as assigned by WSM).

5. g_Alps_merged.gps.txt
Velocities of geodetic benchmarks, with uncertainties.
Obtained from a variety of published sources (as explained and listed in text, 
where full citations may be found).
Converted to the .gps format documented at http://peterbird.name/guide/gps_format.htm
and required for input to NeoKinema.
Note that all velocities are expressed in the velocity reference frame of stable Eurasia.

6. Alps8C.feg.txt
Finite-element grid used in the NeoKinema simulation.
Created with utility programs OrbWin and OrbNumber, available at:
http://peterbird.name/oldFTP/

7. b_Alps8C.nki.txt
Boundary conditions for the finite-element grid listed above.
Note that most boundary nodes are assigned type-4 boundary conditions,
meaning that only the identity of the adjacent rigid plate is specified
(EU = Eurasia, AF = Africa, AT = Anatolia) per Bird (2003).
Program NeoKinema contains Euler poles necessary to compute these velocities.

8. p_Alps200921.nki.txt
Parameter input file for the NeoKinema run which produced preferred model Alps200921.
Meaning of each entry is briefly explained by a text comment on the right side of the line.
Meanings are more fully explained in the NeoKinema "manual" which is built into the source code.