Deriving the kinetic reaction constants associated with hydrogen diffusion and trapping in metals from thermal desorption mass spectroscopy (TDS) spectra proves to be complicated and the existing analysis methods are subject to debate. This article will provide a brief background of several commonly employed analysis techniques and discuss the necessity of a more complex and rigorous analysis method for the determination of the kinetic constants associated with hydrogen trapping interactions. Furthermore, a numerical simulation method will be proposed using the McNabb & Foster equations to fit experimental TDS spectra in order to derive both diffusion and trapping/detrapping parameters, including the respective pre-exponential constants and activation energies associated with these interactions in metals. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.