Physical and Biogeochemical Model Data for the Northern Benguela Upwelling System (1980–2020)

This dataset provides the monthly output fields from the coupled physical–biogeochemical modeling system developed by Salama et al. (2025). The modeling framework consists of the “Nucleus for European Modelling of the Ocean” (NEMO v4.2.2; Madec & the NEMO System Team, 2023) offline coupled to the “Biogeochemical Flux Model” (BFM v5.3; Vichi et al., 2023). The original simulation covers the larger Benguela Upwelling System (BUS) at high resolution (∼7 km), embedded within a global NEMO configuration at ∼25 km resolution and implemented through a two-way nesting strategy. The coupled NEMO–BFM system is forced by a combination of atmospheric, oceanic, and riverine inputs. Atmospheric fluxes of momentum, heat, and freshwater are computed using the ECMWF bulk formulation forced by ERA5 reanalysis fields (Hersbach et al., 2020). Both the global and regional configurations use spatially explicit geothermal heat fluxes from Goutorbe et al. (2011). Riverine freshwater discharge is prescribed from the GloFAS reanalysis (Alfieri et al., 2013), while climatological river nutrient loads are provided by the NEWS2 model (Mayorga et al., 2010). Initial and boundary conditions for temperature and salinity are obtained from the WOA18 climatologies (Locarnini et al., 2018; Zweng et al., 2019). Open boundary conditions for biogeochemical tracers are taken from the CMCC-ESM2 historical and ssp245 simulations (Lovato et al., 2023a,b). For the present dataset, the original model output has been spatially subsetted to the Northern Benguela Upwelling System (NBUS), extending horizontally from 7°E to 14°E and from 15.7°S to 21°S, and vertically from the surface to 1000 m depth. The dataset spans January 1980 to December 2020, providing a 41-year monthly time series of physical, biogeochemical variables in netCDF format. The Zenodo repository contains the following compressed folders (all in .tarformat): • BFM.tar — Biogeochemical variables from the BFM model, including 4D fields of dissolved oxygen (O2o), gross primary production (ruPPYc), net primary production (netPPYc), and total community respiration of phytoplankton (resPPYc), bacteria (resPBAc), and zooplankton (resZOOc). • NBUS_BFM_SPINUP.tar — Spin-up results of the simulated DO in the BFM configuration, over the NBUS domain. • temp_salinity.tar — 4D fields of potential temperature (toce) and practical salinity (soce) from the NEMO model. • ARMOR3D_NBUS.tar — 4D fields of potential temperature (toce) and practical salinity (soce) from the ARMOR3D observations. • CARS2009.tar — Climatological DO data from the CARS2009 observations. • vaisala.tar — Squared Brunt–Väisälä frequency (bn2) derived from NEMO. • MLD.tar — Mixed layer depth (mldkz5) defined as the turbocline depth. • velocity.tar — Three-dimensional velocity components (uo, vo, wo) representing zonal, meridional, and vertical currents. • volume_transport.tar — Corresponding volume transports (uocetr_eff, vocetr_eff, wocetr_eff) expressed in m3 s-1. • wind.tar — Zonal and meridional ERA5 wind speed components at 10 m elevation, interpolated onto the native NEMO model grid. • grid.tar — Contains the domain grid and the mesh-mask .nc file. The included Jupyter notebook (Oxygen.ipynb) in Oxygen_Jupyter_Notebook.tar contains the full workflow used to generate the figures presented in the manuscript submitted to Journal of Geophysical Research: Oceans. Note: *.tarfiles here are tar-only . Extract with tar -xf filename.tar This dataset is released under the Creative Commons Attribution 4.0 International (CC-BY 4.0) license. List of Variables: Symbol Long Name Unit BFM.tar O2o Dissolved Oxygen Concentration mmol O2 m-3 ruPPYc Gross Primary Production mg C m-3 d-1 netPPYc Net Primary Production mg C m-3 d-1 resPPYc Phytoplankton Respiration mg C m-3 d-1 resZOOc Zooplankton Respiration mg C m-3 d-1 resPBAc Bacterial Respiration mg C m-3 d-1 jsurO2o Surface Air-Sea Oxygen Flux mmol O2 m-2 d-1 NBUS_BFM_SPINUP.tar O2o Dissolved Oxygen Concentration mmol O2 m-3 temp_salinity.tar toce Potential Temperature °C soce Practical Salinity PSU vaisala.tar bn2 Squared Brunt–Väisälä Frequency s-2 MLD.tar mldkz5 Mixed Layer Depth m velocity.tar uo Zonal Current Velocity m s-1 vo Meridional Current Velocity m s-1 wo Vertical Current Velocity m s-1 volume_transport.tar uocetr_eff Zonal Volume Transport m3 s-1 vocetr_eff Meridional Volume Transport m3 s-1 wocetr_eff Vertical Volume Transport m3 s-1 wind.tar u10 ERA5 Zonal Wind Speed at 10 m m s-1 v10 ERA5 Meridional Wind Speed at 10 m m s-1