diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/MotazedianAtkinson_2005.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/MotazedianAtkinson_2005.java
index 1024432d595a8884f15f83b1b8ddee51dc7b39e9..6cb0affe6c0ffe402c209c6eda88fcec17921de6 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/MotazedianAtkinson_2005.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/MotazedianAtkinson_2005.java
@@ -242,14 +242,14 @@ public class MotazedianAtkinson_2005 implements GroundMotionModel {
       // Get μ from MA_05_BASE for current IMT (interpolate if necessary)
       double μ = GroundMotions.combine(super.calc(in)).mean();
 
-      double μPga = calcMean(coeffsPGA, in.rRup, in.Mw);
+      double pgaRef = Math.exp(calcMean(coeffsPGA, in.rRup, in.Mw));
       double site = siteAmp.isPresent()
-          ? siteAmp.get().siteAmp(μPga, in.vs30)
-          : calcInterpolatedSite(super.coeffs.imt, μPga, in.vs30);
+          ? siteAmp.get().siteAmp(pgaRef, in.vs30)
+          : calcInterpolatedSite(super.coeffs.imt, pgaRef, in.vs30);
       return GroundMotions.createTree(μ + site, SIGMA);
     }
 
-    private static double calcInterpolatedSite(Imt imt, double μPga, double vs30) {
+    private static double calcInterpolatedSite(Imt imt, double pgaRef, double vs30) {
       Range<Imt> imtRange = INTERPOLATED_SITE_IMTS.get(imt);
       Imt imtLo = imtRange.lowerEndpoint();
       Imt imtHi = imtRange.upperEndpoint();
@@ -263,8 +263,8 @@ public class MotazedianAtkinson_2005 implements GroundMotionModel {
       double tTarget = imt.period();
       BooreAtkinson_2008 ba08lo = (BooreAtkinson_2008) BA_08_BASE.instance(imtLo);
       BooreAtkinson_2008 ba08hi = (BooreAtkinson_2008) BA_08_BASE.instance(imtHi);
-      double siteLo = ba08lo.siteAmp(μPga, vs30);
-      double siteHi = ba08hi.siteAmp(μPga, vs30);
+      double siteLo = ba08lo.siteAmp(pgaRef, vs30);
+      double siteHi = ba08hi.siteAmp(pgaRef, vs30);
       return Interpolator.findY(tLo, siteLo, tHi, siteHi, tTarget);
     }
   }