Saturated soil thermal conductivity (λsat) is the maximum soil thermal conductivity value of a given soil. Although it can be determined accurately with a heat pulse sensor, there are challenges to prepare fully saturated soil samples. Numerous models have been developed to estimate λsat, and among these, the geometric mean method (GMM) generally performs well. The GMM requires soil mineral composition or quartz content information, which is unavailable for most soils. Earlier studies commonly used assumed that quartz content (fquartz) was equal to sand content (fsand) or to 0.5 × fsand, which led to significant λsat estimation errors especially on coarse-textured soils. We derived a novel method to estimate λsat from soil porosity (ϕ) based on a combination of the GMM and differential effective medium theory (DEM). The new DEM-GMM approach has a single parameter, cementation exponent (m). Using a calibration dataset of 43 soils, we determined best fit m values for soils in three groups: 1.66 for Group I (fsand < 0.4), 1.62 for Group II (0.4 <= fsand < 1) and m = -1.34ϕ+1.70 for Group III (fsand = 1). Using best fit m values for different groups, the new model can estimate λsat values from ϕ. Independent validation results on another 46 soils showed that the new model outperformed the GMM method with the assumption that fquartz = fsand or fquartz = 0.5 × fsand. The mean RMSE, Bias and R2 values of the DEM-GMM approach were 0.202 W m−1 K−1, 0.013 W m−1 K−1 and 0.89, respectively, and corresponding values of the GMM with the two assumptions were 0.295 and 0.476 W m−1 K−1, 0.056 and -0.28 W m−1 K−1, 0.80 and 0.82, respectively. The robust performance of the DEM-GMM approach suggests that it can be incorporated into thermal conductivity models to accurately estimate the thermal conductivity of unsaturated soils.