Following the adjustments of thermodynamic parameters according to DFT, LSRs, and GemNet, the kinetics in the model were adjusted to account for the change in the activation energies as the adsorbate enthalpies of formation were varied. A well-accepted method for this is to use BEP relations to linearly adjust the activation energy for a given reaction based on the enthalpy change of the reaction. Doing this after the intermediates' energies are estimated creates a dynamic MKM where changes in the species' thermodynamics have a realistic e↵ect on the activation energies. Unfortunately, the data for BEP relationships is often scarce for catalytic reactions because it relies on the calculation of multiple transition states, which can be a lengthy and expensive process for surface systems. Instead, this study used the Blowers-Masel approximation (BMA) to modify the reaction barriers, 42 as shown in Equation 1 and 2 below. The derivation of BMA parameters only needs the activation energy and reaction enthalpy of one reaction, circumventing issues related to data sparsity.