diff --git a/src/gov/usgs/earthquake/nshmp/gmm/GraizerKalkan_2015.java b/src/gov/usgs/earthquake/nshmp/gmm/GraizerKalkan_2015.java
index df90cd160260577a6c68fc758cf37835aa0ad710..3553b3f98cf0db7ab7de2387b426913848e7283c 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/GraizerKalkan_2015.java
+++ b/src/gov/usgs/earthquake/nshmp/gmm/GraizerKalkan_2015.java
@@ -62,7 +62,7 @@ public final class GraizerKalkan_2015 implements GroundMotionModel {
*
* Basin term is a function of z1p5. Current implementation linearly
* interpolates between z1p0 and z2p5 if they are not NaN, otherwise uses a
- * default value of 0.15 km. This value is
+ * default value of 0.15 km.
*/
static final String NAME = "Graizer & Kalkan (2015)";
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/GroundMotionTables.java b/src/gov/usgs/earthquake/nshmp/gmm/GroundMotionTables.java
index 185b784d2101e1f99fcba1070993a3fb0ae0ecd7..8a8ddc2f3163b16156ec4203b0281ffc3a1685ed 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/GroundMotionTables.java
+++ b/src/gov/usgs/earthquake/nshmp/gmm/GroundMotionTables.java
@@ -46,8 +46,8 @@ import gov.usgs.earthquake.nshmp.internal.Parsing.Delimiter;
*
* Frankel, Atkinson, and Pezeshk tables store ground motion in log10 values.
* Atkinson flavored tables store ground motion in cm/s^2. NGA-East tables
- * contain linear ground motion values. All tables interpolate in log10 distance
- * and linear in magnitude.
+ * contain linear ground motion values but are converted to natural-log when
+ * parsed. All tables interpolate in log10 distance and linear in magnitude.
*
* @author Peter Powers
*/
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java b/src/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java
index 41f4964673cced2aee13b2820711e0f0539f76b3..6cefc0533acc62eb504160d6ff0e81799e74cd3b 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java
+++ b/src/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java
@@ -3,6 +3,9 @@ package gov.usgs.earthquake.nshmp.gmm;
import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.MW;
import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.RJB;
import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.VS30;
+import static java.lang.Math.exp;
+import static java.lang.Math.log;
+import static java.lang.Math.min;
import com.google.common.annotations.Beta;
import com.google.common.collect.Range;
@@ -35,10 +38,16 @@ import gov.usgs.earthquake.nshmp.util.Maths;
* properly processed by {@link ExceedanceModel#NSHM_CEUS_MAX_INTENSITY} at this
* time.
*
- * <p> This class also manages implementations of the 19 'seed' models used to
+ * <p>This class also manages implementations of the 19 'seed' models used to
* generate (via Sammons mapping) the 13 NGA-East models and associated weights.
* Ground motions for the 19 seed and 13 component models are computed via table
* lookups.
+ *
+ * <p>On it's own, NGA-East is a hard rock model returning results for a site
+ * class where Vs30 = 3000 m/s. To accomodate other site classes, the Stewart et
+ * al. (2017) CEUS site amplification model is used. This model is applicable to
+ * 200 ≤ vs30 ≤ 2000 m/s. In the current implementation, for vs30 < 200, vs30 =
+ * 200 m/s; for vs30 > 2000, vs30 = 3000 m/s.
*
* <p><b>Note:</b> Direct instantiation of {@code GroundMotionModel}s is
* prohibited. Use {@link Gmm#instance(Imt)} to retrieve an instance for a
@@ -51,7 +60,16 @@ import gov.usgs.earthquake.nshmp.util.Maths;
* Youngs, R., Graves, R., Atkinson, G., 2017, NGA-East ground-motion models for
* the U.S. Geological Survey national seismic hazard maps: PEER Report No.
* 2017/03, 180 p.
+ *
+ * <p><b>Reference:</b> Stewart, J., Parker, G., Harmon, J., Atkinson, G.,
+ * Boore, D., Darragh, R., Silva, W., and Hashash, Y., 2017, Expert panel
+ * recommendations for ergodic site amplification in central and eastern North
+ * America: PEER Report No. 2017/04, 66 p.
*
+ * <p><b>Reference:</b> Hashash, Y., Harmon, J., Ilhan, O., Parker, G., and
+ * Stewart, 2017, Recommendation for ergodic nonlinear site amplification in
+ * central and eastern North America: PEER Report No. 2017/05, 62 p.
+ *
* <p><b>Component:</b> average horizontal (RotD50)
*
* @author Peter Powers
@@ -69,9 +87,9 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
* least requires more study to determine if it generates an acceptable
* approximation
*
- * Note: the following inconsistency exists: NgaEast GroundMotionTables
- * include SA0P025 but sigma table do not. Sigma tables are used to initialize
- * supported Imts so there is never a problem.
+ * Note: supported periods are derived from sigma coefficient files. Several
+ * supported periods have been commented out because they are not represented
+ * in teh site amplification model.
*
* When supplied with tables for the 13 usgs models, 0.01s was added with
* values distinct from PGA, however the seed models are missing this period.
@@ -87,7 +105,7 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
static final Constraints CONSTRAINTS = Constraints.builder()
.set(MW, Range.closed(4.0, 8.2))
.set(RJB, Range.closed(0.0, 1500.0))
- .set(VS30, Range.singleton(3000.0))
+ .set(VS30, Range.closed(150.0, 3000.0))
.build();
/*
@@ -147,7 +165,8 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
private final Coefficients σCoeffsHi;
private final CoefficientsTotal σCoeffsTotal;
- private final GroundMotionTable[] tables;
+ // private final GroundMotionTable[] tables;
+ // private final GroundMotionTable[] pgaTables;
private final double[] weights;
NgaEastUsgs_2017(final Imt imt) {
@@ -155,7 +174,7 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
σCoeffsMid = new Coefficients(imt, COEFFS_SIGMA_MID);
σCoeffsHi = new Coefficients(imt, COEFFS_SIGMA_HI);
σCoeffsTotal = new CoefficientsTotal(imt, COEFFS_SIGMA_TOTAL);
- tables = GroundMotionTables.getNgaEast(imt);
+ // tables = GroundMotionTables.getNgaEast(imt);
weights = GroundMotionTables.getNgaEastWeights(imt);
}
@@ -226,20 +245,30 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
final int[] models;
final double[] weights;
+ final GroundMotionTable[] tables;
+ final GroundMotionTable[] pgaTables;
+ final SiteAmp siteAmp;
/* Specifiy an array of models ids. */
ModelGroup(Imt imt, int[] models) {
super(imt);
this.models = models;
this.weights = Data.round(8, Data.normalize(selectWeights(super.weights, models)));
+ this.tables = GroundMotionTables.getNgaEast(imt);
+ this.pgaTables = GroundMotionTables.getNgaEast(Imt.PGA);
+ this.siteAmp = new SiteAmp(imt);
}
@Override
public MultiScalarGroundMotion calc(GmmInput in) {
- Position p = super.tables[0].position(in.rRup, in.Mw);
+ Position p = tables[0].position(in.rRup, in.Mw);
double[] μs = new double[models.length];
for (int i = 0; i < models.length; i++) {
- μs[i] = super.tables[models[i] - 1].get(p);
+ int ti = models[i] - 1;
+ double μ = tables[ti].get(p);
+ double μPGA = exp(pgaTables[ti].get(p));
+ double fSite = siteAmp.calc(μPGA, in.vs30);
+ μs[i] = μ + fSite;
}
double[] σs = calcSigmas(in.Mw);
double[] σWts = σs.length > 1 ? SIGMA_WTS : new double[] { 1.0 };
@@ -273,19 +302,25 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
final int id;
final GroundMotionTable table;
+ final GroundMotionTable pgaTable;
+ final SiteAmp siteAmp;
Sammons(int id, Imt imt) {
super(imt);
this.id = id;
- this.table = super.tables[id - 1];
+ this.table = GroundMotionTables.getNgaEast(imt)[id - 1];
+ this.pgaTable = GroundMotionTables.getNgaEast(Imt.PGA)[id - 1];
+ this.siteAmp = new SiteAmp(imt);
}
@Override
public ScalarGroundMotion calc(GmmInput in) {
Position p = table.position(in.rRup, in.Mw);
- return new DefaultScalarGroundMotion(
- table.get(p),
- calcSigmaTotal(in.Mw));
+ double μPGA = exp(pgaTable.get(p));
+ double fSite = siteAmp.calc(μPGA, in.vs30);
+ double μ = table.get(p) + fSite;
+ double σ = calcSigmaTotal(in.Mw);
+ return new DefaultScalarGroundMotion(μ, σ);
}
}
@@ -411,19 +446,25 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
final String id;
final GroundMotionTable table;
+ final GroundMotionTable pgaTable;
+ final SiteAmp siteAmp;
Seed(String id, Imt imt) {
super(imt);
this.id = id;
this.table = GroundMotionTables.getNgaEastSeed(id, imt);
+ this.pgaTable = GroundMotionTables.getNgaEastSeed(id, Imt.PGA);
+ this.siteAmp = new SiteAmp(imt);
}
@Override
public ScalarGroundMotion calc(GmmInput in) {
Position p = table.position(in.rRup, in.Mw);
- return new DefaultScalarGroundMotion(
- table.get(p),
- calcSigmaTotal(in.Mw));
+ double μPGA = exp(pgaTable.get(p));
+ double fSite = siteAmp.calc(μPGA, in.vs30);
+ double μ = table.get(p) + fSite;
+ double σ = calcSigmaTotal(in.Mw);
+ return new DefaultScalarGroundMotion(μ, σ);
}
}
@@ -606,31 +647,18 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
super(ID, imt);
}
}
-
- static abstract class UpdatedSeed extends NgaEastUsgs_2017 {
- static final String NAME = NgaEastUsgs_2017.NAME + ": Seed (updated) : ";
- final String id;
- final GroundMotionTable table;
+ static abstract class UpdatedSeed extends Seed {
+ static final String NAME = NgaEastUsgs_2017.NAME + ": Seed (updated) : ";
UpdatedSeed(String id, Imt imt) {
- super(imt);
- this.id = id;
- this.table = GroundMotionTables.getNgaEastSeed(id, imt);
- }
-
- @Override
- public ScalarGroundMotion calc(GmmInput in) {
- Position p = table.position(in.rRup, in.Mw);
- return new DefaultScalarGroundMotion(
- table.get(p),
- calcSigmaTotal(in.Mw));
+ super(id, imt);
}
}
static final class UpdatedSeed_Graizer16 extends UpdatedSeed {
static final String ID = "Graizer16";
- static final String NAME = Seed.NAME + ID;
+ static final String NAME = UpdatedSeed.NAME + ID;
UpdatedSeed_Graizer16(Imt imt) {
super(ID, imt);
@@ -639,12 +667,167 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
static final class UpdatedSeed_Graizer17 extends UpdatedSeed {
static final String ID = "Graizer17";
- static final String NAME = Seed.NAME + ID;
+ static final String NAME = UpdatedSeed.NAME + ID;
UpdatedSeed_Graizer17(Imt imt) {
super(ID, imt);
}
}
+ /**
+ * Stewart et al. site amplification model.
+ *
+ * The model is applicable to 200 ≤ vs30 ≤ 2000 m/s. In the current
+ * implementation, for vs30 < 200, vs30 = 200 m/s; for vs30 > 2000, vs30 =
+ * 3000 m/s.
+ */
+ static final class SiteAmp {
+
+ static final String NAME = NgaEastUsgs_2017.NAME + ": Site Amplification";
+
+ private static final CoefficientContainer COEFFS = new CoefficientContainer(
+ "nga-east-usgs-siteamp.csv");
+
+ private static final double V_REF = 760.0;
+ private static final double V_REF_NL = 3000.0;
+ private static final double VL = 200.0;
+ private static final double VU = 2000.0;
+
+ private final Coefficients c;
+
+ private static final class Coefficients {
+
+ final double c, v1, v2, vf, σvc, σl, σu, f760, f760σ, f3, f4, f5, vc, σc;
+ final Imt imt;
+ final boolean skipShortPeriod;
+
+ Coefficients(Imt imt, CoefficientContainer cc) {
+ this.imt = imt;
+ skipShortPeriod = (imt == Imt.PGA) || (imt.isSA() && imt.period() < 0.06);
+ Map<String, Double> coeffs = cc.get(imt);
+ c = coeffs.get("c");
+ v1 = coeffs.get("V1");
+ v2 = coeffs.get("V2");
+ vf = coeffs.get("Vf");
+ σvc = coeffs.get("sigma_vc");
+ σl = coeffs.get("sigma_l");
+ σu = coeffs.get("sigma_u");
+ f760 = coeffs.get("f760");
+ f760σ = coeffs.get("f760_sigma");
+ f3 = coeffs.get("f3");
+ f4 = coeffs.get("f4");
+ f5 = coeffs.get("f5");
+ vc = coeffs.get("Vc");
+ σc = coeffs.get("sigma_c");
+ }
+ }
+
+ private SiteAmp(Imt imt) {
+ c = new Coefficients(imt, COEFFS);
+ }
+
+ double calc(double pgaRock, double vs30) {
+
+ /*
+ * Developer notes:
+ *
+ * Waiting on short period guidance; currently doing nothing. IMTs 0.075s
+ * and 0.75s are using coefficients originally supplied for 0.08s and
+ * 0.8s. IMTs 0.15s, 0.25s, and 1.5s are using linearly interpolated
+ * coefficients. Site-amp sigma is implementted below, but currently
+ * commented out and unused.
+ *
+ * ---------------------------------------------------------------------
+ *
+ * Comments from Matlab implementation:
+ *
+ * This function is used to compute site amplification in ln units for
+ * inputs Vs30 (in m/s) and PGAr (peak acceleration for 3000 m/s reference
+ * rock in g). Outputs:
+ *
+ * 1) Linear and nonlinear site amplification (FLin, FNonlin)
+ *
+ * 2) Standard deviation (σLin, σNonlin)
+ *
+ * 3) Total site amplification (fT) and total standard deviation (σT)
+ *
+ * Notes:
+ *
+ * 1) We assume fv and f760 are independent, so variance of linear site
+ * amplication is the sum of variance f760 and f760;
+ *
+ * 2) The third value of σc coefficient, which is used to compute
+ * nonlinear site amplication. 0.12 in report, and 0.1 in excel. Here we
+ * use 0.12.
+ *
+ * 3) Standard deviation of fv calculation is updated by adding a new
+ * coefficient vf instead of using v1.
+ */
+
+ /* Short period filtering */
+
+ if (c.skipShortPeriod) {
+ return 0.0;
+ }
+
+ /* Vs30 filtering */
+
+ if (vs30 > VU) {
+ return 0.0;
+ } else if (vs30 < VL) {
+ vs30 = VL;
+ }
+
+ /* Linear response */
+
+ double fv = 0.0;
+ if (vs30 <= c.v1) {
+ fv = c.c * log(c.v1 / V_REF);
+ } else if (vs30 <= c.v2) {
+ fv = c.c * log(vs30 / V_REF);
+ } else {
+ fv = c.c * log(c.v2 / V_REF) + c.c / 2.0 * log(vs30 / c.v2);
+ }
+
+ // double fvσ = 0.0;
+ // if (vs30 < c.vf) {
+ // double σT = c.σl - c.σvc;
+ // double vT = (vs30 - VL) / (c.vf - VL);
+ // fvσ = c.σl - 2.0 * σT * vT + σT * vT * vT;
+ // } else if (vs30 <= c.v2) {
+ // fvσ = c.σvc;
+ // } else {
+ // double vT = (vs30 - c.v2) / (VU - c.v2);
+ // fvσ = c.σvc + (c.σu - c.σvc) * vT * vT;
+ // }
+
+ double fLin = fv + c.f760;
+ // double σLin = sqrt(fvσ * fvσ + c.f760σ * c.f760σ);
+
+ /* Nonlinear response */
+
+ double fNonlin = 0.0;
+ // double σNonlin = 0.0;
+ if (vs30 < c.vc) {
+
+ double f2 = c.f4 *
+ (exp(c.f5 * (min(vs30, V_REF_NL) - 360.0)) -
+ exp(c.f5 * (V_REF_NL - 360.0)));
+ double fT = log((pgaRock + c.f3) / c.f3);
+ fNonlin = f2 * fT;
+
+ // double σf2 = 0.0;
+ // if (vs30 < 300.0) {
+ // σf2 = c.σc;
+ // } else if (vs30 < 1000.0) {
+ // σf2 = c.σc - c.σc / log (1000.0 / 300.0) * log(vs30 /300.0);
+ // }
+ // σNonlin = σf2 * fT;
+ }
+
+ // double σT = sqrt(σLin * σLin + σNonlin * σNonlin);
+ return fLin + fNonlin;
+ }
+ }
}
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/ShahjoueiPezeshk_2016.java b/src/gov/usgs/earthquake/nshmp/gmm/ShahjoueiPezeshk_2016.java
index 862a1d11644c288235edfde70d71bd33f608f3c7..1e234471f10d3f089c3dc939ea875ef801cbfa33 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/ShahjoueiPezeshk_2016.java
+++ b/src/gov/usgs/earthquake/nshmp/gmm/ShahjoueiPezeshk_2016.java
@@ -49,7 +49,7 @@ public final class ShahjoueiPezeshk_2016 implements GroundMotionModel {
* update was provided as the complete (mean and sigma) functional form and is
* therefore implemented outside the NgaEastUsgs_2017 wrapper class.
*/
- static final String NAME = NgaEastUsgs_2017.Seed.NAME + "SP16";
+ static final String NAME = NgaEastUsgs_2017.UpdatedSeed.NAME + "SP16";
static final Constraints CONSTRAINTS = Constraints.builder()
.set(MW, Range.closed(4.0, 8.0))
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-hi.csv b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-hi.csv
index 50770de2a944dcf7dadf9f8de0fc0ffa359450af..78eef9b1aa549f3816642c8fc4981b07096aae99 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-hi.csv
+++ b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-hi.csv
@@ -1,9 +1,9 @@
T,a,b,tau1,tau2,tau3,tau4
0.01,0.6553,0.4446,0.5706,0.5551,0.5214,0.4618
0.02,0.6537,0.4452,0.5706,0.5551,0.5214,0.4618
-0.025,0.6537,0.4452,0.5706,0.5551,0.5214,0.4618
+# 0.025,0.6537,0.4452,0.5706,0.5551,0.5214,0.4618
0.03,0.6521,0.4459,0.5706,0.5551,0.5214,0.4618
-0.04,0.6503,0.4466,0.5706,0.5551,0.5214,0.4618
+# 0.04,0.6503,0.4466,0.5706,0.5551,0.5214,0.4618
0.05,0.6489,0.4473,0.5706,0.5551,0.5214,0.4618
0.075,0.645,0.4489,0.5706,0.5551,0.5214,0.4618
0.1,0.6412,0.4505,0.5706,0.5551,0.5214,0.4618
@@ -20,7 +20,7 @@ T,a,b,tau1,tau2,tau3,tau4
3,0.4565,0.4556,0.5706,0.5551,0.5214,0.4618
4,0.4436,0.4437,0.5706,0.5551,0.5214,0.4618
5,0.4374,0.4381,0.5706,0.5551,0.5214,0.4618
-7.5,0.4325,0.4337,0.5706,0.5551,0.5214,0.4618
-10,0.4317,0.4329,0.5706,0.5551,0.5214,0.4618
+# 7.5,0.4325,0.4337,0.5706,0.5551,0.5214,0.4618
+# 10,0.4317,0.4329,0.5706,0.5551,0.5214,0.4618
PGA,0.6553,0.4446,0.5706,0.5551,0.5214,0.4618
-PGV,0.601,0.4296,0.4919,0.4845,0.4535,0.4333
+# PGV,0.601,0.4296,0.4919,0.4845,0.4535,0.4333
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-lo.csv b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-lo.csv
index 4cc23a7707d061f4a9c41d06b0839698804c77b9..0f30682f07db421bd667ba13077587aa37009f48 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-lo.csv
+++ b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-lo.csv
@@ -1,9 +1,9 @@
T,a,b,tau1,tau2,tau3,tau4
0.01,0.4367,0.2525,0.328,0.2928,0.2439,0.2343
0.02,0.4357,0.2518,0.328,0.2928,0.2439,0.2343
-0.025,0.4357,0.2514,0.328,0.2928,0.2439,0.2343
+# 0.025,0.4357,0.2514,0.328,0.2928,0.2439,0.2343
0.03,0.4347,0.251,0.328,0.2928,0.2439,0.2343
-0.04,0.4334,0.2503,0.328,0.2928,0.2439,0.2343
+# 0.04,0.4334,0.2503,0.328,0.2928,0.2439,0.2343
0.05,0.4326,0.2496,0.328,0.2928,0.2439,0.2343
0.075,0.4301,0.2478,0.328,0.2928,0.2439,0.2343
0.1,0.4277,0.2461,0.328,0.2928,0.2439,0.2343
@@ -20,7 +20,7 @@ T,a,b,tau1,tau2,tau3,tau4
3,0.2958,0.2969,0.328,0.2928,0.2439,0.2343
4,0.2832,0.2831,0.328,0.2928,0.2439,0.2343
5,0.2764,0.2757,0.328,0.2928,0.2439,0.2343
-7.5,0.2703,0.2691,0.328,0.2928,0.2439,0.2343
-10,0.2692,0.2679,0.328,0.2928,0.2439,0.2343
+# 7.5,0.2703,0.2691,0.328,0.2928,0.2439,0.2343
+# 10,0.2692,0.2679,0.328,0.2928,0.2439,0.2343
PGA,0.4367,0.2525,0.328,0.2928,0.2439,0.2343
-PGV,0.4027,0.285,0.2488,0.237,0.2278,0.2081
+# PGV,0.4027,0.285,0.2488,0.237,0.2278,0.2081
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-mid.csv b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-mid.csv
index f265a95bb6728563274ddc223acd81595a1b3a15..17cf1078aecea8834c65c86acb0ba67ed487a9e0 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-mid.csv
+++ b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-mid.csv
@@ -1,9 +1,9 @@
T,a,b,tau1,tau2,tau3,tau4
0.01,0.5423,0.3439,0.4436,0.4169,0.3736,0.3415
0.02,0.541,0.3438,0.4436,0.4169,0.3736,0.3415
-0.025,0.5404,0.3438,0.4436,0.4169,0.3736,0.3415
+# 0.025,0.5404,0.3438,0.4436,0.4169,0.3736,0.3415
0.03,0.5397,0.3437,0.4436,0.4169,0.3736,0.3415
-0.04,0.5382,0.3436,0.4436,0.4169,0.3736,0.3415
+# 0.04,0.5382,0.3436,0.4436,0.4169,0.3736,0.3415
0.05,0.5371,0.3435,0.4436,0.4169,0.3736,0.3415
0.075,0.5339,0.3433,0.4436,0.4169,0.3736,0.3415
0.1,0.5308,0.3431,0.4436,0.4169,0.3736,0.3415
@@ -20,7 +20,7 @@ T,a,b,tau1,tau2,tau3,tau4
3,0.3733,0.3734,0.4436,0.4169,0.3736,0.3415
4,0.3604,0.3604,0.4436,0.4169,0.3736,0.3415
5,0.3538,0.3537,0.4436,0.4169,0.3736,0.3415
-7.5,0.3482,0.3481,0.4436,0.4169,0.3736,0.3415
-10,0.3472,0.3471,0.4436,0.4169,0.3736,0.3415
+# 7.5,0.3482,0.3481,0.4436,0.4169,0.3736,0.3415
+# 10,0.3472,0.3471,0.4436,0.4169,0.3736,0.3415
PGA,0.5423,0.3439,0.4436,0.4169,0.3736,0.3415
-PGV,0.4985,0.3548,0.3633,0.3532,0.334,0.3136
+# PGV,0.4985,0.3548,0.3633,0.3532,0.334,0.3136
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-total.csv b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-total.csv
index 10fa6714ab9ecd1aae45c51bd1ff20e7183f0e4f..f3e9b8ba42ebeb7a63e1f45753403c1db9cee1c9 100644
--- a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-total.csv
+++ b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-sigma-total.csv
@@ -1,9 +1,9 @@
T,tau_M5,phi_M5,tau_M6,phi_M6,tau_M7,phi_M7
0.01,0.432,0.6269,0.3779,0.5168,0.3525,0.5039
0.02,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
-0.025,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
+# 0.025,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
0.03,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
-0.04,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
+# 0.04,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
0.05,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
0.075,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
0.1,0.471,0.6682,0.4385,0.5588,0.4138,0.5462
@@ -20,7 +20,7 @@ T,tau_M5,phi_M5,tau_M6,phi_M6,tau_M7,phi_M7
3,0.4862,0.5813,0.4186,0.6098,0.3952,0.6098
4,0.4904,0.5726,0.4144,0.6003,0.391,0.6003
5,0.4899,0.5651,0.4182,0.5986,0.3949,0.5986
-7.5,0.4803,0.5502,0.4067,0.5982,0.3835,0.5982
-10,0.4666,0.5389,0.3993,0.5885,0.3761,0.5885
+# 7.5,0.4803,0.5502,0.4067,0.5982,0.3835,0.5982
+# 10,0.4666,0.5389,0.3993,0.5885,0.3761,0.5885
PGA,0.432,0.6269,0.3779,0.5168,0.3525,0.5039
-PGV,0.3925,0.5979,0.3612,0.5218,0.3502,0.509
+# PGV,0.3925,0.5979,0.3612,0.5218,0.3502,0.509
diff --git a/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-siteamp.csv b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-siteamp.csv
new file mode 100644
index 0000000000000000000000000000000000000000..bfa5f433e8af89102d642498fe36a7fb790efd72
--- /dev/null
+++ b/src/gov/usgs/earthquake/nshmp/gmm/coeffs/nga-east-usgs-siteamp.csv
@@ -0,0 +1,21 @@
+T,c,V1,V2,Vf,sigma_vc,sigma_l,sigma_u,f760,f760_sigma,f3,f4,f5,Vc,sigma_c
+PGA,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+0.010,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+0.020,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+0.030,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+0.050,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+0.075,0.359,333, 760,333,0.2679,0.2809,0.4724,0.505,0.332,0.16249,-0.50667,-0.00273,2990,0.12
+0.100,0.344,319, 760,319,0.2700,0.2626,0.4702,0.674,0.366,0.15083,-0.44661,-0.00335,2990,0.12
+0.150,0.376,317, 760,317,0.2607,0.2842,0.4020,0.548,0.290,0.13949,-0.37571,-0.00412,2262,0.12
+0.200,0.408,314, 760,314,0.2513,0.3058,0.3337,0.421,0.214,0.12815,-0.30481,-0.00488,1533,0.12
+0.250,0.412,257, 760,282,0.2382,0.2908,0.3574,0.347,0.174,0.12943,-0.26653,-0.00572,1343,0.14
+0.300,0.416,200, 760,250,0.2250,0.2758,0.3810,0.272,0.134,0.13070,-0.22825,-0.00655,1152,0.15
+0.400,0.418,200, 760,250,0.2250,0.2752,0.3807,0.211,0.112,0.09414,-0.11591,-0.00872,1018,0.15
+0.500,0.446,200, 764,280,0.2250,0.3111,0.3228,0.175,0.105,0.09888,-0.07793,-0.01028, 938,0.15
+0.750,0.527,230, 942,280,0.2250,0.3338,0.3080,0.121,0.134,0.07357,-0.01592,-0.01515, 832,0.10
+1.000,0.554,278,1103,300,0.2250,0.3771,0.3609,0.095,0.124,0.04367,-0.00478,-0.01823, 951,0.06
+1.500,0.567,282,1152,300,0.2419,0.4625,0.3746,0.087,0.121,0.02266,-0.00357,-0.01560, 915,0.05
+2.000,0.580,286,1201,300,0.2587,0.5478,0.3883,0.079,0.118,0.00164,-0.00236,-0.01296, 879,0.04
+3.000,0.593,313, 876,313,0.3063,0.5378,0.5510,0.073,0.110,0.00746,-0.00626,-0.01043, 894,0.04
+4.000,0.597,322, 881,322,0.3400,0.4350,0.5846,0.066,0.120,0.00269,-0.00331,-0.01215, 875,0.03
+5.000,0.582,325, 855,325,0.3400,0.4004,0.5871,0.064,0.108,0.00242,-0.00256,-0.01325, 856,0.02