Soil boundary delineation is an important task in geotechnical site characterization. It can be achieved by either extracting borehole samples, conducting laboratory tests, and classifying them according to a soil classification system such as the Unified Soil Classification System (USCS) or utilizing multiple cone penetration test (CPT) soundings, and identifying soil boundaries at the soundings from the Ic (soil behavior type index) profiles. However, most soil-layer delineation methods can only take a single type of test result as the input. For instance, the well-known Markov random field (MRF) method can only take soil-type data such as sand, silt, or clay at boreholes as the input. Recognizing that soil classifications and soil properties are correlated, this paper proposes a novel coupled MRF-Bayesian framework to infer the spatial variation of USCS classifications (e.g., sand, silt, and clay) as well as soil properties by integrating both CPT and borehole data. This integrated approach leverages both CPT and borehole data to address some main challenges e.g., uncertainties and multivariate soil data input in underground stratification problems by simultaneous sampling of soil properties and soil types. The new unified framework can accommodate multivariate data, hence the new framework is compatible with the geotechnical engineering practice. The uncertainties for the spatial variation of USCS classification at sounding locations are quantified through a “layer-specific” Bayesian updating i.e., updating posterior cross-correlation behaviors for different layers (such as sand, silt, and clay), independently. In this Bayesian updating, soil-type data can provide some information about the soil properties according to the unified soil classification system. Further, the soil boundaries can be identified across the entire domain by the realization of conditional random fields of soil properties once the spatial variation of USCS classification is inferred at sounding locations, followed by a 3-dimensional Markov random field process.