A 3D geological raster model has been constructed of the onshore of the Netherlands. The model displays geological units for the upper 500 m in 3D in an internally consistent way. The units are based on the lithostratigraphical classification of the Netherlands. This classification is used to interpret a selection of boreholes from the national subsurface database. Additional geological information regarding faults, the areal extent of each unit and conceptual genetic models have been combined in an automated workflow to interpolate the basal surfaces of each unit on 100 × 100 metre (x,y dimensions) raster cells. The combination of all interpolated basal surfaces results in a 3D Digital Geological Model (DGM) of the subsurface. A measure of uncertainty of each of these surfaces is also given. The automated workflow ensures an easily updatable subsurface model. The outputs are available for end users through www.dinoloket.nl.