Objectives: The objectives of this paper are to examine the source of multifractality in uterine electromyography (EMG) signals and to study the progression of pregnancy in the term (gestation period > 37 weeks) conditions using multifractal detrending moving average (MFDMA) algorithm. Methods: The signals for the study, considered from an online database, are obtained from the surface of abdomen during the second (T1) and third trimester (T2). The existence of multifractality is tested using Hurst and scaling exponents. With the intention of identifying the origin of multifractality, the preprocessed signals are converted to shuffle and surrogate data. The original and the transformed signals are subjected to MFDMA to extract multifractal spectrum features, namely strength of multifractality, maximum, minimum, and peak singularity exponents. Results: The Hurst and scaling exponents extracted from the signals indicate that uterine EMG signals are multifractal in nature. Further analysis shows that the source of multifractality is mainly owing to the presence of long-range correlation, which is computed as 79.98% in T1 and 82.43% in T2 groups. Among the extracted features, the peak singularity exponent and strength of multifractality show statistical significance in identifying the progression of pregnancy. The corresponding coefficients of variation are found to be low, which show that these features have low intersubject variability. Conclusion: It appears that the multifractal analysis can help in investigating the progressive changes in uterine muscle contractions during pregnancy. © 2013 IEEE.