This study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Yuyao Maternal and Child Health Hospital (Yuyao Second People’s Hospital) (approval No. YMCH-EC-2023KT-0012). All participants provided written informed consent.
Study Participants
Primiparous women with spontaneous onset of labor who successfully achieved vaginal delivery were included. This enabled precise recording of the second-stage duration required for the predictive model. The eligibility criteria were as follows: (1) singleton pregnancy with a gestational age of ≥37 weeks; (2) cephalic presentation; and (3) informed consent obtained from both the participant and her family. The exclusion criteria included the following: (1) absolute contraindications to trial of labor identified during the active phase (e.g., placenta previa, vasa previa, cord prolapse); and/or (2) secondary uterine inertia (defined as failure of cervical dilation or fetal descent to progress for ≥4 h despite oxytocin augmentation or amniotomy). Women undergoing induction of labor were also excluded.
Notably, women who required emergency surgical intervention (C-section or instrumental delivery) due to acute fetal distress during the second stage were managed according to clinical guidelines but were not included in this analysis in order to preserve the completeness of the labor duration records.
The sample size was calculated based on the results of a preliminary trial (n = 30, in which the incidence of second-stage prolongation was 20%). With α set at 0.05 and β at 0.20, assuming an estimated incidence of second-stage prolongation of approximately 20%, a minimum of 150 participants was required to detect a significant predictor with an odds ratio (OR) ≥ 2.5. A total of 159 women were included in the final analysis.
Acquisition of Clinical Data
Maternal variables included age, gestational body mass index (BMI) increase (defined as the difference between pre-pregnancy BMI and pre-delivery BMI), gestational diabetes, gestational hypertension, mental state (assessed using the Self-Rating Anxiety Scale, with a score ≥50 indicating anxiety), uterine inertia (specifically referring to a clinically assessed hypotonic contraction pattern, defined as uterine contractions occurring at intervals >5 min or lasting <30 s, distinguished from the American College of Obstetricians and Gynecologists [ACOG] definition of labor arrest), and parity. Neonatal variables included birth weight (kg), fetal distress (defined as recurrent late or variable decelerations on cardiotocography or a baseline fetal heart rate of <110 beats/min or >160 beats/min), and Apgar score (with a 1-min score of <7 classified as low).
Gestational diabetes was diagnosed according to the International Association of Diabetes and Pregnancy Study Groups criteria (fasting glucose 5.1–6.9 mmol/L or 1 h ≥10.0 mmol/L or 2 h 8.5–11.0 mmol/L on the 75 g Oral Glucose Tolerance Test at 24–28 weeks). Gestational hypertension was defined according to ACOG guidelines (systolic blood pressure ≥140 mmHg or diastolic blood pressure ≥90 mmHg measured twice ≥4 h apart after 20 weeks of gestation, without proteinuria). Fetal distress on cardiotocography was defined as recurrent late decelerations (≥50% of contractions), recurrent variable decelerations lasting >60 s, prolonged deceleration lasting 2–10 min, or a baseline fetal heart rate of <110 or >160 beats/min for ≥10 min (National Institute for Health and Care Excellence [NICE] criteria). Uterine inertia (hypotonic contractions) was defined as contraction intervals >5 min or contraction duration of <30 s over 30 min, as assessed by palpation or tocodynamometry. Mental state (anxiety) was assessed using the Self-Rating Anxiety Scale at admission (cervical dilation < 4 cm), with a raw score of ≥50 (standardized index) classified as abnormal. A low Apgar score was defined as a 1-min score of <7 (WHO).
Because this was an observational study, labor management followed routine clinical practice without a standardized protocol. However, epidural analgesia (yes/no, timing, duration), oxytocin augmentation (yes/no, total dose, duration), artificial rupture of membranes (yes/no, cervical dilation at the procedure), and mode of delivery (spontaneous vaginal delivery only, according to the inclusion criteria) were systematically recorded. No significant interaction between these variables and the primary model predictors was observed.
Transabdominal and Transperineal Ultrasound Examination at the Early First Stage of Labor
Ultrasound scanning was conducted at the beginning of the active phase of labor (cervical dilation ≥4 cm). The parturient was placed in the supine position with a 15° left lateral tilt after bladder emptying. For transperineal scanning, the thighs were slightly abducted and flexed. All images were captured at the peak of a uterine contraction, which was identified by clinical palpation in conjunction with external tocodynamometry (cardiotocography) to confirm the peak of each contraction. Each measurement was repeated during three separate contractions, with the mean value used for analysis. Ultrasound examinations were performed using a portable diagnostic ultrasound system (see Table of Materials). For transabdominal imaging, a convex probe (2–5 MHz) was used; for transperineal imaging, a high-frequency linear or small-footprint convex probe (5–8 MHz) was used. Imaging settings were standardized (depth 6–10 cm, gain optimized, and a single focal zone positioned at the pubic symphysis).
For fetal position assessment, a transabdominal transverse scan identified the fetal orbit, cerebellum, and thalamic midline relative to the maternal pelvis to determine occiput anterior (OA), transverse, or posterior position. For transperineal measurements, the transducer was covered with a sterile sheath and placed on the perineum in the midsagittal plane without excessive pressure. Anatomical landmarks included the entire pubic symphysis, the fetal skull contour, and the anorectal junction.
The following parameters were measured (Figure 1A–D): angle of progression (AOP), defined as the angle between the long axis of the pubic symphysis and a line extending from its inferior edge to the deepest point of the fetal skull; midline angle (MLA), defined as the angle between the fetal skull midline (falx cerebri) and the maternal pelvic anteroposterior axis; head–symphysis distance (HSD), defined as the distance from the inferior edge of the pubic symphysis to the nearest point on the fetal skull; and head–perineum distance (HPD), defined as the shortest distance from the outer fetal skull to the perineal skin7.

