Application of the belief propagation (BP) algorithm is proposed for pilot-assisted reception of multi-carrier code division multiple access systems where the input symbols are binary phase shift keying symbols. The proposed method utilizes a Bayesian network graph that includes channel parameters as explicit nodes representing random variables. Channel parameters in both frequency and time domains are included in the graph. The messages on the graph are obtained by making large system assumptions for the transmitted symbol messages and Gaussian assumptions for the messages related to the channel parameters. It is observed via simulations that the proposed approach achieves performance very close to that of the previously proposed BP method with known channel parameters.