The Alps are one of the best studied mountain ranges in the world, yet significant unknowns remain regarding their crustal structure and density distribution at depth. Previous published interpretations of crustal features within the orogen have been primarily based upon 2D seismic sections, and those that do integrate multiple geo-scientific datasets in 3D, have either focused on smaller sub-sections of the Alps or included the Alps, in low resolution, as part of a much larger study area. Therefore the generation of a 3D, crustal scale, gravity constrained, structural model of the Alps and their forelands at an appropriate resolution has been created here to more accurately describe crustal heterogeneity in the region. The study area of this work focuses on a region of 660 km x 620 km covering the vast majority of the Alps and their forelands are included, with the Central and Eastern Alps and the northern foreland being the best covered regions. Surface Generation All referenced data was integrated to constrain sub-surface lithospheric features including: previous regional scale gravitationally and seismically constrained models of the TRANSALP study area, the Upper Rhine Graben, the Mollasse Basin and the Po Basin; continental scale integrative best fit models (EuCRUST-07 and EPcrust); and seismic reflection depths from numerous published deep seismic surveys (e.g. ALP’75, EGT’86, ALP 2002 and EASI). The software package Petrel was used for the creation and visualisation of the modelled surfaces in 3D, representing the key structural and density contrasts within the region. All surfaces were generated with a grid resolution of 20 km x 20 km using Petrel’s convergent interpolation algorithm. During interpolation, a hierarchy of data source types was used in the case of contradiction between the different data sources and was as follows: 1. regional scale, gravitationally and seismically constrained models; 2. regional scale, seismically constrained models; 3. individual seismic reflection surfaces and interpreted sections; 4. continental scale, seismically constrained, integrative best fit models. No subduction interfaces were modelled. Topography and bathymetry were taken from ETOPO1 and the LAB from Geissler (2010). Gravity Modelling The generated surfaces and the calculated free-air anomaly from the global gravity model EIGEN-6C4, at a fixed height of 6 km above the datum were used in the 3D gravity modelling software IGMAS+ for the constrained of lithospheric density distribution. The layers of the generated model were split laterally into domains of different density, to reflect the heterogeneous nature of the crust within the region. Densities used in the initial structural model were derived from empirical P wave velocity to density conversions (Brocher, 2005) from the input seismic data sources. The densities were then modified, through multiple iterations, until the resulting model produced a gravity field within ± 25 mGal of the observed one. Surfaces generated as part of the integration work were not modified. Files The surface depths, thicknesses and densities of the model can be found as tab separated text files for each individual layer of the model (Unconsolidated Sediments, Consolidated Sediments, Upper Crust, Lower Crust and Lithospheric Mantle). The columns in each file are identical: the Easting is given in the “X COORD (UTM Zone 32N)”, Northing in the “Y COORD (UTM Zone 32N)”, the top surface depth of each layer is given as TOP (m.a.s.l), the thickness of each layer is given as THICKNESS (m), and the bulk density of that layer is given as DENSITY (Kg/m^3).