  We provide a systematic study of charmless $B_s \to PP, PV, VV$ decays ($P$
and $V$ denote pseudoscalar and vector mesons, respectively) based on an
approximate six-quark operator effective Hamiltonian from QCD. The calculation
of the relevant hard-scattering kernels is carried out, the resulting
transition form factors are consistent with the results of QCD sum rule
calculations. By taking into account important classes of power corrections
involving "chirally-enhanced" terms and the vertex corrections as well as weak
annihilation contributions with non-trivial strong phase, we present
predictions for the branching ratios and CP asymmetries of $B_s$ decays into
PP, PV and VV final states, and also for the corresponding polarization
observables in VV final states. It is found that the weak annihilation
contributions with non-trivial strong phase have remarkable effects on the
observables in the color-suppressed and penguin-dominated decay modes. In
addition, we discuss the SU(3) flavor symmetry and show that the symmetry
relations are generally respected.
