Subsurface engineering relies on sparsely distributed data points of the stress state of the earth's crust. 3D geomechanical--numerical modelling is applied to estimate the stress state in the entire volume of a large area. We present a multi-stage approach of differently sized models which provide the stress state in an area of interest derived from few and widely scattered data records. Furthermore we demonstrate the changes in reliability of the model depending on different input parameters.