Background Bronchopulmonary dysplasia (BPD) is a significant cause of morbidity in preterm infants, yet its development and severity at high altitudes (>1500m) remain poorly understood. This study aimed to identify altitude-specific risk factors and develop robust, interpretable predictive models for BPD in this unique population.Methods In this retrospective matched cohort study, 378 preterm infants (<32 weeks gestation, <1500 g birth weight) admitted to a high-altitude (1500m) NICU between 2019 and 2023 were analysed. The cohort included 189BPD cases (91 mild, 61 moderate, 37 severe) and 189 matched controls. Maternal, perinatal and postnatal data were collected. Machine learning models (XGBoost, logistic regression, random forest) were developed and rigorously evaluated using comprehensive performance metrics to predict BPD occurrence and severity. SHAP (SHapley Additive exPlanations) analysis was employed to interpret the best-performing model. Results Key risk factors for BPD development included maternal hypertension (OR 2.31, 95%CI 1.56 to 3.42), initial oxygen requirement >30% (OR 3.15, 95%CI 2.13 to 4.65) and lack of exclusive breast milk feeding (OR 1.89, 95%CI 1.28 to 2.79). Severe BPD was independently associated with prolonged invasive ventilation (>7 days) (OR 4.12, 95%CI 2.78 to 6.11), elevated C reactive protein (>10mg/L) (OR 2.87, 95%CI 1.93 to 4.26) and patent ductus arteriosus (OR 2.53, 95%CI 1.71 to 3.74). Machine learning models demonstrated strong predictive performance; the optimal XGBoost model achieved an area under the curve of 0.89 (95% CI 0.85 to 0.93), an F1 score of 0.82, a Matthews Correlation Coefficient of 0.73 and a balanced accuracy of 0.85. SHAP analysis identified initial FiO2 >30%, mechanical ventilation and maternal hypertension as the top three most influential predictors for the XGBoost model.Conclusions This study provides the first comprehensive analysis of BPD risk factors at a specific high altitude and validates effective, interpretable machine learning models for its prediction. These findings highlight the critical importance of altitude-specific adjustments in risk assessment and emphasise the potential for modelguided early interventions to improve outcomes for this vulnerable population.