diff --git a/include/cantera/kinetics/AqueousKinetics.h b/include/cantera/kinetics/AqueousKinetics.h
index 4b410950a..bec0b5a89 100644
--- a/include/cantera/kinetics/AqueousKinetics.h
+++ b/include/cantera/kinetics/AqueousKinetics.h
@@ -139,8 +139,6 @@ protected:
 
     std::vector<size_t> m_irrev;
 
-    std::vector<size_t> m_fwdOrder;
-
     size_t m_nirrev;
     size_t m_nrev;
 
diff --git a/include/cantera/kinetics/GasKinetics.h b/include/cantera/kinetics/GasKinetics.h
index d0c9a2999..be286a1de 100644
--- a/include/cantera/kinetics/GasKinetics.h
+++ b/include/cantera/kinetics/GasKinetics.h
@@ -123,8 +123,6 @@ protected:
     Rate1<Plog> m_plog_rates;
     Rate1<ChebyshevRate> m_cheb_rates;
 
-    std::vector<size_t> m_fwdOrder;
-
     size_t m_nirrev;
     size_t m_nrev;
 
diff --git a/src/kinetics/AqueousKinetics.cpp b/src/kinetics/AqueousKinetics.cpp
index 06d1064d6..ddba2ca73 100644
--- a/src/kinetics/AqueousKinetics.cpp
+++ b/src/kinetics/AqueousKinetics.cpp
@@ -56,8 +56,6 @@ AqueousKinetics& AqueousKinetics::operator=(const AqueousKinetics& right)
     m_nfall = right.m_nfall;
     m_rates = right.m_rates;
     m_irrev = right.m_irrev;
-
-    m_fwdOrder = right.m_fwdOrder;
     m_nirrev = right.m_nirrev;
     m_nrev = right.m_nrev;
     m_rxntype = right.m_rxntype;
@@ -357,9 +355,6 @@ void AqueousKinetics::addElementaryReaction(ReactionData& r)
 
     // add constant term to rate coeff value vector
     m_rfn.push_back(r.rateCoeffParameters[0]);
-
-    // forward rxn order equals number of reactants
-    m_fwdOrder.push_back(r.reactants.size());
 }
 
 void AqueousKinetics::installReagents(const ReactionData& r)
diff --git a/src/kinetics/GasKinetics.cpp b/src/kinetics/GasKinetics.cpp
index 437b145c0..8b3cfc337 100644
--- a/src/kinetics/GasKinetics.cpp
+++ b/src/kinetics/GasKinetics.cpp
@@ -68,7 +68,6 @@ GasKinetics& GasKinetics::operator=(const GasKinetics& right)
     m_plog_rates = right.m_plog_rates;
     m_cheb_rates = right.m_cheb_rates;
 
-    m_fwdOrder = right.m_fwdOrder;
     m_nirrev = right.m_nirrev;
     m_nrev = right.m_nrev;
     m_rrxn = right.m_rrxn;
@@ -499,10 +498,6 @@ void GasKinetics::addFalloffReaction(ReactionData& r)
     m_falloffn.install(m_nfall, r.falloffType, r.reactionType,
                        r.falloffParameters);
 
-    // forward rxn order equals number of reactants, since rate
-    // coeff is defined in terms of the high-pressure limit
-    m_fwdOrder.push_back(r.reactants.size());
-
     // increment the falloff reaction counter
     ++m_nfall;
 }
@@ -514,9 +509,6 @@ void GasKinetics::addElementaryReaction(ReactionData& r)
 
     // add constant term to rate coeff value vector
     m_rfn.push_back(r.rateCoeffParameters[0]);
-
-    // forward rxn order equals number of reactants
-    m_fwdOrder.push_back(r.reactants.size());
 }
 
 void GasKinetics::addThreeBodyReaction(ReactionData& r)
@@ -527,9 +519,6 @@ void GasKinetics::addThreeBodyReaction(ReactionData& r)
     // add constant term to rate coeff value vector
     m_rfn.push_back(r.rateCoeffParameters[0]);
 
-    // forward rxn order equals number of reactants + 1
-    m_fwdOrder.push_back(r.reactants.size() + 1);
-
     m_3b_concm.install(nReactions(), r.thirdBodyEfficiencies,
                        r.default_3b_eff);
 }
@@ -541,8 +530,6 @@ void GasKinetics::addPlogReaction(ReactionData& r)
 
     // add a dummy entry in m_rfn, where computed rate coeff will be put
     m_rfn.push_back(0.0);
-
-    m_fwdOrder.push_back(r.reactants.size());
 }
 
 void GasKinetics::addChebyshevReaction(ReactionData& r)
@@ -552,8 +539,6 @@ void GasKinetics::addChebyshevReaction(ReactionData& r)
 
     // add a dummy entry in m_rfn, where computed rate coeff will be put
     m_rfn.push_back(0.0);
-
-    m_fwdOrder.push_back(r.reactants.size());
 }
 
 void GasKinetics::installReagents(const ReactionData& r)