We have produced a 5 year mean secular variation (SV) of the geomagnetic field for the period 2020-2025. We use the NASA Geomagnetic Ensemble Modeling System (GEMS), which consists of the NASA Goddard geodynamo model and ensemble Kalman filter (EnKF) with 512 ensemble members. Geomagnetic field models are used as observations for the assimilation, including gufm1 (1590-1960), CM4 (1961-2000) and CM6 (2001-2019). The forecast involves a bias correction scheme that assumes that the model bias changes on timescales much longer than the forecast period, so that they can be removed by successive forecasts. The algorithm was validated on the time period 2010-2015 by comparing with the 2015 IGRF before being applied to the 2020-2025 time period. This forecast has been submitted as a candidate model for IGRF 2025.