Figure 1. Ultrasound measurements obtained during the first stage of labor. (A) Measurement of the angle of progression (AOP), defined as the angle between the long axis of the pubic symphysis and a line extending from its inferior edge to the deepest point of the fetal skull. (B) Measurement of the midline angle (MLA), defined as the angle between the fetal skull midline (falx cerebri) and the maternal pelvic anteroposterior axis. (C) Measurement of the head–symphysis distance (HSD), defined as the distance from the inferior edge of the pubic symphysis to the nearest point on the fetal skull. (D) Measurement of the head–perineum distance (HPD), defined as the shortest distance from the outer surface of the fetal skull to the maternal perineal skin. Please click here to view a larger version of this figure.
Each measurement was performed three times during three separate contraction peaks, and the average value was used for analysis. All examinations were conducted by three experienced sonographers (each with ≥10 years of experience in obstetric ultrasound), who completed standardized training in intrapartum transperineal ultrasound measurements before the study. Competency was confirmed by achieving an intraclass correlation coefficient (ICC) of ≥0.90 against a gold-standard assessor (senior author). Inter-observer and intra-observer reproducibility were assessed during a pilot phase involving 10 patients (not included in the main study), yielding ICCs ranging from 0.89 to 0.96 for all parameters (AOP, HSD, HPD, and MLA).
Image acceptance required clear visualization of the pubic symphysis, fetal skull contour, and anorectal junction. Images with excessive transducer pressure or poorly visualized anatomical landmarks were discarded and reacquired. Sonographers were blinded to group assignment (normal versus prolonged second stage) and to clinical data during image acquisition and measurement. Data were double-entered with range checks, and out-of-range values triggered image re-review. Missing data (<5%) were handled using complete-case analysis after confirming that the data were missing at random (Little’s missing completely at random test, p > 0.05).
Recording and Grouping of Labor Duration
Throughout the labor process, the course of labor for each participant was carefully documented. The first stage of labor was defined as the period from the onset of regular uterine contractions to complete cervical dilation (10 cm). The active phase of labor was defined as cervical dilation ≥4 cm accompanied by intensified, regular uterine contractions. The second stage of labor was defined as the period from complete cervical dilation (10 cm) to complete fetal expulsion. Ultrasound examinations were performed during the active phase of the first stage of labor to assess fetal head position, fetal position, and related ultrasound parameters. According to the World Health Organization (WHO) guidelines, a normal second stage of labor should not exceed 120 min2. To maintain a standardized research threshold for all primiparas in this study, participants were categorized according to this 120-min criterion into two groups: a normal second-stage group (≤120 min) and a prolonged second-stage group (>120 min).
Statistical Analysis
All statistical analyses were performed using SPSS version 26.0 (RRID:SCR_002865) and R version 4.1.0 (RRID:SCR_001905) (see Table of Materials). The following R packages were used: glmnet (version 4.1-2) for least absolute shrinkage and selection operator (LASSO) regression, rms (version 6.2-0) for nomogram construction and calibration, pROC (version 1.18.0) for receiver operating characteristic (ROC) analysis, and rmda (version 1.6) for decision curve analysis (DCA). Normality of continuous variables was assessed using the Shapiro–Wilk test to guide the choice of descriptive statistics and comparison methods. Normally distributed data were presented as mean ± standard deviation (x̄ ± s) and compared using the independent-samples t-test, which is robust to mild deviations from normality when sample sizes are adequate. Non-normally distributed data were presented as the median with interquartile range (M [P25, P75]) and compared using the Mann–Whitney U test, which evaluates whether one distribution tends to have larger values than the other rather than strictly testing differences in medians. Categorical variables were presented as counts and percentages (n [%]) and compared using the chi-square (χ2) test or Fisher’s exact test when the expected frequency was <5. All tests were two-tailed, and a p-value <0.05 was considered statistically significant. Variables with p <0.1 in the univariate analysis were entered into the LASSO regression for variable selection. LASSO regression was performed using the glmnet package with 10-fold cross-validation to determine the optimal regularization parameter (λ). The one-standard-error criterion (the largest λ within one standard error of the minimum cross-validation error) was applied to select the most parsimonious set of predictors with non-zero coefficients. Before LASSO fitting, all continuous predictors (maternal age, gestational BMI increase, AOP, HSD, HPD, and neonatal birth weight) were centered and standardized to z-scores using the default glmnet setting (standardize = TRUE) so that coefficients were penalized equally regardless of measurement units. The selected variables were subsequently entered into a multivariable logistic regression model (enter method) to calculate adjusted odds ratios (aORs) and 95% confidence intervals (CIs). Model significance was assessed using the likelihood ratio test.
Based on the logistic regression coefficients, a predictive model was constructed as Logit(P) = β0 + Σ(βiXi). This model was visualized as a nomogram using the rms package, in which each predictor was assigned a point value proportional to its regression coefficient, and the total points were mapped to predicted probabilities. Discrimination was evaluated by plotting the ROC curve and calculating the area under the curve (AUC) with its 95% CI using the pROC package (bootstrap method with 1,000 resamples). The optimal cutoff probability was determined by maximizing the Youden index (sensitivity + specificity − 1), with the corresponding sensitivity and specificity reported. Calibration was assessed using a calibration curve generated from 1,000 bootstrap resamples (rms package), supplemented by the Hosmer–Lemeshow goodness-of-fit test (a non-significant p-value indicates good calibration). Clinical utility was evaluated using DCA with the rmda package, which calculates the net benefit across a range of threshold probabilities (0%–100%). Net benefit was defined as (true positives / n) − (false positives / n) × (threshold probability / [1 − threshold probability]). The overall study workflow is described above and illustrated in Supplementary Figure 1. This model is exploratory; the predictor neonatal birth weight refers to actual weight measured after delivery, and clinical prediction would require model refitting using estimated fetal weight (EFW) obtained before delivery.