From 954b7bb000a9c2ecf2e8e24022cb35f16a8d80a5 Mon Sep 17 00:00:00 2001
From: Peter Powers <pmpowers@usgs.gov>
Date: Mon, 9 Dec 2024 12:18:11 -0700
Subject: [PATCH] replaced bias with adjusted or similar
---
.../nshmp/gmm/AbrahamsonGulerce_2020.java | 30 +++++++--------
.../earthquake/nshmp/gmm/KuehnEtAl_2020.java | 38 +++++++++----------
.../earthquake/nshmp/gmm/ParkerEtAl_2020.java | 36 ++++++++----------
...sv => nga-sub-ak-interface-adjustment.csv} | 2 +-
4 files changed, 47 insertions(+), 59 deletions(-)
rename src/main/resources/gmm/coeffs/{nga-sub-ak-interface-bias-corr.csv => nga-sub-ak-interface-adjustment.csv} (96%)
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonGulerce_2020.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonGulerce_2020.java
index f8ad6f86..f4639525 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonGulerce_2020.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/AbrahamsonGulerce_2020.java
@@ -124,8 +124,8 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
static final CoefficientContainer COEFFS = new CoefficientContainer("AG20.csv");
static final CoefficientContainer COEFFS_PRVI = new CoefficientContainer("AG20_PRVI.csv");
- static final CoefficientContainer COEFFS_AK_BIAS =
- new CoefficientContainer("nga-sub-ak-interface-bias-corr.csv");
+ static final CoefficientContainer COEFFS_AK_ADJUSTED =
+ new CoefficientContainer("nga-sub-ak-interface-adjustment.csv");
private static final double C1S = 7.5;
private static final double C4 = 10.0;
@@ -147,8 +147,8 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
final double a1Adj;
final double ε; // epistemic uncertainty in mean
- // USGS bias correction coefficient for Gobal interface models in Alaska
- final double correctionAK;
+ // USGS data adjustment for NGA-Sub Gobal interface models in Alaska
+ final double ak_adjust;
Coefficients(
Imt imt,
@@ -215,11 +215,7 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
vlin = coeffs.get("vlin");
b = coeffs.get("b");
- /*
- * Average bias-correction for NGA-Subduction global interface models in
- * Alaska NSHM.
- */
- correctionAK = COEFFS_AK_BIAS.get(imt).get("bias_ak");
+ ak_adjust = COEFFS_AK_ADJUSTED.get(imt).get("adj_ak");
}
}
@@ -253,8 +249,8 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
return false;
}
- /* Apply USGS bias adjustment for Alaska */
- boolean akBiasAdjusted() {
+ /* Apply USGS Alaska data adjustment to global interface. */
+ boolean akAdjusted() {
return false;
}
@@ -269,11 +265,11 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
double pgaRock = 0.0;
if (in.vs30 < coeffs.vlin) {
pgaRock = exp(calcMean(
- coeffsPGA, slab(), usgsBasin(), adjusted(), akBiasAdjusted(),
+ coeffsPGA, slab(), usgsBasin(), adjusted(), akAdjusted(),
0.0, in.Mw, in.rRup, in.zTor, VS30_ROCK, Double.NaN));
}
double μ = calcMean(
- coeffs, slab(), usgsBasin(), adjusted(), akBiasAdjusted(),
+ coeffs, slab(), usgsBasin(), adjusted(), akAdjusted(),
pgaRock, in.Mw, in.rRup, in.zTor, in.vs30, in.z2p5);
double σ = calcSigma(coeffs, coeffsPGA, pgaRock, in.rRup, in.vs30);
@@ -290,7 +286,7 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
boolean slab,
boolean usgsBasin,
boolean adjust,
- boolean akBiasAdjusted,
+ boolean akAdjusted,
double pgaRock,
double Mw,
double rRup,
@@ -365,8 +361,8 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
μ += c.a1Adj;
}
- if (!slab && akBiasAdjusted) {
- μ += c.correctionAK;
+ if (!slab && akAdjusted) {
+ μ += c.ak_adjust;
}
return μ;
@@ -488,7 +484,7 @@ public abstract class AbrahamsonGulerce_2020 implements GroundMotionModel {
}
@Override
- boolean akBiasAdjusted() {
+ boolean akAdjusted() {
return true;
}
}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/KuehnEtAl_2020.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/KuehnEtAl_2020.java
index 33832fe6..e3a74453 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/KuehnEtAl_2020.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/KuehnEtAl_2020.java
@@ -147,8 +147,8 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
static final CoefficientContainer COEFFS = new CoefficientContainer("KBCG20.csv");
static final CoefficientContainer COEFFS_PRVI = new CoefficientContainer("KBCG20_PRVI.csv");
- static final CoefficientContainer COEFFS_AK_BIAS =
- new CoefficientContainer("nga-sub-ak-interface-bias-corr.csv");
+ static final CoefficientContainer COEFFS_AK_ADJUSTED =
+ new CoefficientContainer("nga-sub-ak-interface-adjustment.csv");
private static final double δM = 0.1; // p. 20
private static final double MREF = 6.0; // p. 20
@@ -188,8 +188,8 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
final double zBreak; // dzb_if dzb_slab
final double zRef;
- // USGS bias correction coefficient for Gobal interface models in Alaska
- final double correctionAK;
+ // USGS data adjustment for NGA-Sub Gobal interface models in Alaska
+ final double ak_adjust;
Coefficients(
Imt imt,
@@ -256,11 +256,7 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
Ï„ = coeffs.get("tau");
Ï = coeffs.get("rho");
- /*
- * Average bias-correction for NGA-Subduction global interface models in
- * Alaska NSHM.
- */
- correctionAK = COEFFS_AK_BIAS.get(imt).get("bias_ak");
+ ak_adjust = COEFFS_AK_ADJUSTED.get(imt).get("adj_ak");
}
}
@@ -306,8 +302,8 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
return false;
}
- /* Apply USGS bias adjustment for Alaska */
- boolean akBiasAdjusted() {
+ /* Apply USGS Alaska data adjustment to global interface. */
+ boolean akAdjusted() {
return false;
}
@@ -315,17 +311,17 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
public LogicTree<GroundMotion> calc(GmmInput in) {
double pgaRef = exp(calcPgaRef(
- coeffsPGA, slab(), akBiasAdjusted(),
+ coeffsPGA, slab(), akAdjusted(),
in.Mw, in.rRup, in.zTor));
double μ = calcMean(
- coeffs, slab(), usgsBasin(), seattle(), m9(), akBiasAdjusted(),
+ coeffs, slab(), usgsBasin(), seattle(), m9(), akAdjusted(),
in.Mw, in.rRup, in.zTor, in.vs30, in.z2p5, pgaRef);
// short periods can't be lower than PGA
if (coeffs.imt.isSA() && coeffs.imt.period() <= 0.1) {
double μPga = calcMean(
- coeffsPGA, slab(), usgsBasin(), seattle(), m9(), akBiasAdjusted(),
+ coeffsPGA, slab(), usgsBasin(), seattle(), m9(), akAdjusted(),
in.Mw, in.rRup, in.zTor, in.vs30, in.z2p5, pgaRef);
μ = max(μ, μPga);
}
@@ -347,7 +343,7 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
private static double calcPgaRef(
Coefficients c,
boolean slab,
- boolean akBiasAdjusted,
+ boolean akAdjusted,
double Mw,
double rRup,
double zTor) {
@@ -359,8 +355,8 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
double fSite = siteAmpTerm(c, VSROCK, 0);
double μ = c.θ1 + fMag + fGeom + fDepth + fAtten + fSite;
- if (!slab && akBiasAdjusted) {
- μ += c.correctionAK;
+ if (!slab && akAdjusted) {
+ μ += c.ak_adjust;
}
return μ;
@@ -372,7 +368,7 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
boolean usgsBasin,
boolean seattle,
boolean m9,
- boolean akBiasAdjusted,
+ boolean akAdjusted,
double Mw,
double rRup,
double zTor,
@@ -388,8 +384,8 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
double fBasin = basinTerm(c, vs30, z2p5, slab, usgsBasin, seattle, m9);
double μ = c.θ1 + fMag + fGeom + fDepth + fAtten + fSite + fBasin;
- if (!slab && akBiasAdjusted) {
- μ += c.correctionAK;
+ if (!slab && akAdjusted) {
+ μ += c.ak_adjust;
}
return μ;
@@ -611,7 +607,7 @@ public abstract class KuehnEtAl_2020 implements GroundMotionModel {
}
@Override
- boolean akBiasAdjusted() {
+ boolean akAdjusted() {
return true;
}
}
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java
index db5b4f48..d6bbcf22 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java
@@ -117,7 +117,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
static final CoefficientContainer COEFFS_INTERFACE;
static final CoefficientContainer COEFFS_SLAB;
- static final CoefficientContainer COEFFS_AK_BIAS;
+ static final CoefficientContainer COEFFS_AK_ADJUSTED;
static final CoefficientContainer COEFFS_PRVI;
static {
@@ -128,7 +128,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
*/
COEFFS_INTERFACE = new CoefficientContainer("PSBAH20_interface.csv");
COEFFS_SLAB = new CoefficientContainer("PSBAH20_slab.csv");
- COEFFS_AK_BIAS = new CoefficientContainer("nga-sub-ak-interface-bias-corr.csv");
+ COEFFS_AK_ADJUSTED = new CoefficientContainer("nga-sub-ak-interface-adjustment.csv");
COEFFS_PRVI = new CoefficientContainer("PSBAH20_PRVI.csv");
}
@@ -162,8 +162,8 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
final double Ï„, Ï•21, Ï•22, Ï•2v;
final double σε1, σε2, t1, t2;
- // USGS bias correction coefficient for Gobal interface models in Alaska
- final double correctionAK;
+ // USGS data adjustment for NGA-Sub Gobal interface models in Alaska
+ final double ak_adjust;
// inlined
// final double b4, v1, vRef;
@@ -259,11 +259,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
throw new IllegalArgumentException("Subduction zone [" + zone + "] not supported");
}
- /*
- * Average bias-correction for NGA-Subduction global interface models in
- * Alaska NSHM.
- */
- correctionAK = COEFFS_AK_BIAS.get(imt).get("bias_ak");
+ ak_adjust = COEFFS_AK_ADJUSTED.get(imt).get("adj_ak");
}
}
@@ -295,8 +291,8 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
return false;
}
- /* Apply USGS bias adjustment for Alaska */
- boolean akBiasAdjusted() {
+ /* Apply USGS Alaska data adjustment to global interface. */
+ boolean akAdjusted() {
return false;
}
@@ -304,10 +300,10 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
public LogicTree<GroundMotion> calc(GmmInput in) {
double pgaRef = exp(calcMean(
- coeffsPGA, slab(), akBiasAdjusted(), in));
+ coeffsPGA, slab(), akAdjusted(), in));
double μ = calcMean(
- coeffs, slab(), usgsBasin(), m9(), akBiasAdjusted(), pgaRef, in);
+ coeffs, slab(), usgsBasin(), m9(), akAdjusted(), pgaRef, in);
double Ï•Total = phiTotal(coeffs, in.rRup, in.vs30);
double σ = Maths.hypot(coeffs.τ, ϕTotal);
@@ -324,7 +320,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
private static double calcMean(
Coefficients c,
boolean slab,
- boolean akBiasAdjusted,
+ boolean akAdjusted,
GmmInput in) {
double h = nearSourceSaturation(in.Mw, c.mc, slab);
@@ -333,8 +329,8 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
double Fd = depthTerm(c, slab, in.zTor);
double μ = c.c0 + Fp + Fm + Fd;
- if (!slab && akBiasAdjusted) {
- μ += c.correctionAK;
+ if (!slab && akAdjusted) {
+ μ += c.ak_adjust;
}
return μ;
@@ -346,7 +342,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
boolean slab,
boolean usgsBasin,
boolean m9,
- boolean akBiasAdjusted,
+ boolean akAdjusted,
double pgaRef,
GmmInput in) {
@@ -363,8 +359,8 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
double μ = c.c0 + Fp + Fm + Fd + Fs;
- if (!slab && akBiasAdjusted) {
- μ += c.correctionAK;
+ if (!slab && akAdjusted) {
+ μ += c.ak_adjust;
}
return μ;
@@ -600,7 +596,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
}
@Override
- boolean akBiasAdjusted() {
+ boolean akAdjusted() {
return true;
}
}
diff --git a/src/main/resources/gmm/coeffs/nga-sub-ak-interface-bias-corr.csv b/src/main/resources/gmm/coeffs/nga-sub-ak-interface-adjustment.csv
similarity index 96%
rename from src/main/resources/gmm/coeffs/nga-sub-ak-interface-bias-corr.csv
rename to src/main/resources/gmm/coeffs/nga-sub-ak-interface-adjustment.csv
index 9879ac02..a485a326 100644
--- a/src/main/resources/gmm/coeffs/nga-sub-ak-interface-bias-corr.csv
+++ b/src/main/resources/gmm/coeffs/nga-sub-ak-interface-adjustment.csv
@@ -1,4 +1,4 @@
-T,bias_ak
+T,adj_ak
PGA,-0.26167
PGV,0.00833
0.01,-0.30351
--
GitLab