

************************************************************************************************************
**************************************** Introduction ******************************************************
************************************************************************************************************

This repository contains the model input mesh file and model data output files for the figures in the main body of the paper "Bridging Paleoseismic and Geodetic Observations for the AD 1700 Cascadia Earthquake with an Earthquake Cycle Model", by Weilun Qin, Rob Govers, Natasha L.M. Barlow, and Riccardo E.M. Riva submitted for publication in Geophysical Research Letters.



************************************************************************************************************
**************************************** File Content ******************************************************
************************************************************************************************************

- [Cascadia_EQC_300_20to30_80.msh]		
3D meshes of the Cascadia subduction earthquake cycle model. It has a model depth of 300 km with the seismogenic zone varies from 20 to 30 km and shear zone terminates at 80 km. Recommend to view using Gmsh.


- [vel_surf_tot_3180_GPSsites_xy.dat]	
Horizontal velocity output 310 years after the AD 1700 Cascadia earthquake (i.e., step 3180) at GPS sites in McKenzie & Furlong (2021)
File format: Longitude | Latitude | Azimuth (in degrees east of north) | Horizontal velocity (in mm/year)


- [residual_vel_surf_tot_3180_GPSsites_xy.dat]	
Residual between horizontal velocity output 310 years after the AD 1700 Cascadia earthquake compared to GPS velocities used in McKenzie & Furlong (2021)
File format: Longitude | Latitude | Azimuth (in degrees east of north) | Horizontal velocity Residual (in mm/year)


- [dsp_surf_tot_3039_3049_z_deg.dat]	
Coseismic vertical surface deformation for the AD 1700 Cascadia earthquake (i.e., difference between step 3049 and 3039)
File format: Longitude | Latitude | Vertical displacement (in m)


- [dsp_surf_tot_3039_3059(50%)_z_deg.dat]	
Coseismic + 50% Primary afterslip surface deformation for the AD 1700 Cascadia earthquake (i.e., 50% difference between step 3059 and 3049 plus the coseismic deformation)
File format: Longitude | Latitude | Vertical displacement (in m)


- [conn_surf.dat]				
Connectivity file for the displacement output


