We examine supergravity models with grand unification at M G possessing R parity invariance. Current data has begun to significantly constrain the parameter space. Thus for mSUGRA, accelerator data places a lower bound on m 1/2 of m 1/2 > ∼ 300 GeV while astronomical data on the amount of relic dark matter narrowly determines m 0 in terms of m 1/2 (for fixed value of tan β and A 0 ) due to co-annihilation effects. Additional new data could fix the parameters further. Thus the parameter space is sensitive to the muon magnetic moment anomaly, δa µ , and if δa µ lies 1σ above its current central value, it would exclude mSUGRA, while if it lies 1σ below (but is still positive) it pushes the SUSY spectrum into the TeV domain. The B s → µ + µ − decay is seen to be accessible to the Tevatron RunII with branching ratio sensitivity of Br[B s → µ + µ − ] > 6.5 × 10 −9 with 15fb −1 /detector, and a value of 7(14) × 10 −8 obtainable with 2fb −1 would be sufficient to exclude mSUGRA for tan β < 50(55). Measurements of B s → µ + µ − can cover the full mSUGRA parameter space for tan β > 40 if δa µ > 11 × 10 −10 , and combined measurements of B s → µ + µ − , a µ and m h (or alternately the gluino mass) would effectively determine the mSUGRA parameters for µ > 0. Detector cross sections are then within the range of planned future dark matter experiments. Non-universal models are also discussed, and it is seen that detector cross sections there can be much larger, and can be in the DAMA data region.