The interest over of the decomposition of sulfates has increased due to its capacity of generating clean H2 through the thermochemical cycles. Understanding the decomposition mechanism is relevant to future industrial design and applications. Moreover, the modeling of these processes gives the information needed to know how much energy is required for the occurrence of the reactions. Among the different reaction systems, it is observed a range of complexity, with the presence of intermediate phases, and multiple consecutive or simultaneous reactions. Therefore, the present work proposed to develop a modeling methodology for the thermal decomposition of sulfates systems with different complexity levels: aluminum sulfate, potassium alum, mixture of aluminum sulfate and potassium sulfate, zinc sulfate, and iron (II) sulfate. The experiments were performed using thermogravimetric analysis (TGA) to understand the decomposition stages and use the data in the modeling step. The developed model consisted of a system of differential equations to describe every reaction taking place in the decomposition. The kinetic parameters estimation was made by using particle swarm optimization. The results indicate that potassium sulfate catalyzes the decomposition of aluminum sulfate. In the case of zinc, the desulfation of anhydrous zinc sulfate occurs in two stages, with the presence zinc oxysulfate as an intermediate phase. Iron (II) sulfate also shows a complex decomposition system, as it first decomposes into iron (III) sulfate before it is completely converted into hematite. All the modeling results displayed an excellent agreement with the experimental data, with R2 values above 0.98 for all cases.