  An estimate of the NNLO supersymmetric QCD effects for Higgs production at
hadron colliders is given. Assuming an effective gluon-Higgs interaction, these
corrections enter only in terms of process-independent, factorizable terms. We
argue that the current knowledge of these terms up to NLO is sufficient to
derive the NNLO hadronic cross section within the limitations of the standard
theoretical uncertainties arising mainly from renormalization and factorization
scale variations. The SUSY contributions are small with respect to the QCD
effects, which means that the NNLO corrections to Higgs production are very
similar in the Standard Model and the MSSM.
