Recipe of plotting the figures

Notebooks and codes to reproduce the figures

This page shows the list of the link to the notebooks and codes to plot the figures. We cleaned up the output of the cells showing the figures in the notebooks due to minimizing the size of repositories and the copyrights. To replot them, please download the data and follow the instructions documented in the directories.

Figure 1. Map of Parkfield. Download and generate the input files such as .grd files following Maps/BPnetwork_introduction/README.md. Then, run plot_map_introduction_Parkfield.sh

Figure 2. Number of station-channel pairs. Run Post/ModelFit/code/plotfigure_plotdvvpsd.ipynb. See Post/ModelFit/README.md for the detail.

Figure 3. All dv/v time history. Run Post/ModelFit/code/plotfigure_plotdvvpsd.ipynb.

Figure 4. Dv/v with different station-channel pairs. Run Post/ModelFit/code/plotfigure_plotdvvpsd.ipynb.

Figure 5. Dv/v with different frequency bands. Run Post/ModelFit/code/plotfigure_plotdvvpsd.ipynb.

Figure 6. Schematic of dv/v model components. Run Appx/modelschematic/plot_model_schematics_dvv.ipynb to generate the time histories of model components. The data is imported from the Post/ModelFit/data/interped_tempandprecip_longterm.csv.

Figure 7. Scatter matrix of MCMC using stretching associated with LCCB-SCYB. Run Post/ModelFit/code/plotfigure_MCMCscattermatrix_stretching.ipynb. You need to download the output of the MCMC sampler. See the details in Post/ModelFit/README.md.

Figure 8. Scatter matrix of MCMC using MWCS associated with LCCB-SCYB. Run Post/ModelFit/code/plotfigure_MCMCscattermatrix_mwcs.ipynb.

Figure 9. Fitting of dv/v with stretching. Run Post/ModelFit/code/plotfigure_MCMCdvvmodelfit.ipynb. You can select the dvvmethod as stretching or MWCS.

Figure 10. Fitting of dv/v with MWCS. Run Post/ModelFit/code/plotfigure_MCMCdvvmodelfit.ipynb.

Figure 11. Stats plot of model parameters. Run Post/ModelFit/code/modelfit_07_MCMC_plotstats.ipynb.

Figure 12. Model fit with the residual healing model. Run Appx/casestudy_residual_healing/code/resheal_plot_MCMCdvvmodelfit.ipynb.

Figure 13. Map of cumulative dilation and shear strains. Run Post/CumulativeStrain/code/02_postprocess_cumulativestrain_prepareGMTinput.ipynb to dump the data of cumulative strain associated with the triangular elements. Them, run /Post/CumulativeStrain/code/plotmap/plot_gps_cumulativestrain_v4.sh to plot the maps with the strain.

Figure 14. Sensitivity of dv/v to the dilation and shear strains. Run Post/CumulativeStrain/code/06_plotcomparison_dvvandstrain.ipynb.

Figure 15. Rotated axial strain and its sensitivity. Run Post/CumulativeStrain/code/06_plotcomparison_dvvandstrain.ipynb.

Figure S1. Map with the station names and the approximated planar fault. Run Maps/BPnetwork_Faultdist/BPnetwork_GMT_local/plot_BPnetwork_GMTlocal.sh.

Figure S2. PPSD associated with the EADB. Followng Post/Spectrogram/README.md. The dataset of PPSD is available in SeisMonitoring_PPSDdata.tar.gz, and run Post/Spectrogram/code/plot_PPSD_BP_3cols.ipynb to plot the figures.

Figure S3. Removal of earthquakes and tremors. Download the data and the remove the transient signals by running the SeisMonitoring processing in Appx/removal_eqandtremor/seisremoveeq_demo. Then run Appx/removal_eqandtremor/plot_removalschematic/Tremor_signal_convert_Juliadata.ipynb and Appx/removal_eqandtremor/plot_removalschematic/Tremor_signal_plot_removal.ipynb to plot the example of the removal of earthquakes and tremors from the raw data.

Figure S4. Data availability. The datasheet of data availability is obtained by the Post/DataAvailability/compute_DataAvailability.jl in the cluster. Then, run Post/DataAvailability/plot_DataAvailability_BP.ipynb to plot the figures.

Figure S5. Cross-correlation functions associated with the LCCB-SCYB. Run Appx/plot_CCF/code/plot_ccf_master_v04_medianmute.ipynb to plot the CCFs over 20 years for the 9 components associated with a given station pair. As the data size is large, we plotted the all station pairs using Appx/plot_CCF/code/plot_ccf_master_v04_medianmute.py in the work stations. In the notebook, you can download an example of the set of CCF data from here. NOTE: To recreate the .npz file, we can also download the cross-correlation file with .jld2 from dasway and locate it in Appx/plot_CCF/cc_channel_collection/. Then, run convert_ccf_threads_uwcascadia.jl to convert the data to the .npz format. We can plot the CFs using the dataset in the .npz format.
Figure S6 Comparison of dv/v and LFE activity. Run Appx/comparison_ParkfieldLFErate/code/comparison_ParkfieldLFErate.ipynb.
Figure S7. Dv/v time histories with channel-weighted station pairs for the case with the stretching. Run Post/ModelFit/code/plotfigure_alldvvmodelfit.ipynb.

Figure S8. Dv/v time histories with channel-weighted station pairs for the case with the MWCS. Run Post/ModelFit/code/plotfigure_alldvvmodelfit.ipynb.

Figure S9. The long-term increase in dv/v with different frequency bands. Run Post/ModelFit/code/plotfigure_plotdvvpsd.ipynb.

Figure S10. Sensitivity kernel with depth. Run Appx/BP_sensitivity_kernel/plot_sensitivity_kernel_Parkfield.ipynb. See the details in Appx/BP_sensitivity_kernel/README.md.

Figure S11. Trade-off between S, τmin and τmax. Run Appx/tradeoff_logheal_SandTmax/code/tradeoff_logheal_SandTminTmax_v2.ipynb to conduct the MCMC parameter search associated with the S, taumin and taumax. Then, run tradeoff_logheal_SandTmax/code/tradeoff_logheal_SandTminTmax_plotmaster.ipynb to plot the summary figure.

Figure S12. Scatter matrix with MWCS, residual healing model. Run Appx/casestudy_residual_healing/code/resheal_plot_MCMCscattermatrix_mwcs.ipynb

Figure S13. The coefficient of $b_0$ with the fault-normal distance. Run Post/ModelFit/code/modelfit_07_MCMC_plotstats.ipynb.

Figure S14. The scaling of parallelization. Run Appx/Scaling_Frontera/plot_scaling_parallelization_master.ipynb.
Figure S15. The spatial distribution of dilatational strain at Parkfield. . Run Others/coseismic_strain/plot_strain_gmt/plot_BP_dilation.sh and Others/coseismic_strain/compute_strain_BP/code/02_plot_dilation_cross-section_BP.ipynb to plot the horizontal distribution and cross-section of the strain. See the details in Others/coseismic_strain/README.md.
Table S1. Run Others/interstation_distance/compute_stationdistances.ipynb.