Here the authors developed a two-dimensional two-temperature chemical non-equilibrium (2T-NCE) model of Ar-CO2-H-2 inductively coupled thermal plasmas (ICTP) around atmospheric pressure (760 torr). Assuming 22 different particles in this model and by solving mass conservation equations for each particle, considering diffusion, convection and net production terms resulting from 198 chemical reactions, chemical non-equilibrium effects were taken into account. Species density of each particle or simply particle composition was also derived from the mass conservation equation of each one taking the non chemical equilibrium effect into account. Transport and thermodynamic properties of Ar-CO2-H-2 thermal plasmas were self-consistently calculated using the first order approximation of the Chapman-Enskog method at each iteration point implementing the local particle composition and temperature. Calculations at reduced pressure (500 and 300 torr) were also done to investigate the effect of pressure on non-equilibrium condition. Results obtained by the present model were compared with results from one temperature chemical equilibrium (1T-CE) model, one-temperature chemically non equilibrium (1T-NCE) model and finally with 2T-NCE model of Ar-N-2-H-2 plasmas. Investigation shows that consideration of non-chemical equilibrium causes the plasma volume radially wider than CE model due to particle diffusion. At low pressure with same input power, presence of diffusion is relatively stronger than at high pressure. Comparison of present reactive model with non-reactive Ar-N-2-H-2 plasmas shows that maximum temperature reaches higher in reactive C-H-O molecular system than non-reactive plasmas due to extra contribution of reaction heat.