AbstractWe present two new empirical models of radiation belt electron flux at geostationary orbit. GOES‐15 measurements of 0.8 MeV electrons were used to train a Nonlinear Autoregressive with Exogenous input (NARX) neural network for both modeling GOES‐15 flux values and an upper boundary condition scaling factor (BF). The GOES‐15 flux model utilizes an input and feedback delay of 2 and 2 time steps (i.e., 5 min time steps) with the most efficient number of hidden layers set to 10. Magnetic local time, Dst, Kp, solar wind dynamic pressure, AE, and solar wind velocity were found to perform as predicative indicators of GOES‐15 flux and therefore were used as the exogenous inputs. The NARX‐derived upper boundary condition scaling factor was used in conjunction with the Versatile Electron Radiation Belt (VERB) code to produce reconstructions of the radiation belts during the period of July–November 1990, independent of in‐situ observations. Here, Kp was chosen as the sole exogenous input to be more compatible with the VERB code. This Combined Release and Radiation Effects Satellite‐era reconstruction showcases the potential to use these neural network‐derived boundary conditions as a method of hindcasting the historical radiation belts. This study serves as a companion paper to another recently published study on reconstructing the radiation belts during Solar Cycles 17–24 (Saikin et al., 2021, https://doi.org/10.1029/2020sw002524), for which the results featured in this paper were used.