Increasing groundwater levels (GWLs) may become one of the most serious issues for the city of Odessa, Ukraine. This study investigated the spatial distribution characteristics and multifractal scaling behaviour of the groundwater-level/-depth fluctuations for a Quaternary aquifer in the city of Odessa using a geostatistical approach and multifractal detrended fluctuation analysis (MF-DFA). These two methods were applied to monthly GWL fluctuation time series from 1970 to 2020 to monitor 72 hydrogeological wells situated in different parts of the city of Odessa. The spatial distribution of the GWLs revealed an overall trend of decline and recovery from 1970 to 2020 in the study area, except for most of the southern region, where a persistent recovery of the groundwater depth was observed. The MF-DFA results suggest that the dynamics of the GWL fluctuations have multifractal characteristics in the Odessa area. In addition, both long-range correlations and fat-tail probability distribution contribute to the multifractality. However, long-range correlations among the fluctuations made a major contribution to the observed multifractality of the GWL fluctuation time series. The generalised Hurst exponents show a wide range of change (0.20 < h(q) < 2.85), indicating the sensitivity of the GWL fluctuations to changes in small-scale factors and large-scale factors. Regarding the long-range correlations of the GWL depths, the Hurst exponents (q = 2) demonstrated the positive persistence of groundwater-depth recovery in the southern region and the persistence of groundwater-depth variation in the other regions of the study area. The dynamic changes in the GWL depths in the Odessa area may be influenced by both natural and anthropogenic factors.