The Rhodes Gyre is a cyclonic persistent feature of the general circulation of the Levantine Basin in the eastern Mediterranean Sea. Although it is located in the most oligotrophic basin of the Mediterranean Sea, it is a relatively high primary production area due to strong winter nutrient supply associated with the formation of Levantine Intermediate Water. In this study, a 3D coupled hydrodynamic-biogeochemical model (SYMPHONIE/Eco3M-S) was used to characterize the seasonal and interannual variability of the Rhodes Gyre’s ecosystem and to estimate an annual organic carbon budget over the 2013–2020 period. Comparisons of model outputs with satellite data and compiled in situ data from cruises and BioGeoChemical-Argo floats revealed the ability of the model to reconstruct the main seasonal and spatial biogeochemical dynamics of the Levantine Basin. The model results indicated that during the winter mixing period, phytoplankton first progressively grow sustained by nutrient supply. Then, short episodes of convection driven by heat loss and wind events, favoring nutrient injections, organic carbon export, and inducing light limitation on primary production, alternate with short episodes of phytoplankton growth. The estimate of the annual organic carbon budget indicated that the Rhodes Gyre is an autotrophic area with a positive net community production in the upper layer (0–150 m) amounting to 31.2 ± 6.9 g C m-2 year-1. Net community production in the upper layer is almost balanced over the seven-year period by physical transfers, (1) via downward export (16.8 ± 6.2 g C m-2 year-1) and (2) through lateral transport towards the surrounding regions (14.1 ± 2.1 g C m-2 year-1). The intermediate layer (150–400 m) also appears to be a source of organic carbon for the surrounding Levantine Sea (7.5 ± 2.8 g C m-2 year-1) mostly through the subduction of Levantine Intermediate Water following winter mixing. The Rhodes Gyre shows high interannual variability with enhanced primary production, net community production, and exports during years marked by intense heat losses and deep mixed layers. However, annual primary production appears to be only partially driven by winter vertical mixing. Based on our results, we can speculate that future increase of temperature and stratification could strongly impact the carbon fluxes in this region.