  In this paper we use the conformal anomaly in QCD to derive the coupling of
light radion to gluons in the Randall-Sundrum model and use it to compute the
radion production cross section at hadron colliders by gluon fusion. We find
that the radion production cross section by gluon fusion at LHC would exceed
that of the higgs boson by a factor that lies between 7 and 8 over most of the
range. The decay modes of the radion are similar to that of the SM higgs boson.
But the striking feature is the enhancement of radion to 2-photon and radion to
2-gluon branching ratio over the SM case. Utilising this, we then discuss the
possible search strategies of such scalars at Tevatron and LHC. Using the
$\gamma \gamma$ decay mode one can explore/exclude radion mass upto 1 TeV. Even
with the current collected data at the Tevatron, one can exclude radion mass
upto 120 GeV for $\vphi$= 1 TeV.
