Decline of bone mineral density (BMD) during menopause is related to increased risk of fractures in postmenopausal women, however, this relationship in premenopausal women has not been established. To quantify this relationship, real‐world data (RWD) from the National Health and Nutrition Examination Survey (NHANES), and longitudinal data from the elagolix phase III clinical trials were modeled across a wide age range, and covariates were evaluated. The natural changes in femoral neck BMD (FN‐BMD) were well‐described by a bi‐exponential relationship with first‐order BMD formation (k1) and resorption (k2) rate constants. Body mass index (BMI) and race (i.e., Black) were significant predictors indicating that patients with high BMI or Black race experience a relatively lower BMD loss. Simulations suggest that untreated premenopausal women with uterine fibroids (UFs) from elagolix phase III clinical trials (median age 43 years [minimum 25–maximum 53]) lose 0.6% FN‐BMD each year up to menopausal age. For clinical relevance, the epidemiological FRAX model was informed by the simulation results to predict the 10‐year risk of major osteoporotic fracture (MOF). Premenopausal women with UFs, who received placebo only in the elagolix phase III trials, have a projected FN‐BMD of 0.975 g/cm2 at menopause, associated with a 10‐year risk of MOF of 2.3%. Integration of modeling, RWD, and clinical trials data provides a quantitative framework for projecting long‐term postmenopausal risk of fractures, based on natural history of BMD changes in premenopausal women. This framework enables quantitative evaluation of the future risk of MOF for women receiving medical therapies (i.e., GnRH modulators) that adversely affect BMD.