Electric resistivity surveys and borehole collection of resistivity data are one of the oldest geophysical tools for characterization of the vadose zone. The current trend is to conduct such surveys in a tomographic manner, which requires significant computational resources. We present a simple semi-analytical approach to delineate multiple layers in partially saturated soils from resistivity and saturation measurements taken at several depths along a borehole. The number of layers and their hydraulic properties are assumed to be known. The proposed inversion algorithm is computationally efficient and can serve either as a stand-alone tool for layers delineation or as an autonomous module in a more comprehensive geophysical survey. It is most robust when each layer is sampled at least once. When one or more layers have not been sampled, the algorithm’s robustness (convergence) is dependent on the accuracy of an initial guess (e.g., expert knowledge and other hard and/or soft data). We provide a detailed analysis of the algorithm’s convergence and identify potential pitfalls.