The SUDETES 2003 wide-angle refraction/reflection experiment covered the area of the south-western Poland and the northern Bohemian Massif. The good quality data that were gathered combined with the data from previous experiments (POLONAISE'97, CELEBRATION 2000) allowed us to prepare a 3D seismic model of the crust and uppermost mantle for this area. We inverted travel times of both refracted and reflected P waves using the JIVE3D package. This allowed us to obtain a model of P-wave velocity distribution as well as the shape of major boundaries in the crust. We also present a detailed uncertainty analysis for both the boundary depths and the velocity field. In doing the uncertainty analysis we found an interesting, strong dependence between uncertainty and inversion scheme (order of used phases). We also compared the model with surface geology and found good correlation between velocity inhomogeneities in the uppermost crust (down to 2 km) and major geological units. The higher velocity lower crust (6.9–7.2 km/s) could result from remelting of the lower crust or magmatic underplating.