This paper considers an integrated sensing and communication system, where some radar targets also serve as communication scatterers. A location domain channel modeling method is proposed based on the position of targets and scatterers in the scattering environment, and the resulting radar and communication channels exhibit a two-dimensional (2-D) joint burst sparsity. We propose a joint scattering environment sensing and channel estimation scheme to enhance the target/scatterer localization and channel estimation performance simultaneously, where a spatially non-stationary Markov random field (MRF) model is proposed to capture the 2-D joint burst sparsity. Specifically, the base station (BS) first transmits downlink pilots to sense the targets in the scattering environment. Then the user transmits uplink pilots to estimate the communication channel. Finally, joint scattering environment sensing and channel estimation are performed at the BS based on the reflected downlink pilot signal and received uplink pilot signal. A message passing based algorithm is designed based on the turbo approach. Furthermore, a variation of the expectation maximization method, based on pseudo-likelihood approximation, is proposed to learn the unknown parameters in the non-stationary MRF model that capture the random sizes and locations of the scattering clusters. The advantages of our proposed scheme are verified in the simulations.