Joint distribution on the WL bias and scatter parameters $(b_l^{\text{WL}},s_l=\ln {({\sigma }_l^{\text{WL}})}^2)$, and their mass slopes ($b_{\mathrm{M}}^{\mathrm{WL}}$, σ_M), along with the principal sources of systematic uncertainty: the lensing efficiencies Σ_crit, l to the snapshots l, in units of 10¹⁵ h M_⊙ Mpc⁻², and the scale of the intrinsic mis-centering σ_intr. The high redshift biases anticorrelate with the respective lensing efficiencies, indicating that in this regime, we are limited by the DES Y3 photo-z accuracy. Conversely, the low redshift biases and the scatter values correlate among themselves, as they are dominated by the hydro uncertainty floor we added. Insert: Redshift interpolation of the WL lensing bias in cyan (filled region 1 sigma, faded lines 2 sigma), together with the 1 sigma uncertainty due to hydro modeling (green), explorations of the mean impact of the multiplicative shear bias m modeling (red), and the additional omission of the high mis-centering tail f = 0 (dashed red), and the posteriors of the WL bias in the four snap-shots. Hydrodynamical modeling uncertainties are dominant at low redshift. When considering our multiplicative shear bias modeling choices, the bias is consistent with zero, while the presence of the highly mis-centered tail does not affect the WL bias.