diff --git a/include/materials/GrainBoundaryCavitation.h b/include/materials/GrainBoundaryCavitation.h
index 7439da1..559f98f 100644
--- a/include/materials/GrainBoundaryCavitation.h
+++ b/include/materials/GrainBoundaryCavitation.h
@@ -77,7 +77,9 @@ class GrainBoundaryCavitation : public CZMComputeLocalTractionTotalBase
 
   /// @{ Initial conditions
   const MaterialProperty<Real> & _a0;
+  const MaterialProperty<Real> & _a0_old;
   const MaterialProperty<Real> & _b0;
+  const MaterialProperty<Real> & _b0_old;
   /// @}
 
   /// @{ Model parameters. See MooseDocs for their descriptions.
diff --git a/src/materials/GrainBoundaryCavitation.C b/src/materials/GrainBoundaryCavitation.C
index a3d6f23..e78bd18 100644
--- a/src/materials/GrainBoundaryCavitation.C
+++ b/src/materials/GrainBoundaryCavitation.C
@@ -82,7 +82,9 @@ GrainBoundaryCavitation::GrainBoundaryCavitation(const InputParameters & paramet
     _D_old(getMaterialPropertyOldByName<Real>(_base_name + "damage")),
     // initial conditions
     _a0(getMaterialPropertyByName<Real>(_base_name + getParam<MaterialPropertyName>("a0"))),
+    _a0_old(getMaterialPropertyOldByName<Real>(_base_name + getParam<MaterialPropertyName>("a0"))),
     _b0(getMaterialPropertyByName<Real>(_base_name + getParam<MaterialPropertyName>("b0"))),
+    _b0_old(getMaterialPropertyOldByName<Real>(_base_name + getParam<MaterialPropertyName>("b0"))),
     // model parameters
     _psi(getParam<Real>("psi") / 180 * M_PI),
     _D_GB(getMaterialPropertyByName<Real>(_base_name + getParam<MaterialPropertyName>("D_GB"))),
diff --git a/tests/materials/czm/GBCavitation-new/gold/Tx_150_Ty_0.e b/tests/materials/czm/GBCavitation-new/gold/Tx_150_Ty_0.e
index d77080f..fb1c894 100644
Binary files a/tests/materials/czm/GBCavitation-new/gold/Tx_150_Ty_0.e and b/tests/materials/czm/GBCavitation-new/gold/Tx_150_Ty_0.e differ