The determination of high-resolution geoid for marine regions requires the integration of gravity data provided by different sources, e.g. global geopotential models, satellite altimetry, and shipborne gravimetric observations. Shipborne gravity data, acquired over a long time, comprises the short-wavelengths gravitation signal. This paper aims to produce a consistent gravity field over the Red Sea region to be used for geoid modelling. Both, the leave-one-out cross-validation and Kriging prediction techniques were chosen to ensure that the observed shipborne gravity data are consistent as well as free of gross-errors. A confidence level equivalent to 95.4% was decided to filter the observed shipborne data, while the cross-validation algorithm was repeatedly applied until the standard deviation of the residuals between the observed and estimated values are less than 1.5 mGal, which led to the elimination of about 17.7% of the shipborne gravity dataset. A comparison between the shipborne gravity data with DTU13 and SSv23.1 satellite altimetry-derived gravity models is done and reported. The corresponding results revealed that altimetry models almost have identical data content when compared one another, where the DTU13 gave better results with a mean and standard deviation of -2.40 and 8.71 mGal, respectively. A statistical comparison has been made between different global geopotential models (GGMs) and shipborne gravity data. The Spectral Enhancement Method was applied to overcome the existing spectral gap between the GGMs and shipborne gravity data. EGM2008 manifested the best results with differences characterised with a mean of 1.35 mGal and a standard deviation of 11.11 mGal. Finally, the least-squares collocation (LSC) was implemented to combine the shipborne gravity data with DTU13 in order to create a unique and consistent gravity field over the Red Sea with no data voids. The combined data were independently tested using a total number of 95 randomly-chosen shipborne gravity stations. The comparison between the extracted shipborne gravity data and DTU13 altimetry anomalies before and after applying the LSC revealed that a significant improvement is procurable from the combined dataset, in which the mean and standard deviation of the differences dropped from -3.60 and 9.31 mGal to -0.39 and 2.04 mGal, respectively.