An efficient and stable computational scheme for parameterization of atmospheric chemistry is described. The 24-hour-average concentration of OH is represented as a set of high-order polynomials in variables such as temperature, densities of H2O, CO, O3, and NO1 (defined at NO+NO2+NO3+2N2O5+HNO2+HNO4) as well as variables determining solar irradiance: cloud cover, density of the overhead ozone column, surface albedo, latitude, and solar declination. This parameterization of OH chemistry was used in the three-dimensional study of global distribution of CH3CCl3 (Spivakovsky et al., this issue). The proposed computational scheme can be used for parameterization of rates of chemical production and loss or of any other output of a full chemical model. Coefficients for the polynomials are computed to provide the least squares fit to results of the full chemical model. Highly overdetermined systems are used with the sets of independent variables selected randomly in accordance with the distributions expected in the atmosphere. The least squares problem is solved using the Householder method of triangularization (by orthogonal transformations). The method allows detection and rectification of ill-defined conditions (i.e., linear dependence among terms), as well as evaluation of the individual contribution of each term the polynomial in reducing the residual vector. On the basis of that information the terms that have little bearing on the residual norm are discarded. Once the domain and the statistical distributions of independent variables are chosen, the entire parameterization procedure is implemented as a complete sequence of computer programs requiring no subjective analysis. The output of the procedure includes estimates of accuracy of the approximation against an independent sample of points, and computer written FORTRAN subroutines to compute the polynomials. ¿ American Geophysical Union 1990