We present the numerical investigation of the fermionic two-body decays of squarks in the Minimal Supersymmetric Standard Model with complex parameters. In the analysis we particularly take into account the cosmological bounds imposed by WMAP data. We find that the phase dependences of the decay widths of the third family squarks, as well as those of the first and second families, are quite significant which can provide viable probes of additional CP sources. We plot the CP phase dependences for each fermionic two-body decay channel of squarksq i (i=1,2 , q=u,d;c,s;t,b) and speculate about the branching ratios and total (two-body) decay widths. PACS numbers: 14.80.Ly, 12.60.Jv Recently, in the light of this cosmological constraint an extensive analysis of the neutralino relic density in the presence of SUSY-CP phases has been given by Bélanger et al. [10].In this study we present the numerical investigation of the fermionic two-body decays of squarks in MSSM with complex SUSY parameters. Actually, we had performed a 3 In our calculation, we have used WMAP-allowed bands in the plane M 1 −ϕ which are based on 1st year data. Now the WMAP 3rd year data is also available [11], but the new WMAP + SDSS combined value for relic density of dark matter does not change the numerical results in Ref.[10], namely the WMAP-allowed bands. See "Note added" section of Ref. [10].