A high dimensional model representation (HDMR) technique is introduced to capture the input-output behavior of chemical kinetic models. The HDMR expresses the output chemical species concentrations as a rapidly convergent hierarchical correlated function expansion in the input variables. In this paper, the input variables are taken as the species concentrations at time t(i) and the output is the concentrations at time t(i) + 6, where delta can be much larger than conventional integration time steps. A specially designed set: of model runs is performed to determine the correlated functions making up the HDMR. The resultant HDMR can be used to (i) identify the key input variables acting independently or cooperatively on the output, and (ii) create a high speed fully equivalent operational model (FEOM) serving to replace the original kinetic model and its differential equation salver. A demonstration of the HDMR technique is presented for stratospheric chemical kinetics. The FEOM proved to give accurate and stable chemical concentrations out to long times of many years. In addition, the FEOM was found to be orders of magnitude faster than a conventional stiff equation solver. This computational acceleration should have significance in many chemical kinetic applications.