From 740fd4832e7a8087981de58cebe7b8d40358d64c Mon Sep 17 00:00:00 2001
From: Peter Powers <pmpowers@usgs.gov>
Date: Thu, 9 Nov 2023 15:18:54 -0700
Subject: [PATCH] alaska 2007 gmm mods
---
.../nshmp/gmm/AbrahamsonSilva_1997.java | 68 +--
.../nshmp/gmm/BooreEtAl_1997_Alaska.java | 106 +++++
.../nshmp/gmm/CampbellBozorgnia_2003.java | 54 +--
.../gmm/CampbellBozorgnia_2003_Alaska.java | 151 +++++++
.../earthquake/nshmp/gmm/CombinedGmm.java | 16 +-
.../gov/usgs/earthquake/nshmp/gmm/Gmm.java | 61 ++-
.../earthquake/nshmp/gmm/SadighEtAl_1997.java | 61 ++-
.../earthquake/nshmp/model/ModelLoader.java | 4 +-
src/main/resources/gmm/coeffs/BJF97-AK.csv | 8 +
src/main/resources/gmm/coeffs/CB03-AK.csv | 8 +
.../earthquake/nshmp/gmm/LegacySoftRock.java | 4 +-
.../resources/gmm/legacy-vs760-results.csv | 400 +++++++++---------
12 files changed, 630 insertions(+), 311 deletions(-)
create mode 100644 src/main/java/gov/usgs/earthquake/nshmp/gmm/BooreEtAl_1997_Alaska.java
create mode 100644 src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003_Alaska.java
create mode 100644 src/main/resources/gmm/coeffs/BJF97-AK.csv
create mode 100644 src/main/resources/gmm/coeffs/CB03-AK.csv
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonSilva_1997.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonSilva_1997.java
index 0231495f..53910b94 100755
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonSilva_1997.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonSilva_1997.java
@@ -18,6 +18,8 @@ import gov.usgs.earthquake.nshmp.gmm.GmmInput.Constraints;
import gov.usgs.earthquake.nshmp.tree.LogicTree;
/**
+ * Modified to work for Alaska. (11/9/2023 pmpowers)
+ *
* Implementation of the Abrahamson & Silva (1997) ground motion model for
* shallow earthquakes in active continental crust. In keeping with prior NSHMP
* implementations of this older model, only soft rock sites are supported (Vs30
@@ -41,7 +43,7 @@ import gov.usgs.earthquake.nshmp.tree.LogicTree;
*/
public class AbrahamsonSilva_1997 implements GroundMotionModel {
- static final String NAME = "Abrahamson & Silva (1997)";
+ static final String NAME = "Abrahamson & Silva (1997) : Alaska";
static final Constraints CONSTRAINTS = Constraints.builder()
.set(MW, Range.closed(4.0, 8.0))
@@ -103,12 +105,12 @@ public class AbrahamsonSilva_1997 implements GroundMotionModel {
return GroundMotions.createTree(μ, σ);
}
- private static final double calcMean(Coefficients c, GmmInput in) {
+ double calcMean(Coefficients c, GmmInput in) {
// frequently used method locals
double Mw = in.Mw;
double rRup = in.rRup;
- double rake = in.rake;
+ // double rake = in.rake;
// base model
double R = sqrt(rRup * rRup + c.c4 * c.c4);
@@ -120,41 +122,55 @@ public class AbrahamsonSilva_1997 implements GroundMotionModel {
// style-of-faulting term - use rake to determine F; the values below
// were derived from the OpenSHA implementation for which guidance was
// supplied by Norm Abrahamson
- double f3 = 0.0;
- if (rake > 67.5 && rake < 112.5) {
- f3 = calcFaultStyle(c, Mw);
- } else if (rake > 22.5 && rake < 157.5) {
- f3 = 0.5 * calcFaultStyle(c, Mw);
- }
+ // double f3 = 0.0;
+ // if (rake > 67.5 && rake < 112.5) {
+ // f3 = calcFaultStyle(c, Mw);
+ // } else if (rake > 22.5 && rake < 157.5) {
+ // f3 = 0.5 * calcFaultStyle(c, Mw);
+ // }
// hanging-wall term
- double f4 = 0.0;
- if (in.rJB == 0.0 && in.dip < 90.0) {
- double hwM = (Mw >= 6.5) ? 1.0 : (Mw > 5.5) ? Mw - 5.5 : 0.0;
- double hwR = 0.0;
- if (rRup < 4.0) {
- hwR = 0.0;
- } else if (rRup < 8.0) {
- hwR = c.a9 * (rRup - 4.0) / 4.0;
- } else if (rRup < 18.0) {
- hwR = c.a9;
- } else if (rRup < 25.0) {
- hwR = c.a9 * (1.0 - (rRup - 18.0) / 7.0);
- }
- f4 = hwM * hwR;
- }
+ // double f4 = 0.0;
+ // if (in.rJB == 0.0 && in.dip < 90.0) {
+ // double hwM = (Mw >= 6.5) ? 1.0 : (Mw > 5.5) ? Mw - 5.5 : 0.0;
+ // double hwR = 0.0;
+ // if (rRup < 4.0) {
+ // hwR = 0.0;
+ // } else if (rRup < 8.0) {
+ // hwR = c.a9 * (rRup - 4.0) / 4.0;
+ // } else if (rRup < 18.0) {
+ // hwR = c.a9;
+ // } else if (rRup < 25.0) {
+ // hwR = c.a9 * (1.0 - (rRup - 18.0) / 7.0);
+ // }
+ // f4 = hwM * hwR;
+ // }
// no site response term required for rock (f5)
- return f1 + f3 + f4;
+ return f1; // + f3 + f4;
}
private static double calcFaultStyle(Coefficients c, double Mw) {
- return (Mw <= 5.8) ? c.a5 : (Mw < C1) ? c.a5 + (c.a6 - c.a5) / (C1 - 5.8) : c.a6;
+ return (Mw <= 5.8) ? c.a5 : (Mw < C1) ? c.a5 + (c.a6 - c.a5) * (Mw - 5.8) / (C1 - 5.8) : c.a6;
}
private static double calcStdDev(Coefficients c, double Mw) {
return (Mw <= 5.0) ? c.b5 : (Mw < 7.0) ? c.b5 - c.b6 * (Mw - 5.0) : c.b5 - 2 * c.b6;
}
+ static class Gridded extends AbrahamsonSilva_1997 {
+
+ static final String NAME = AbrahamsonSilva_1997.NAME + " : Gridded";
+
+ Gridded(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ double calcMean(Coefficients c, GmmInput in) {
+ return super.calcMean(c, in) + 0.5 * calcFaultStyle(c, in.Mw);
+ }
+ }
+
}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/BooreEtAl_1997_Alaska.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/BooreEtAl_1997_Alaska.java
new file mode 100644
index 00000000..fa57df73
--- /dev/null
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/BooreEtAl_1997_Alaska.java
@@ -0,0 +1,106 @@
+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.RAKE;
+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.log;
+import static java.lang.Math.sqrt;
+
+import java.util.Map;
+
+import com.google.common.collect.Range;
+
+import gov.usgs.earthquake.nshmp.Faults;
+import gov.usgs.earthquake.nshmp.gmm.GmmInput.Constraints;
+import gov.usgs.earthquake.nshmp.tree.LogicTree;
+
+/**
+ * Alaska-specific implementation of the Boore, Joyner & Fumal (1997) ground
+ * motion model for shallow earthquakes in active continental crust. This
+ * implementation was created to match fortran implementation with custom
+ * coefficients for different tectonic settings and styles of faulting. It's not
+ * clear what the relationship between the published 'b1' coefficients is and
+ * the values included in the fortran input files. This implementation is only
+ * for soft rock sites (Vs30 = 760 m/s).
+ *
+ * <p><b>Component:</b> random horizontal
+ *
+ * @author U.S. Geological Survey
+ * @see Gmm#BJF_97_AK
+ */
+public class BooreEtAl_1997_Alaska implements GroundMotionModel {
+
+ static final String NAME = "Boore et al. (1997) : Alaska";
+
+ static final Constraints CONSTRAINTS = Constraints.builder()
+ .set(MW, Range.closed(4.0, 8.0))
+ .set(RJB, Range.closed(0.0, 80.0))
+ .set(RAKE, Faults.RAKE_RANGE)
+ .set(VS30, Range.singleton(760.0))
+ .build();
+
+ static final CoefficientContainer COEFFS = new CoefficientContainer("BJF97-AK.csv");
+
+ private static final class Coefficients {
+
+ final Imt imt;
+ final double b1, b1mt, b2, b3, b5, h, σ;
+
+ Coefficients(Imt imt, CoefficientContainer cc) {
+ this.imt = imt;
+ Map<String, Double> coeffs = cc.get(imt);
+ b1 = coeffs.get("b1");
+ b1mt = coeffs.get("b1mt");
+ b2 = coeffs.get("b2");
+ b3 = coeffs.get("b3");
+ b5 = coeffs.get("b5");
+ h = coeffs.get("h");
+ σ = coeffs.get("sigma");
+ }
+ }
+
+ private final Coefficients coeffs;
+
+ BooreEtAl_1997_Alaska(Imt imt) {
+ coeffs = new Coefficients(imt, COEFFS);
+ }
+
+ boolean interfaceZone() {
+ return false;
+ }
+
+ @Override
+ public Imt imt() {
+ return coeffs.imt;
+ }
+
+ @Override
+ public LogicTree<GroundMotion> calc(GmmInput in) {
+ double μ = calcMean(coeffs, in);
+ return GroundMotions.createTree(μ, coeffs.σ);
+ }
+
+ private final double calcMean(Coefficients c, GmmInput in) {
+ double b1 = interfaceZone() ? c.b1mt : c.b1;
+ double mFac = in.Mw - 6.0;
+ double r = sqrt(in.rJB * in.rJB + c.h * c.h);
+ double μ = b1 + c.b2 * mFac + c.b3 * mFac * mFac + c.b5 * log(r);
+ return μ;
+ }
+
+ static final class InterfaceZone extends BooreEtAl_1997_Alaska {
+
+ static final String NAME = BooreEtAl_1997_Alaska.NAME + " : Interface";
+
+ InterfaceZone(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ boolean interfaceZone() {
+ return true;
+ }
+ }
+
+}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003.java
index 52ac7268..b9eb05a6 100755
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003.java
@@ -136,7 +136,7 @@ public class CampbellBozorgnia_2003 implements GroundMotionModel {
return GroundMotions.createTree(μ, σ);
}
- double calcMean(
+ private final double calcMean(
Coefficients c,
GmmInput in) {
@@ -173,56 +173,4 @@ public class CampbellBozorgnia_2003 implements GroundMotionModel {
return c16 - ((Mw < 7.4) ? 0.07 * Mw : 0.518);
}
- /**
- * Custom version of Campbell & Bozorgnia (2003) designed to match the
- * fortran hazFX version used with the Alaska 2007 NSHM. All crustal faults in
- * the AK NSHM are strike slip so we don't need to worry about style of
- * faulting and hanging wall terms. Also removes dependence on BJF97 for site
- * scaling.
- */
- public class Alaska extends CampbellBozorgnia_2003 {
- static final String NAME = CampbellBozorgnia_2003.NAME + " : Alaska";
-
- Alaska(Imt imt) {
- super(imt);
- }
-
- @Override
- final double calcMean(
- Coefficients c,
- GmmInput in) {
-
- double Mw = in.Mw;
-
- // magnitude scaling
- double mTerm = (8.5 - Mw) * (8.5 - Mw);
- double f1 = c.c2 * Mw + c.c3 * mTerm;
-
- // distance scaling - per verbal communication with Ken Campbell it
- // is ok to approximate rSeis by constraining rRup to be >= 3.0 km.
- // (see also Campbell 1997)
- double rSeis = max(in.rRup, 3.0);
- double rTerm = gS * exp(c.c8 * Mw + c.c9 * mTerm);
- double f2 = rSeis * rSeis + rTerm * rTerm;
-
- // focal mechanism
- // FaultStyle style = GmmUtils.rakeToFaultStyle_NSHMP(in.rake);
- // double f3 = (style == REVERSE) ? revStyleTerm : (style == UNKNOWN) ?
- // unkStyleTerm : 0.0;
-
- // far-source effect of local site conditions
- // (pre-calculated)
-
- // // hanging-wall effect
- // double HW = (in.rJB < 5.0 && in.dip <= 70.0) ? (5.0 - in.rJB) / 5.0 :
- // 0.0;
- // double hwM = (Mw > 6.5) ? 1.0 : (Mw >= 5.5) ? Mw - 5.5 : 0.0;
- // double hwR = (rSeis < 8.0) ? c.c15 * (rSeis / 8.0) : c.c15;
- // double f5 = HW * f3 * hwM * hwR;
-
- // return c.c1 + f1 + c.c4 * log(sqrt(f2)) + f3 + f4 + f5;
- return c.c1 + f1 + c.c4 * log(sqrt(f2)) + c.c14;
- }
-
- }
}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003_Alaska.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003_Alaska.java
new file mode 100644
index 00000000..320671a4
--- /dev/null
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/CampbellBozorgnia_2003_Alaska.java
@@ -0,0 +1,151 @@
+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.RAKE;
+import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.RRUP;
+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.max;
+import static java.lang.Math.sqrt;
+
+import java.util.Map;
+
+import com.google.common.collect.Range;
+
+import gov.usgs.earthquake.nshmp.Faults;
+import gov.usgs.earthquake.nshmp.gmm.GmmInput.Constraints;
+import gov.usgs.earthquake.nshmp.tree.LogicTree;
+
+/**
+ * Alaska-specific implementation of the Boore, Joyner & Fumal (1997) ground
+ * motion model for shallow earthquakes in active continental crust. This
+ * implementation was created to match fortran implementation with custom
+ * coefficients such as 'csite'.
+ *
+ * <p><b>Component:</b> average horizontal
+ *
+ * @author U.S. Geological Survey
+ * @see Gmm#CB_03
+ */
+public class CampbellBozorgnia_2003_Alaska implements GroundMotionModel {
+
+ static final String NAME = "Campbell & Bozorgnia (2003) : Alaska";
+
+ static final Constraints CONSTRAINTS = Constraints.builder()
+ .set(MW, Range.closed(5.8, 8.0))
+ .set(RRUP, Range.closed(0.0, 300.0))
+ .set(RAKE, Faults.RAKE_RANGE)
+ .set(VS30, Range.singleton(760.0))
+ .build();
+
+ static final CoefficientContainer COEFFS;
+
+ static {
+ COEFFS = new CoefficientContainer("CB03-AK.csv");
+ }
+
+ private static final class Coefficients {
+
+ final Imt imt;
+ final double c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11, csite, c15, c16;
+
+ Coefficients(Imt imt, CoefficientContainer cc) {
+ this.imt = imt;
+ Map<String, Double> coeffs = cc.get(imt);
+ c1 = coeffs.get("c1");
+ c2 = coeffs.get("c2");
+ c3 = coeffs.get("c3");
+ c4 = coeffs.get("c4");
+ c5 = coeffs.get("c5");
+ c6 = coeffs.get("c6");
+ c7 = coeffs.get("c7");
+ c8 = coeffs.get("c8");
+ c9 = coeffs.get("c9");
+ c10 = coeffs.get("c10");
+ c11 = coeffs.get("c11");
+ csite = coeffs.get("csite");
+ c15 = coeffs.get("c15");
+ c16 = coeffs.get("c16");
+ }
+ }
+
+ private final Coefficients coeffs;
+ private final double gS;
+
+ CampbellBozorgnia_2003_Alaska(Imt imt) {
+ coeffs = new Coefficients(imt, COEFFS);
+ gS = coeffs.c5 + 0.5 * coeffs.c6 + 0.5 * coeffs.c7;
+ }
+
+ @Override
+ public Imt imt() {
+ return coeffs.imt;
+ }
+
+ @Override
+ public LogicTree<GroundMotion> calc(GmmInput in) {
+ double μ = calcMean(coeffs, in);
+ double σ = calcStdDev(coeffs.c16, in.Mw);
+ return GroundMotions.createTree(μ, σ);
+ }
+
+ double calcMean(
+ Coefficients c,
+ GmmInput in) {
+
+ double Mw = in.Mw;
+
+ // magnitude scaling
+ double mTerm = (8.5 - Mw) * (8.5 - Mw);
+ double f1 = c.c2 * Mw + c.c3 * mTerm;
+
+ // distance scaling - per verbal communication with Ken Campbell it
+ // is ok to approximate rSeis by constraining rRup to be >= 3.0 km.
+ // (see also Campbell 1997)
+ double rSeis = max(in.rRup, 3.0);
+ double rTerm = gS * exp(c.c8 * Mw + c.c9 * mTerm);
+ double f2 = rSeis * rSeis + rTerm * rTerm;
+
+ // // focal mechanism
+ // FaultStyle style = GmmUtils.rakeToFaultStyle_NSHMP(in.rake);
+ // double f3 = (style == REVERSE) ? revStyleTerm : (style == UNKNOWN) ?
+ // unkStyleTerm : 0.0;
+
+ // far-source effect of local site conditions
+ // (pre-calculated)
+
+ // // hanging-wall effect
+ // double HW = (in.rJB < 5.0 && in.dip <= 70.0) ? (5.0 - in.rJB) / 5.0 :
+ // 0.0;
+ // double hwM = (Mw > 6.5) ? 1.0 : (Mw >= 5.5) ? Mw - 5.5 : 0.0;
+ // double hwR = (rSeis < 8.0) ? c.c15 * (rSeis / 8.0) : c.c15;
+ // double f5 = HW * f3 * hwM * hwR;
+
+ return c.c1 + f1 + c.c4 * log(sqrt(f2)) + c.csite; // + f3 + f5;
+ }
+
+ private static final double calcStdDev(double c16, double Mw) {
+ return c16 - ((Mw < 7.4) ? 0.07 * Mw : 0.518);
+ }
+
+ public static class Gridded extends CampbellBozorgnia_2003_Alaska {
+
+ static final String NAME = CampbellBozorgnia_2003_Alaska.NAME + " : Gridded";
+
+ Gridded(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ double calcMean(
+ Coefficients c,
+ GmmInput in) {
+
+ double unkStyleTerm = 0.25 * c.c10 + 0.25 * c.c11;
+ return super.calcMean(c, in) + unkStyleTerm;
+ }
+
+ }
+
+}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/CombinedGmm.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/CombinedGmm.java
index f46045f7..1963e21d 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/CombinedGmm.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/CombinedGmm.java
@@ -17,18 +17,18 @@ import static gov.usgs.earthquake.nshmp.gmm.Gmm.AM_09_INTERFACE_BASIN_M9;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.ASK_14;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.ASK_14_BASIN;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.ASK_14_CYBERSHAKE;
-import static gov.usgs.earthquake.nshmp.gmm.Gmm.AS_97;
+import static gov.usgs.earthquake.nshmp.gmm.Gmm.AS_97_AK;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.ATKINSON_08_PRIME;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.ATKINSON_10;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BCHYDRO_12_INTERFACE_BASIN;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BCHYDRO_12_SLAB;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BCHYDRO_12_SLAB_BASIN;
-import static gov.usgs.earthquake.nshmp.gmm.Gmm.BJF_97;
+import static gov.usgs.earthquake.nshmp.gmm.Gmm.BJF_97_AK;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BSSA_14;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BSSA_14_BASIN;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BSSA_14_CYBERSHAKE;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.CAMPBELL_03;
-import static gov.usgs.earthquake.nshmp.gmm.Gmm.CB_03;
+import static gov.usgs.earthquake.nshmp.gmm.Gmm.CB_03_AK_FAULT;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.CB_14;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.CB_14_BASIN;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.CB_14_CYBERSHAKE;
@@ -65,7 +65,7 @@ import static gov.usgs.earthquake.nshmp.gmm.Gmm.PSBAH_20_CASCADIA_SLAB_BASIN;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.PSBAH_20_GLOBAL_INTERFACE;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.PSBAH_20_GLOBAL_SLAB;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.SADIGH_97;
-import static gov.usgs.earthquake.nshmp.gmm.Gmm.SADIGH_97_INTERFACE_AK;
+import static gov.usgs.earthquake.nshmp.gmm.Gmm.SADIGH_97_AK_INTERFACE;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.SILVA_02;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.SOMERVILLE_01;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.TORO_97_MW;
@@ -494,10 +494,10 @@ class CombinedGmm implements GroundMotionModel {
}
private static final Map<Gmm, Double> AK_ACTIVE_CRUST_2007 = Map.of(
- AS_97, 0.25,
- BJF_97, 0.25,
+ AS_97_AK, 0.25,
+ BJF_97_AK, 0.25,
SADIGH_97, 0.25,
- CB_03, 0.25);
+ CB_03_AK_FAULT, 0.25);
/* 2.0 */
static final class AkActiveCrust2007 extends CombinedGmm {
@@ -525,7 +525,7 @@ class CombinedGmm implements GroundMotionModel {
private static final Map<Gmm, Double> AK_INTERFACE_2007 = Map.of(
YOUNGS_97_INTERFACE, 0.5,
- SADIGH_97_INTERFACE_AK, 0.5);
+ SADIGH_97_AK_INTERFACE, 0.5);
/* 2.0 Sadigh only used out to 70km; only Youngs 70 - 1000 km. */
static final class AkInterface2007 extends CombinedGmm {
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java
index 5f65c4f0..5275126a 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java
@@ -259,12 +259,19 @@ public enum Gmm {
/* Active Continent AK 2007, HI 1998 */
/** @see AbrahamsonSilva_1997 */
- AS_97(
+ AS_97_AK(
AbrahamsonSilva_1997.class,
AbrahamsonSilva_1997.NAME,
AbrahamsonSilva_1997.COEFFS,
AbrahamsonSilva_1997.CONSTRAINTS),
+ /** @see AbrahamsonSilva_1997 */
+ AS_97_AK_GRID(
+ AbrahamsonSilva_1997.Gridded.class,
+ AbrahamsonSilva_1997.Gridded.NAME,
+ AbrahamsonSilva_1997.COEFFS,
+ AbrahamsonSilva_1997.CONSTRAINTS),
+
/** @see BooreEtAl_1997 */
BJF_97(
BooreEtAl_1997.class,
@@ -272,6 +279,20 @@ public enum Gmm {
BooreEtAl_1997.COEFFS,
BooreEtAl_1997.CONSTRAINTS),
+ /** @see BooreEtAl_1997 */
+ BJF_97_AK(
+ BooreEtAl_1997_Alaska.class,
+ BooreEtAl_1997_Alaska.NAME,
+ BooreEtAl_1997_Alaska.COEFFS,
+ BooreEtAl_1997_Alaska.CONSTRAINTS),
+
+ /** @see BooreEtAl_1997 */
+ BJF_97_AK_INTERFACE_ZONE(
+ BooreEtAl_1997_Alaska.InterfaceZone.class,
+ BooreEtAl_1997_Alaska.InterfaceZone.NAME,
+ BooreEtAl_1997_Alaska.COEFFS,
+ BooreEtAl_1997_Alaska.CONSTRAINTS),
+
/** @see Campbell_1997 */
CAMPBELL_97(
Campbell_1997.class,
@@ -287,11 +308,18 @@ public enum Gmm {
CampbellBozorgnia_2003.CONSTRAINTS),
/** @see CampbellBozorgnia_2003 */
- CB_03_ALASKA(
- CampbellBozorgnia_2003.Alaska.class,
- CampbellBozorgnia_2003.Alaska.NAME,
- CampbellBozorgnia_2003.COEFFS,
- CampbellBozorgnia_2003.CONSTRAINTS),
+ CB_03_AK_FAULT(
+ CampbellBozorgnia_2003_Alaska.class,
+ CampbellBozorgnia_2003_Alaska.NAME,
+ CampbellBozorgnia_2003_Alaska.COEFFS,
+ CampbellBozorgnia_2003_Alaska.CONSTRAINTS),
+
+ /** @see CampbellBozorgnia_2003 */
+ CB_03_AK_GRID(
+ CampbellBozorgnia_2003_Alaska.Gridded.class,
+ CampbellBozorgnia_2003_Alaska.Gridded.NAME,
+ CampbellBozorgnia_2003_Alaska.COEFFS,
+ CampbellBozorgnia_2003_Alaska.CONSTRAINTS),
/** @see MunsonThurber_1997 */
MT_97(
@@ -308,9 +336,16 @@ public enum Gmm {
SadighEtAl_1997.CONSTRAINTS),
/** @see SadighEtAl_1997 */
- SADIGH_97_INTERFACE_AK(
- SadighEtAl_1997.InterfaceAlaska.class,
- SadighEtAl_1997.InterfaceAlaska.NAME,
+ SADIGH_97_AK_INTERFACE(
+ SadighEtAl_1997.AlaskaInterface.class,
+ SadighEtAl_1997.AlaskaInterface.NAME,
+ SadighEtAl_1997.COEFFS_BC_HI,
+ SadighEtAl_1997.CONSTRAINTS),
+
+ /** @see SadighEtAl_1997 */
+ SADIGH_97_AK_INTERFACE_ZONE(
+ SadighEtAl_1997.AlaskaInterfaceZone.class,
+ SadighEtAl_1997.AlaskaInterfaceZone.NAME,
SadighEtAl_1997.COEFFS_BC_HI,
SadighEtAl_1997.CONSTRAINTS),
@@ -1863,9 +1898,9 @@ public enum Gmm {
AK_07_ACTIVE_CRUST(
"2007 Active Crust (AK)",
List.of(
- AS_97,
- BJF_97,
- CB_03,
+ AS_97_AK,
+ BJF_97_AK,
+ CB_03_AK_FAULT,
SADIGH_97)),
AK_23_ACTIVE_CRUST(
@@ -1880,7 +1915,7 @@ public enum Gmm {
"2007 Subduction Interface (AK)",
List.of(
YOUNGS_97_INTERFACE,
- SADIGH_97_INTERFACE_AK)),
+ SADIGH_97_AK_INTERFACE)),
AK_07_SLAB(
"2007 Subduction Intraslab (AK)",
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/SadighEtAl_1997.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/SadighEtAl_1997.java
index 6258c601..2d62b5a6 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/SadighEtAl_1997.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/SadighEtAl_1997.java
@@ -140,8 +140,11 @@ public class SadighEtAl_1997 implements GroundMotionModel {
return GroundMotions.createTree(μ, σ);
}
- private static final double calcRockMean(Coefficients c, double Mw,
- double rRup, FaultStyle style) {
+ double calcRockMean(
+ Coefficients c,
+ double Mw,
+ double rRup,
+ FaultStyle style) {
/*
* Rock site coeffs are not dependent on style-of-faulting so we just use
@@ -159,8 +162,11 @@ public class SadighEtAl_1997 implements GroundMotionModel {
return (style == REVERSE) ? lnY + 0.18232 : lnY;
}
- private static final double calcSoilMean(Coefficients c, double Mw,
- double rRup, FaultStyle style) {
+ private static final double calcSoilMean(
+ Coefficients c,
+ double Mw,
+ double rRup,
+ FaultStyle style) {
double c1 = (style == REVERSE) ? c.c1r : c.c1ss;
double c6 = (style == REVERSE) ? c.c6r : c.c6ss;
@@ -187,13 +193,13 @@ public class SadighEtAl_1997 implements GroundMotionModel {
* motions for YOUNGS_97_INTERFACE when rRup > 70 km, obviating the need for a
* distance dependent GMM logic tree.
*/
- static final class InterfaceAlaska extends SadighEtAl_1997 {
+ static final class AlaskaInterface extends SadighEtAl_1997 {
- static final String NAME = SadighEtAl_1997.NAME + " : Interface : Alaska";
+ static final String NAME = SadighEtAl_1997.NAME + " : Alaska : Interface";
private final Geomatrix_1993.Interface delegate;
- InterfaceAlaska(Imt imt) {
+ AlaskaInterface(Imt imt) {
super(imt);
delegate = (Geomatrix_1993.Interface) GEOMATRIX_93_INTERFACE.instance(imt);
}
@@ -205,6 +211,47 @@ public class SadighEtAl_1997 implements GroundMotionModel {
}
return super.calc(in);
}
+
+ @Override
+ double calcRockMean(
+ Coefficients c,
+ double Mw,
+ double rRup,
+ FaultStyle style) {
+
+ /*
+ * We know interface events have a reverse rake. The Sadigh GMM as
+ * published includes a scale factor for reverse events, but it is not
+ * applied in the AK 2007 NSHM on the megathrust so we subtract it here.
+ */
+ return super.calcRockMean(c, Mw, rRup, style) - 0.18232;
+ }
+ }
+
+ /*
+ * Custom version of model for shallow, low magnitude interface zones. These
+ * are modeled as gridded sources with uniform rates. We know that the focal
+ * mechanisms coming in will be 100% strike-slip (and the reverse scale factor
+ * in the published implementation, above, will not be added) so we add half
+ * the reverse scale factor here.
+ */
+ static final class AlaskaInterfaceZone extends SadighEtAl_1997 {
+
+ static final String NAME = SadighEtAl_1997.NAME + " : Alaska : Interface (zone)";
+
+ AlaskaInterfaceZone(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ double calcRockMean(
+ Coefficients c,
+ double Mw,
+ double rRup,
+ FaultStyle style) {
+
+ return super.calcRockMean(c, Mw, rRup, style) + 0.18232 / 2.0;
+ }
}
}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/model/ModelLoader.java b/src/main/java/gov/usgs/earthquake/nshmp/model/ModelLoader.java
index f506a61d..e2be9870 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/model/ModelLoader.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/model/ModelLoader.java
@@ -104,9 +104,9 @@ abstract class ModelLoader {
public static void main(String[] args) throws IOException {
// Path testModel = Paths.get("../nshm-conus-2018-5.2.x-maint");
- Path testModel = Paths.get("../nshm-conus");
+ // Path testModel = Paths.get("../nshm-conus");
// Path testModel = Paths.get("../nshm-alaska");
- // Path testModel = Paths.get("../nshm-alaska-2007-2.2.x-maint");
+ Path testModel = Paths.get("../nshm-alaska-2007-2.2.x-maint");
// Path testModel = Paths.get("../nshm-hawaii");
// Allow ModelLoader to be run from the command line for checking models
diff --git a/src/main/resources/gmm/coeffs/BJF97-AK.csv b/src/main/resources/gmm/coeffs/BJF97-AK.csv
new file mode 100644
index 00000000..f287f6aa
--- /dev/null
+++ b/src/main/resources/gmm/coeffs/BJF97-AK.csv
@@ -0,0 +1,8 @@
+T, b1, b1mt, b2, b3, b5, h, sigma
+PGA, -0.08741, 0.01059, 0.527, 0.000, -0.778, 5.57, 0.520
+0.1, 1.08669, 1.12719, 0.753, -0.226, -0.934, 6.27, 0.479
+0.2, 1.29827, 1.38377, 0.711, -0.207, -0.924, 7.02, 0.502
+0.3, 1.01182, 1.11432, 0.769, -0.161, -0.893, 5.94, 0.522
+0.5, 0.34925, 0.45375, 0.884, -0.090, -0.846, 4.13, 0.556
+1, -0.70360, -0.64160, 1.036, -0.032, -0.798, 2.90, 0.613
+2, -1.13607, -1.18707, 1.085, -0.085, -0.812, 5.85, 0.672
diff --git a/src/main/resources/gmm/coeffs/CB03-AK.csv b/src/main/resources/gmm/coeffs/CB03-AK.csv
new file mode 100644
index 00000000..27fe17d6
--- /dev/null
+++ b/src/main/resources/gmm/coeffs/CB03-AK.csv
@@ -0,0 +1,8 @@
+T, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11, csite, c15, c16
+PGA, -4.033, 0.812, 0.036, -1.061, 0.041, -0.005, -0.018, 0.766, 0.034, 0.343, 0.351, -0.289, 0.370, 0.920
+0.1, -2.661, 0.812, 0.060, -1.308, 0.166, -0.009, -0.068, 0.621, 0.046, 0.224, 0.313, -0.299, 0.370, 0.958
+0.2, -2.771, 0.812, 0.030, -1.153, 0.098, -0.014, -0.038, 0.704, 0.026, 0.296, 0.342, -0.330, 0.370, 0.981
+0.3, -2.999, 0.812, 0.007, -1.080, 0.059, -0.007, -0.022, 0.752, 0.007, 0.359, 0.385, -0.453, 0.370, 0.984
+0.5, -3.556, 0.812, -0.035, -0.964, 0.023, -0.002, -0.004, 0.842, -0.036, 0.406, 0.479, -0.528, 0.370, 0.990
+1, -3.867, 0.812, -0.101, -0.964, 0.019, 0.000, 0.000, 0.842, -0.105, 0.329, 0.338, -0.607, 0.281, 1.021
+2, -4.311, 0.812, -0.180, -0.964, 0.019, 0.000, 0.000, 0.842, -0.187, 0.060, 0.064, -0.649, 0.160, 1.021
diff --git a/src/test/java/gov/usgs/earthquake/nshmp/gmm/LegacySoftRock.java b/src/test/java/gov/usgs/earthquake/nshmp/gmm/LegacySoftRock.java
index 9d80aee3..8a27fa04 100644
--- a/src/test/java/gov/usgs/earthquake/nshmp/gmm/LegacySoftRock.java
+++ b/src/test/java/gov/usgs/earthquake/nshmp/gmm/LegacySoftRock.java
@@ -1,6 +1,6 @@
package gov.usgs.earthquake.nshmp.gmm;
-import static gov.usgs.earthquake.nshmp.gmm.Gmm.AS_97;
+import static gov.usgs.earthquake.nshmp.gmm.Gmm.AS_97_AK;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.BJF_97;
import static gov.usgs.earthquake.nshmp.gmm.Gmm.CB_03;
import static gov.usgs.earthquake.nshmp.gmm.Imt.PGA;
@@ -42,7 +42,7 @@ class LegacySoftRock implements ArgumentsProvider {
}
/* Result generation sets */
- private static Set<Gmm> gmms = EnumSet.of(AS_97, BJF_97, CB_03);
+ private static Set<Gmm> gmms = EnumSet.of(AS_97_AK, BJF_97, CB_03);
private static Set<Imt> imts = EnumSet.of(PGA, SA0P2, SA1P0, SA2P0);
public static void main(String[] args) throws IOException {
diff --git a/src/test/resources/gmm/legacy-vs760-results.csv b/src/test/resources/gmm/legacy-vs760-results.csv
index 7dfed57c..ea180d0b 100644
--- a/src/test/resources/gmm/legacy-vs760-results.csv
+++ b/src/test/resources/gmm/legacy-vs760-results.csv
@@ -1,203 +1,203 @@
-0-AS_97-PGA,0.0912346815,0.7000000000
-1-AS_97-PGA,0.0664055757,0.7000000000
-2-AS_97-PGA,0.1388495079,0.7000000000
-3-AS_97-PGA,0.1027419465,0.7000000000
-4-AS_97-PGA,0.0418660886,0.7000000000
-5-AS_97-PGA,0.2113142229,0.7000000000
-6-AS_97-PGA,0.1589611632,0.7000000000
-7-AS_97-PGA,0.0680369976,0.7000000000
-8-AS_97-PGA,0.3592072340,0.7000000000
-9-AS_97-PGA,0.2615721705,0.7000000000
-10-AS_97-PGA,0.1165211652,0.7000000000
-11-AS_97-PGA,0.0464803589,0.7000000000
-12-AS_97-PGA,0.5429719391,0.6325000000
-13-AS_97-PGA,0.4031719061,0.6325000000
-14-AS_97-PGA,0.1887503165,0.6325000000
-15-AS_97-PGA,0.0796679124,0.6325000000
-16-AS_97-PGA,0.4962584824,0.5650000000
-17-AS_97-PGA,0.3917444428,0.5650000000
-18-AS_97-PGA,0.1863551774,0.5650000000
-19-AS_97-PGA,0.0855729261,0.5650000000
-20-AS_97-PGA,0.0175349559,0.5650000000
-21-AS_97-PGA,0.9079301614,0.4975000000
-22-AS_97-PGA,0.7656635286,0.4975000000
-23-AS_97-PGA,0.3651683908,0.4975000000
-24-AS_97-PGA,0.1770816384,0.4975000000
-25-AS_97-PGA,0.0405492780,0.4975000000
-26-AS_97-PGA,1.1890855989,0.4300000000
-27-AS_97-PGA,0.9756178802,0.4300000000
-28-AS_97-PGA,0.5183018371,0.4300000000
-29-AS_97-PGA,0.2734335445,0.4300000000
-30-AS_97-PGA,0.0664554868,0.4300000000
-31-AS_97-PGA,0.0420638415,0.4300000000
-32-AS_97-PGA,0.0078553603,0.4300000000
-33-AS_97-PGA,1.5037915207,0.4300000000
-34-AS_97-PGA,1.2741791096,0.4300000000
-35-AS_97-PGA,0.7201792097,0.4300000000
-36-AS_97-PGA,0.4050046633,0.4300000000
-37-AS_97-PGA,0.1082819875,0.4300000000
-38-AS_97-PGA,0.0703941728,0.4300000000
-39-AS_97-PGA,0.0346926389,0.4300000000
-40-AS_97-PGA,0.0147324541,0.4300000000
-41-AS_97-PGA,1.8709096570,0.4300000000
-42-AS_97-PGA,1.6087189314,0.4300000000
-43-AS_97-PGA,0.9564798077,0.4300000000
-44-AS_97-PGA,0.5660765353,0.4300000000
-45-AS_97-PGA,0.1701386980,0.4300000000
-46-AS_97-PGA,0.1149151228,0.4300000000
-47-AS_97-PGA,0.0603036242,0.4300000000
-48-AS_97-PGA,0.0276302038,0.4300000000
-49-AS_97-PGA,0.4962584824,0.5650000000
-0-AS_97-SA0P2,0.1621752434,0.7700000000
-1-AS_97-SA0P2,0.1149595051,0.7700000000
-2-AS_97-SA0P2,0.2621950197,0.7700000000
-3-AS_97-SA0P2,0.1892710756,0.7700000000
-4-AS_97-SA0P2,0.0752372131,0.7700000000
-5-AS_97-SA0P2,0.4209860313,0.7700000000
-6-AS_97-SA0P2,0.3094760058,0.7700000000
-7-AS_97-SA0P2,0.1295197974,0.7700000000
-8-AS_97-SA0P2,0.7593927170,0.7700000000
-9-AS_97-SA0P2,0.5365762515,0.7700000000
-10-AS_97-SA0P2,0.2333987779,0.7700000000
-11-AS_97-SA0P2,0.0937358684,0.7700000000
-12-AS_97-SA0P2,1.1933067226,0.7025000000
-13-AS_97-SA0P2,0.8617733636,0.7025000000
-14-AS_97-SA0P2,0.3949779986,0.7025000000
-15-AS_97-SA0P2,0.1679849187,0.7025000000
-16-AS_97-SA0P2,1.1383857858,0.6350000000
-17-AS_97-SA0P2,0.8711647548,0.6350000000
-18-AS_97-SA0P2,0.4047163141,0.6350000000
-19-AS_97-SA0P2,0.1869851122,0.6350000000
-20-AS_97-SA0P2,0.0396649729,0.6350000000
-21-AS_97-SA0P2,2.1334911339,0.5675000000
-22-AS_97-SA0P2,1.7489578410,0.5675000000
-23-AS_97-SA0P2,0.8170022623,0.5675000000
-24-AS_97-SA0P2,0.3990014112,0.5675000000
-25-AS_97-SA0P2,0.0946148206,0.5675000000
-26-AS_97-SA0P2,2.8628598407,0.5000000000
-27-AS_97-SA0P2,2.2869133758,0.5000000000
-28-AS_97-SA0P2,1.1866345100,0.5000000000
-29-AS_97-SA0P2,0.6298389975,0.5000000000
-30-AS_97-SA0P2,0.1586868604,0.5000000000
-31-AS_97-SA0P2,0.1017531523,0.5000000000
-32-AS_97-SA0P2,0.0199392334,0.5000000000
-33-AS_97-SA0P2,3.6644006155,0.5000000000
-34-AS_97-SA0P2,3.0375900570,0.5000000000
-35-AS_97-SA0P2,1.6753527250,0.5000000000
-36-AS_97-SA0P2,0.9476721579,0.5000000000
-37-AS_97-SA0P2,0.2627855730,0.5000000000
-38-AS_97-SA0P2,0.1731229333,0.5000000000
-39-AS_97-SA0P2,0.0872235836,0.5000000000
-40-AS_97-SA0P2,0.0380459844,0.5000000000
-41-AS_97-SA0P2,4.5699406822,0.5000000000
-42-AS_97-SA0P2,3.8538903618,0.5000000000
-43-AS_97-SA0P2,2.2446397564,0.5000000000
-44-AS_97-SA0P2,1.3377146556,0.5000000000
-45-AS_97-SA0P2,0.4171866340,0.5000000000
-46-AS_97-SA0P2,0.2855516994,0.5000000000
-47-AS_97-SA0P2,0.1531911750,0.5000000000
-48-AS_97-SA0P2,0.0720962311,0.5000000000
-49-AS_97-SA0P2,1.1383857858,0.6350000000
-0-AS_97-SA1P0,0.0069161718,0.8300000000
-1-AS_97-SA1P0,0.0047547677,0.8300000000
-2-AS_97-SA1P0,0.0167512327,0.8300000000
-3-AS_97-SA1P0,0.0117950721,0.8300000000
-4-AS_97-SA1P0,0.0050383156,0.8300000000
-5-AS_97-SA1P0,0.0385548262,0.8300000000
-6-AS_97-SA1P0,0.0278049983,0.8300000000
-7-AS_97-SA1P0,0.0125864269,0.8300000000
-8-AS_97-SA1P0,0.0974550960,0.8300000000
-9-AS_97-SA1P0,0.0664845917,0.8300000000
-10-AS_97-SA1P0,0.0312389105,0.8300000000
-11-AS_97-SA1P0,0.0147384202,0.8300000000
-12-AS_97-SA1P0,0.2002513846,0.7710000000
-13-AS_97-SA1P0,0.1408018321,0.7710000000
-14-AS_97-SA1P0,0.0702245107,0.7710000000
-15-AS_97-SA1P0,0.0351567509,0.7710000000
-16-AS_97-SA1P0,0.1976614170,0.7120000000
-17-AS_97-SA1P0,0.1436031615,0.7120000000
-18-AS_97-SA1P0,0.0728466928,0.7120000000
-19-AS_97-SA1P0,0.0393333804,0.7120000000
-20-AS_97-SA1P0,0.0118931199,0.7120000000
-21-AS_97-SA1P0,0.4807551183,0.6530000000
-22-AS_97-SA1P0,0.3739944180,0.6530000000
-23-AS_97-SA1P0,0.1945673889,0.6530000000
-24-AS_97-SA1P0,0.1113073102,0.6530000000
-25-AS_97-SA1P0,0.0376511118,0.6530000000
-26-AS_97-SA1P0,0.7472312819,0.5940000000
-27-AS_97-SA1P0,0.5734005920,0.5940000000
-28-AS_97-SA1P0,0.3246600006,0.5940000000
-29-AS_97-SA1P0,0.2003710889,0.5940000000
-30-AS_97-SA1P0,0.0729804845,0.5940000000
-31-AS_97-SA1P0,0.0528187854,0.5940000000
-32-AS_97-SA1P0,0.0161509580,0.5940000000
-33-AS_97-SA1P0,1.0463237546,0.5940000000
-34-AS_97-SA1P0,0.8422820986,0.5940000000
-35-AS_97-SA1P0,0.5050381646,0.5940000000
-36-AS_97-SA1P0,0.3314040352,0.5940000000
-37-AS_97-SA1P0,0.1336432600,0.5940000000
-38-AS_97-SA1P0,0.0996795945,0.5940000000
-39-AS_97-SA1P0,0.0616005380,0.5940000000
-40-AS_97-SA1P0,0.0344095632,0.5940000000
-41-AS_97-SA1P0,1.3566129620,0.5940000000
-42-AS_97-SA1P0,1.1234207279,0.5940000000
-43-AS_97-SA1P0,0.7201107429,0.5940000000
-44-AS_97-SA1P0,0.4992424467,0.5940000000
-45-AS_97-SA1P0,0.2266598578,0.5940000000
-46-AS_97-SA1P0,0.1756518927,0.5940000000
-47-AS_97-SA1P0,0.1155873036,0.5940000000
-48-AS_97-SA1P0,0.0696644178,0.5940000000
-49-AS_97-SA1P0,0.1976614170,0.7120000000
-0-AS_97-SA2P0,0.0011482389,0.8500000000
-1-AS_97-SA2P0,0.0008042754,0.8500000000
-2-AS_97-SA2P0,0.0033241893,0.8500000000
-3-AS_97-SA2P0,0.0023869825,0.8500000000
-4-AS_97-SA2P0,0.0010885972,0.8500000000
-5-AS_97-SA2P0,0.0089730201,0.8500000000
-6-AS_97-SA2P0,0.0066053073,0.8500000000
-7-AS_97-SA2P0,0.0031951568,0.8500000000
-8-AS_97-SA2P0,0.0259019805,0.8500000000
-9-AS_97-SA2P0,0.0180986501,0.8500000000
-10-AS_97-SA2P0,0.0091019753,0.8500000000
-11-AS_97-SA2P0,0.0046347095,0.8500000000
-12-AS_97-SA2P0,0.0600522251,0.7975000000
-13-AS_97-SA2P0,0.0433095968,0.7975000000
-14-AS_97-SA2P0,0.0231431191,0.7975000000
-15-AS_97-SA2P0,0.0125077914,0.7975000000
-16-AS_97-SA2P0,0.0633796674,0.7450000000
-17-AS_97-SA2P0,0.0465358484,0.7450000000
-18-AS_97-SA2P0,0.0256784007,0.7450000000
-19-AS_97-SA2P0,0.0149355537,0.7450000000
-20-AS_97-SA2P0,0.0052399435,0.7450000000
-21-AS_97-SA2P0,0.1692150497,0.6925000000
-22-AS_97-SA2P0,0.1313028504,0.6925000000
-23-AS_97-SA2P0,0.0755293244,0.6925000000
-24-AS_97-SA2P0,0.0465581018,0.6925000000
-25-AS_97-SA2P0,0.0182750996,0.6925000000
-26-AS_97-SA2P0,0.2787290161,0.6400000000
-27-AS_97-SA2P0,0.2180722785,0.6400000000
-28-AS_97-SA2P0,0.1337493055,0.6400000000
-29-AS_97-SA2P0,0.0887026440,0.6400000000
-30-AS_97-SA2P0,0.0377098861,0.6400000000
-31-AS_97-SA2P0,0.0286833674,0.6400000000
-32-AS_97-SA2P0,0.0105248113,0.6400000000
-33-AS_97-SA2P0,0.4067378204,0.6400000000
-34-AS_97-SA2P0,0.3349654927,0.6400000000
-35-AS_97-SA2P0,0.2169398716,0.6400000000
-36-AS_97-SA2P0,0.1528418888,0.6400000000
-37-AS_97-SA2P0,0.0721281691,0.6400000000
-38-AS_97-SA2P0,0.0566113023,0.6400000000
-39-AS_97-SA2P0,0.0380395895,0.6400000000
-40-AS_97-SA2P0,0.0235138706,0.6400000000
-41-AS_97-SA2P0,0.5400456399,0.6400000000
-42-AS_97-SA2P0,0.4586031555,0.6400000000
-43-AS_97-SA2P0,0.3181144043,0.6400000000
-44-AS_97-SA2P0,0.2368735366,0.6400000000
-45-AS_97-SA2P0,0.1258653984,0.6400000000
-46-AS_97-SA2P0,0.1026421377,0.6400000000
-47-AS_97-SA2P0,0.0734409892,0.6400000000
-48-AS_97-SA2P0,0.0489816392,0.6400000000
-49-AS_97-SA2P0,0.0633796674,0.7450000000
+0-AS_97_AK-PGA,0.0495724435,0.7000000000
+1-AS_97_AK-PGA,0.0360815273,0.7000000000
+2-AS_97_AK-PGA,0.0754440008,0.7000000000
+3-AS_97_AK-PGA,0.0558249259,0.7000000000
+4-AS_97_AK-PGA,0.0227479756,0.7000000000
+5-AS_97_AK-PGA,0.1148177667,0.7000000000
+6-AS_97_AK-PGA,0.0863716862,0.7000000000
+7-AS_97_AK-PGA,0.0369679618,0.7000000000
+8-AS_97_AK-PGA,0.1951755627,0.7000000000
+9-AS_97_AK-PGA,0.1421254662,0.7000000000
+10-AS_97_AK-PGA,0.0633118764,0.7000000000
+11-AS_97_AK-PGA,0.0252551434,0.7000000000
+12-AS_97_AK-PGA,0.2950242750,0.6325000000
+13-AS_97_AK-PGA,0.2190638056,0.6325000000
+14-AS_97_AK-PGA,0.1025576485,0.6325000000
+15-AS_97_AK-PGA,0.0432876294,0.6325000000
+16-AS_97_AK-PGA,0.4831998124,0.5650000000
+17-AS_97_AK-PGA,0.3641963076,0.5650000000
+18-AS_97_AK-PGA,0.1814513806,0.5650000000
+19-AS_97_AK-PGA,0.0833211387,0.5650000000
+20-AS_97_AK-PGA,0.0170735367,0.5650000000
+21-AS_97_AK-PGA,0.7000609907,0.4975000000
+22-AS_97_AK-PGA,0.5382067691,0.4975000000
+23-AS_97_AK-PGA,0.2815636668,0.4975000000
+24-AS_97_AK-PGA,0.1365390781,0.4975000000
+25-AS_97_AK-PGA,0.0312655851,0.4975000000
+26-AS_97_AK-PGA,0.9168463367,0.4300000000
+27-AS_97_AK-PGA,0.7384629958,0.4300000000
+28-AS_97_AK-PGA,0.3996374534,0.4300000000
+29-AS_97_AK-PGA,0.2108313681,0.4300000000
+30-AS_97_AK-PGA,0.0512406085,0.4300000000
+31-AS_97_AK-PGA,0.0324333917,0.4300000000
+32-AS_97_AK-PGA,0.0060568880,0.4300000000
+33-AS_97_AK-PGA,1.1595008368,0.4300000000
+34-AS_97_AK-PGA,0.9824578231,0.4300000000
+35-AS_97_AK-PGA,0.5552953217,0.4300000000
+36-AS_97_AK-PGA,0.3122794879,0.4300000000
+37-AS_97_AK-PGA,0.0834909982,0.4300000000
+38-AS_97_AK-PGA,0.0542775386,0.4300000000
+39-AS_97_AK-PGA,0.0267498143,0.4300000000
+40-AS_97_AK-PGA,0.0113594821,0.4300000000
+41-AS_97_AK-PGA,1.4425678579,0.4300000000
+42-AS_97_AK-PGA,1.2404052832,0.4300000000
+43-AS_97_AK-PGA,0.7374952725,0.4300000000
+44-AS_97_AK-PGA,0.4364742102,0.4300000000
+45-AS_97_AK-PGA,0.1311857129,0.4300000000
+46-AS_97_AK-PGA,0.0886054877,0.4300000000
+47-AS_97_AK-PGA,0.0464972051,0.4300000000
+48-AS_97_AK-PGA,0.0213043124,0.4300000000
+49-AS_97_AK-PGA,0.4831998124,0.5650000000
+0-AS_97_AK-SA0P2,0.0881180595,0.7700000000
+1-AS_97_AK-SA0P2,0.0624633470,0.7700000000
+2-AS_97_AK-SA0P2,0.1424638918,0.7700000000
+3-AS_97_AK-SA0P2,0.1028406034,0.7700000000
+4-AS_97_AK-SA0P2,0.0408802051,0.7700000000
+5-AS_97_AK-SA0P2,0.2287431260,0.7700000000
+6-AS_97_AK-SA0P2,0.1681540567,0.7700000000
+7-AS_97_AK-SA0P2,0.0703746945,0.7700000000
+8-AS_97_AK-SA0P2,0.4126166927,0.7700000000
+9-AS_97_AK-SA0P2,0.2915491726,0.7700000000
+10-AS_97_AK-SA0P2,0.1268174288,0.7700000000
+11-AS_97_AK-SA0P2,0.0509314655,0.7700000000
+12-AS_97_AK-SA0P2,0.6483842448,0.7025000000
+13-AS_97_AK-SA0P2,0.4682453060,0.7025000000
+14-AS_97_AK-SA0P2,0.2146116388,0.7025000000
+15-AS_97_AK-SA0P2,0.0912747515,0.7025000000
+16-AS_97_AK-SA0P2,1.1084300169,0.6350000000
+17-AS_97_AK-SA0P2,0.8099029682,0.6350000000
+18-AS_97_AK-SA0P2,0.3940665076,0.6350000000
+19-AS_97_AK-SA0P2,0.1820647391,0.6350000000
+20-AS_97_AK-SA0P2,0.0386212188,0.6350000000
+21-AS_97_AK-SA0P2,1.6450317221,0.5675000000
+22-AS_97_AK-SA0P2,1.2293924338,0.5675000000
+23-AS_97_AK-SA0P2,0.6299508900,0.5675000000
+24-AS_97_AK-SA0P2,0.3076506708,0.5675000000
+25-AS_97_AK-SA0P2,0.0729529074,0.5675000000
+26-AS_97_AK-SA0P2,2.2074126201,0.5000000000
+27-AS_97_AK-SA0P2,1.7310065109,0.5000000000
+28-AS_97_AK-SA0P2,0.9149564207,0.5000000000
+29-AS_97_AK-SA0P2,0.4856383579,0.5000000000
+30-AS_97_AK-SA0P2,0.1223557554,0.5000000000
+31-AS_97_AK-SA0P2,0.0784569294,0.5000000000
+32-AS_97_AK-SA0P2,0.0153741775,0.5000000000
+33-AS_97_AK-SA0P2,2.8254419056,0.5000000000
+34-AS_97_AK-SA0P2,2.3421386305,0.5000000000
+35-AS_97_AK-SA0P2,1.2917833754,0.5000000000
+36-AS_97_AK-SA0P2,0.7307041202,0.5000000000
+37-AS_97_AK-SA0P2,0.2026212328,0.5000000000
+38-AS_97_AK-SA0P2,0.1334867123,0.5000000000
+39-AS_97_AK-SA0P2,0.0672538825,0.5000000000
+40-AS_97_AK-SA0P2,0.0293354166,0.5000000000
+41-AS_97_AK-SA0P2,3.5236600101,0.5000000000
+42-AS_97_AK-SA0P2,2.9715482750,0.5000000000
+43-AS_97_AK-SA0P2,1.7307330438,0.5000000000
+44-AS_97_AK-SA0P2,1.0314470065,0.5000000000
+45-AS_97_AK-SA0P2,0.3216724157,0.5000000000
+46-AS_97_AK-SA0P2,0.2201750907,0.5000000000
+47-AS_97_AK-SA0P2,0.1181182984,0.5000000000
+48-AS_97_AK-SA0P2,0.0555899133,0.5000000000
+49-AS_97_AK-SA0P2,1.1084300169,0.6350000000
+0-AS_97_AK-SA1P0,0.0042370294,0.8300000000
+1-AS_97_AK-SA1P0,0.0029128962,0.8300000000
+2-AS_97_AK-SA1P0,0.0102622473,0.8300000000
+3-AS_97_AK-SA1P0,0.0072259725,0.8300000000
+4-AS_97_AK-SA1P0,0.0030866051,0.8300000000
+5-AS_97_AK-SA1P0,0.0236197041,0.8300000000
+6-AS_97_AK-SA1P0,0.0170340759,0.8300000000
+7-AS_97_AK-SA1P0,0.0077107773,0.8300000000
+8-AS_97_AK-SA1P0,0.0597035640,0.8300000000
+9-AS_97_AK-SA1P0,0.0407302157,0.8300000000
+10-AS_97_AK-SA1P0,0.0191377811,0.8300000000
+11-AS_97_AK-SA1P0,0.0090291452,0.8300000000
+12-AS_97_AK-SA1P0,0.1226792837,0.7710000000
+13-AS_97_AK-SA1P0,0.0862589187,0.7710000000
+14-AS_97_AK-SA1P0,0.0430213888,0.7710000000
+15-AS_97_AK-SA1P0,0.0215379535,0.7710000000
+16-AS_97_AK-SA1P0,0.2681524204,0.7120000000
+17-AS_97_AK-SA1P0,0.1880915237,0.7120000000
+18-AS_97_AK-SA1P0,0.0988256448,0.7120000000
+19-AS_97_AK-SA1P0,0.0533606473,0.7120000000
+20-AS_97_AK-SA1P0,0.0161345039,0.7120000000
+21-AS_97_AK-SA1P0,0.4745457501,0.6530000000
+22-AS_97_AK-SA1P0,0.3441201508,0.6530000000
+23-AS_97_AK-SA1P0,0.1920543828,0.6530000000
+24-AS_97_AK-SA1P0,0.1098696800,0.6530000000
+25-AS_97_AK-SA1P0,0.0371648151,0.6530000000
+26-AS_97_AK-SA1P0,0.7375801436,0.5940000000
+27-AS_97_AK-SA1P0,0.5580980066,0.5940000000
+28-AS_97_AK-SA1P0,0.3204667358,0.5940000000
+29-AS_97_AK-SA1P0,0.1977831229,0.5940000000
+30-AS_97_AK-SA1P0,0.0720378784,0.5940000000
+31-AS_97_AK-SA1P0,0.0521365851,0.5940000000
+32-AS_97_AK-SA1P0,0.0159423544,0.5940000000
+33-AS_97_AK-SA1P0,1.0328095782,0.5940000000
+34-AS_97_AK-SA1P0,0.8314032967,0.5940000000
+35-AS_97_AK-SA1P0,0.4985151599,0.5940000000
+36-AS_97_AK-SA1P0,0.3271236654,0.5940000000
+37-AS_97_AK-SA1P0,0.1319171417,0.5940000000
+38-AS_97_AK-SA1P0,0.0983921463,0.5940000000
+39-AS_97_AK-SA1P0,0.0608049138,0.5940000000
+40-AS_97_AK-SA1P0,0.0339651339,0.5940000000
+41-AS_97_AK-SA1P0,1.3390911321,0.5940000000
+42-AS_97_AK-SA1P0,1.1089107775,0.5940000000
+43-AS_97_AK-SA1P0,0.7108098898,0.5940000000
+44-AS_97_AK-SA1P0,0.4927942987,0.5940000000
+45-AS_97_AK-SA1P0,0.2237323497,0.5940000000
+46-AS_97_AK-SA1P0,0.1733831965,0.5940000000
+47-AS_97_AK-SA1P0,0.1140943936,0.5940000000
+48-AS_97_AK-SA1P0,0.0687646416,0.5940000000
+49-AS_97_AK-SA1P0,0.2681524204,0.7120000000
+0-AS_97_AK-SA2P0,0.0007696875,0.8500000000
+1-AS_97_AK-SA2P0,0.0005391220,0.8500000000
+2-AS_97_AK-SA2P0,0.0022282707,0.8500000000
+3-AS_97_AK-SA2P0,0.0016000422,0.8500000000
+4-AS_97_AK-SA2P0,0.0007297085,0.8500000000
+5-AS_97_AK-SA2P0,0.0060147952,0.8500000000
+6-AS_97_AK-SA2P0,0.0044276699,0.8500000000
+7-AS_97_AK-SA2P0,0.0021417776,0.8500000000
+8-AS_97_AK-SA2P0,0.0173626168,0.8500000000
+9-AS_97_AK-SA2P0,0.0121318880,0.8500000000
+10-AS_97_AK-SA2P0,0.0061012365,0.8500000000
+11-AS_97_AK-SA2P0,0.0031067387,0.8500000000
+12-AS_97_AK-SA2P0,0.0402542103,0.7975000000
+13-AS_97_AK-SA2P0,0.0290312909,0.7975000000
+14-AS_97_AK-SA2P0,0.0155132966,0.7975000000
+15-AS_97_AK-SA2P0,0.0083842233,0.7975000000
+16-AS_97_AK-SA2P0,0.0967834918,0.7450000000
+17-AS_97_AK-SA2P0,0.0696551252,0.7450000000
+18-AS_97_AK-SA2P0,0.0392120278,0.7450000000
+19-AS_97_AK-SA2P0,0.0228072360,0.7450000000
+20-AS_97_AK-SA2P0,0.0080016202,0.7450000000
+21-AS_97_AK-SA2P0,0.1858928420,0.6925000000
+22-AS_97_AK-SA2P0,0.1385881369,0.6925000000
+23-AS_97_AK-SA2P0,0.0829734754,0.6925000000
+24-AS_97_AK-SA2P0,0.0511468565,0.6925000000
+25-AS_97_AK-SA2P0,0.0200762887,0.6925000000
+26-AS_97_AK-SA2P0,0.3062004771,0.6400000000
+27-AS_97_AK-SA2P0,0.2376565491,0.6400000000
+28-AS_97_AK-SA2P0,0.1469316030,0.6400000000
+29-AS_97_AK-SA2P0,0.0974451540,0.6400000000
+30-AS_97_AK-SA2P0,0.0414265628,0.6400000000
+31-AS_97_AK-SA2P0,0.0315103928,0.6400000000
+32-AS_97_AK-SA2P0,0.0115621340,0.6400000000
+33-AS_97_AK-SA2P0,0.4468257966,0.6400000000
+34-AS_97_AK-SA2P0,0.3679796065,0.6400000000
+35-AS_97_AK-SA2P0,0.2383214103,0.6400000000
+36-AS_97_AK-SA2P0,0.1679059466,0.6400000000
+37-AS_97_AK-SA2P0,0.0792371032,0.6400000000
+38-AS_97_AK-SA2P0,0.0621908979,0.6400000000
+39-AS_97_AK-SA2P0,0.0417887618,0.6400000000
+40-AS_97_AK-SA2P0,0.0258313917,0.6400000000
+41-AS_97_AK-SA2P0,0.5932724009,0.6400000000
+42-AS_97_AK-SA2P0,0.5038029660,0.6400000000
+43-AS_97_AK-SA2P0,0.3494676791,0.6400000000
+44-AS_97_AK-SA2P0,0.2602197322,0.6400000000
+45-AS_97_AK-SA2P0,0.1382706601,0.6400000000
+46-AS_97_AK-SA2P0,0.1127585207,0.6400000000
+47-AS_97_AK-SA2P0,0.0806793144,0.6400000000
+48-AS_97_AK-SA2P0,0.0538092571,0.6400000000
+49-AS_97_AK-SA2P0,0.0967834918,0.7450000000
0-BJF_97-PGA,0.0784693168,0.4950000000
1-BJF_97-PGA,0.0726772362,0.4950000000
2-BJF_97-PGA,0.1021260341,0.4950000000
--
GitLab