We present measurements of the inclusive production cross sections of the Upsilon(1S) bottomonium state in ppbar collisions at sqrt(s) = 1.96 TeV. Using the Upsilon(1S) to mu+mu- decay mode for a data sample of 159 +- 10 pb^-1 collected by the D0 detector at the Fermilab Tevatron collider, we determine the differential cross sections as a function of the Upsilon(1S) transverse momentum for three ranges of the Upsilon(1S) rapidity: 0 < |y| < 0.6, 0.6 < |y| < 1.2, and 1.2 < |y| < 1.8.
The proton-nucleon cross section ratio $R=Br(\Upsilon\to l^+l^-) d\sigma(\Upsilon)/dy|_{y=0} / {\sigma(J/\psi)}$ has been measured with the HERA-B spectrometer in fixed-target proton-nucleus collisions at 920 GeV proton beam energy corresponding to a proton-nucleon cms energy of sqrt{s}=41.6 GeV. The combined results for the Upsilon decay channels Upsilon $\to e^+e^-$ and Upsilon $\to\mu^+\mu^-$ yield a ratio $R=(9.0 \pm 2.1) 10^{-6}$. The corresponding Upsilon production cross section per nucleon at mid-rapidity (y=0) has been determined to be $Br(\Upsilon\to{}l^+l^-) {d\sigma(\Upsilon)/dy}|_{y=0}= 4.5 \pm 1.1 $ pb/nucleon.