Not yet reviewed by Pith; the record is open.

The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for ~250 million unique stars. By combining BDBS photometry with the latest Gaia astrometry, we characterize the chemo-dynamics of red clump stars across the BDBS footprint, using an unprecedented sample size and sky coverage. We construct a sample of ~2.3 million red clump giants in the bulge with photometric metallicities, BDBS photometric distances, and proper motions. We study the kinematics of the red clump stars as a function of sky position and metallicity, by investigating proper motion rotation curves, velocity dispersions, and proper motion correlations across the southern Galactic bulge. We find that metal-poor red clump stars exhibit lower rotation amplitudes, at ~29 km s$^{-1}$ kpc^{-1}. The peak of the angular velocity is ~39 km s^{-1} kpc^{-1} for [Fe/H] ~ -0.2 dex, exhibiting declining rotation at higher [Fe/H]. The velocity dispersion is higher for metal-poor stars, while metal-rich stars show a steeper gradient with Galactic latitude, with a maximum dispersion at low latitudes along the bulge minor axis. Only metal-rich stars ([Fe/H] >~ -0.5 dex) show clear signatures of the bar in their kinematics, while the metal-poor population exhibits isotropic motions with an axisymmetric pattern around Galactic longitude l = 0. This work reports the largest sample of bulge stars with distance, metallicity, and astrometry and shows clear kinematic differences with metallicity. The global kinematics over the bulge agrees with earlier studies. However, we see striking changes with increasing metallicity and for the first time, see kinematic differences for stars with [Fe/H]>-0.5, suggesting that the bar itself may have kinematics that depends on metallicity.