Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Ground motion model (GMM) response spectra plots
%
% Author: Peter Powers (pmpowers@usgs.gov)
%
% Updated: 01/19/2023
%
% This script demonstrates how to call web services for the ground motion
% models implemented in USGS National Sesmic Hazard Models and plot the
% results. A JSON web service response is automatically converted to a
% struct by the Matlab webread() function. This script shows how to plot
% the GMMs requested and optionally the underlying models of epistemic
% uncertainty.
%
% Users are STRONGLY encouraged to run the web services locally so as to
% not overburden USGS servers. See the repository README for details on how
% to do so: https://code.usgs.gov/ghsc/nshmp/nshmp-ws
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Plot a figure of GMM reponse spectra and epistemic uncertainty
clearvars
% The root url for the response spectra web service
urlbase = "https://earthquake.usgs.gov/ws/nshmp/data/gmm/spectra?";
% If you are running nshmp-haz locally, use this URL:
% urlbase = "http://localhost:8080/gmm/spectra?";
% Struct of ground motion model parameters
in.Mw = 7.5;
in.dip = 90;
in.rake = 0;
in.width = 15;
in.rJB = 50;
in.rRup = 50;
in.rX = 50;
in.vs30 = 760;
in.zHyp = 7.5;
in.zTor = 10;
figure
plot_handles = [];
% Set to true to show epistemic branches
epi = false;
gmms_2023 = ["ASK_14" "BSSA_14" "CB_14" "CY_14"];
means = readSpectra(urlbase, gmms_2023, in).response.means.data;
plotGmms(means, in, epi);
%% Shared Functions
function plotGmms(gmm_data, in, epi)
plot_handles = [];
legend_labels = {};
% loop means_new response array
for j = 1 : length(gmm_data)
gmm_xs = gmm_data(j).data.xs;
gmm_ys = gmm_data(j).data.ys;
gmm_tree = gmm_data(j).tree;
% if 0.01 absent use PGA (0.001 is used for PGA)
if gmm_xs(1) == 0.001 && gmm_xs(2) ~= 0.01
gmm_xs(1) = 0.01;
end
% Plot the total spectrum
ph = loglog(gmm_xs, gmm_ys, ...
'LineWidth',3, ...
'LineStyle','-');
color = get(ph,'Color');
hold on;
plot_handles = [plot_handles ph];
legend_labels{end+1} = gmm_data(j).label;
% Plot epistemic spectra, if present
if epi && ~isempty(gmm_tree)
for k = 1 : length(gmm_tree)
id = gmm_tree(k).id;
epi_ys = gmm_tree(k).values;
phe = loglog(gmm_xs, epi_ys, ...
'LineWidth',1, ...
'LineStyle',':', ...
'Color',color);
plot_handles = [plot_handles phe];
legend_labels{end+1} = [gmm_data(j).label ' ' id];
end
end
end
grid on
xlabel('Spectral Period (s)','FontSize',16)
ylabel('Median Ground Motion (g)','FontSize',16)
title(['M',num2str(in.Mw),', ',num2str(in.rRup),'km'],'FontSize',20)
xlim([0.01 10]);
set(gca,'FontSize',20);
set(gca,'XTick',[0.01 0.1 1 10]);
set(gca,'XTickLabel',{'0.01','0.1','1','10'});
legend( ...
plot_handles,legend_labels, ...
'Location','southwest', ...
'Interpreter', 'none');
end
% Read reponse spectra from a NSHM web service for the supplied
% ground motion models (GMMs) and GMM input parameters
function response = readSpectra(urlbase, gmms, in)
url = urlbase;
for i = 1 : size(gmms, 2)
if i == 1
url = url + "gmm=" + gmms(i);
else
url = url + "&gmm=" + gmms(i);
end
end
url = url + ...
"&Mw=" + num2str(in.Mw) + ...
"&dip=" + num2str(in.dip) + ...
"&rake=" + num2str(in.rake) + ...
"&width=" + num2str(in.width) + ...
"&rJB=" + num2str(in.rJB) + ...
"&rRup=" + num2str(in.rRup) + ...
"&rX=" + num2str(in.rX) + ...
"&vs30=" + num2str(in.vs30) + ...
"&zHyp=" + num2str(in.zHyp) + ...
"&zTop=" + num2str(in.zTor);
if isfield(in, 'z2p5')
url = url + num2str(in.z2p5);
end
if isfield(in, 'z1p0')
url = url + num2str(in.z1p0);
end
if isfield(in, 'zSed')
url = url + num2str(in.zSed);
end
response = webread(url);
end