This study investigates the solubility of neptunium (Np) in the deep natural groundwater of the Korea Atomic Energy Research Institute Underground Research Tunnel (KURT). According to a Pourbaix diagram (pH-$E_h$ diagram) that was calculated using the geochemical modeling program PHREEQC 2.0, the redox potential and the carbonate ion concentration both control the solubility of neptunium. The carbonate effect becomes pronounced when the total carbonate concentration is higher than $1.5;{ imes};10^{-2}$ M at $E_h$ = -200 mV and the pH value is 10. Given the assumption that the solubility-limiting stable solid phase is $Np(OH)_4(am)$ under the reducing condition relevant to KURT, the soluble neptunium concentrations were in the range of $1;{ imes};10^{-9}$ M to $3;{ imes};10^{-9}$ M under natural groundwater conditions. However, the solubility of neptunium, which was calculated with the formation constants of neptunium complexes selected in an OECD-NEA TDB review, strongly deviates from the value measured in natural groundwater. Thus, it is highly recommended that a prediction of neptunium solubility is based on the formation constants of ternary Np(IV) hydroxo-carbonato complexes, even though the presence of those complexes is deficient in terms of the characterization of neptunium species. Based on a comparison of the measurements and calculations of geochemical modeling, the formation constants for the "upper limit" of the Np(IV) hydroxo-carbonato complexes, namely $Np(OH)_y(CO_3)_z^{4-y-2z}$, were appraised as follows: log $K^{circ}_{122};=;-3.0{pm}0.5$ for $Np(OH)_2(CO_3)_2^{2-}$, log $K^{circ}_{131};=;-5.0{pm}0.5$ for $Np(OH)_3(CO_3)^-$, and log $K^{circ}_{141};=;-6.0{pm}0.5$ for $Np(OH)_4(CO_3)^{2-}$.