Frequency-domain electromagnetic induction (FDEMI) methods are frequently used in non-invasive, area-wise mapping of the subsurface electromagnetic soil properties. A crucial part of data analysis is the geophysical inversion of the data, resulting in either conductivity and/or magnetic susceptibility subsurface distributions. We present a novel 1D stochastic optimization approach that combines dimension-adapting reversible jump Markov chain Monte Carlo (MCMC) with artificial bee colony (ABC) optimization for geophysical inversion, with specific application to frequency-domain electromagnetic induction (FDEMI) data. Several solution models of simplified model geometry and a variable number of model knots, which are found by the inversion method, are used to create re-sampled resulting average models. We present synthetic test inversions using conductivity models based on 14 direct-push (DP) EC logs from Greece, Italy, and Germany, as well as field data applications using multi-coil FDEMI devices from three sites in Azerbaijan and Germany. These examples show that the method can effectively lead to solutions that resemble the known DP input models or image reasonable stratigraphic and archaeological features in the field data. Neighboring 1D solutions on field data examples show high coherence along profiles even though each 1D inversion is independently handled. The computational effort for one 1D inversion is less than 120,000 forward calculations, which is much less than usually needed in MCMC inversions, whereas the resulting models show more plausible solutions due to the dimension-adapting properties of the inversion method.