Recent data from the CDF Collaboration on the production of spin-triplet bottomonium states at the Fermilab Tevatron pp¯ collider are analyzed within the NRQCD factorization formalism. The color-singlet matrix elements are determined from electromagnetic decays and from potential models. The color-octet matrix elements are determined by fitting the CDF data on the cross sections for Υ(1S), Υ(2S), and Υ(3S) at large pT and the fractions of Υ(1S) coming from χb(1P) and χb(2P). We use the resulting matrix elements to predict the cross sections at the Tevatron for the spin-singlet states ηb(nS) and hb(nP). We argue that ηb(1S) may be observable in run II through the decay ηb-->J/ψ+J/ψ.