We demonstrate that CP violation in the Higgs sector, e.g. of a multi-doublet model, can be directly probed using gluon-gluon collisions at the SSC.Understanding the Higgs sector is one of the fundamental missions of future high energy colliders such as the SSC and LHC. In particular, it will be important to know if CP violation is present in the Higgs sector. Generally, either spontaneous or explicit CP violation can be present if the Higgs sector consists of more than the single doublet field of the Standard Model (SM). (For a review of this and other issues summarized below, see Ref.[1], and references therein.) However, important classes of models with extended Higgs sectors either do not allow for Higgs sector CP violation or are inconsistent with current experiment if significant CP violation in the Higgs sector is present. Among such models, supersymmetric theories are the most important example. There, a phase for a Higgs field vacuum expectation value in excess of about 10 −2 would imply imaginary components for slepton, squark, chargino and neutralino propagators that would result in electric dipole moments of the electron and neutron in excess of experimental limits. Thus, once a Higgs boson is discovered, it will be crucial to determine whether or not it is a pure CP eigenstate.Although there are a variety of experimental observables that are indirectly sensitive to CP violation in the Higgs sector (such as EDM's, top quark production and decay distributions, etc.), CP-violating contributions typically first appear at one-loop, or are otherwise suppressed, and will be very difficult to detect in a realistic experimental environment. In addition, if CP violation in this class of observables is detected, it could easily arise from sources other than the Higgs sector. In this letter, we shall show that the CP nature of a neutral Higgs boson (φ) is directly probed by the difference between its production rates through gluon-gluon fusion processes for colliding proton beams of opposite polarizations. (The proposed asymmetry is closely analogous to that developed previously for collisions of polarized back scattered laser beams at a future linear e + e − collider.[2] ) We compute the ⋆