Modeling a very broad spectrum of the Earth's gravity field needs observations from various measurement techniques with different spectral sensitivities. Typically, high‐resolution regional gravity data are combined with low‐resolution global observations. To exploit the gravitational information as optimally as possible, we set up a regional modeling approach using radial spherical basis functions, emphasizing the strengths of various data sets by the flexible combination of high‐ and middle‐resolution terrestrial, airborne, shipborne, and altimetry measurements. The basis functions are defined and located in the region of interest in such a manner, which the highest measure of information of the input data is captured. Any functional of the Earth's gravity field can be derived, as, e.g., quasi‐geoid heights or gravity anomalies. Here we present results of a study area in Northern Germany. A comprehensive cross validation to external observation data delivers standard deviations less than 5 cm. Differences to an existing regional quasi‐geoid model count on average ±6 cm and proof the plausibility of our solution. The comparison with existing global models reaches higher standard deviations for the more sensitive gravity anomalies as for quasi‐geoid heights, showing the additional value of our solution in the high frequency domain. Covering a broad frequency spectrum, our regional models can be used as basis for various applications, such as refinement of global models, national geoid determination, and detection of mass anomalies in the Earth's interior.