The Central Andean orogeny is caused by the subduction of the Nazca oceanic plate beneath the South-American continental plate. In Particular, the Southern Central Andes (SCA, 27°-40°S) are characterized by a strong N-S and E-W variation in the crustal deformation style and intensity. Despite being the surface geology relatively well known, the information on the deep structure of the upper plate in terms of its thickness and density configurations is still scarcely constrained. Previous seismic studies have focused on the crustal structure of the northern part of the SCA (~27°-33°S) based upon 2D cross-sections, while 3D crustal models centred on the South-American or the Nazca Plate have been published with lower resolution.

To gain insight into the present-day state of the lithosphere in the area, we derived a 3D model that is consistent with both the available geological and seismic data and with the observed gravity field. The model consists on a continental plate with sediments, a two-layer crust and the lithospheric mantle being subducted by an oceanic plate. The model extension covers an area of 700 km x 1100 km, including the orogen, the forearc and the forelands.

Different data sets were integrated to derive the lithospheric features:

• We used the global relief model of ETOPO1 (Amante and Eakins 2009) for the topography and bathymetry.

• The sub-surface structures were defined by integrating seismically constrained models, including the South-American crustal thickness of Assumpção et al. (2013; model A; 0.5 degree resolution), the sediment thickness of CRUST1 (Laske et al. 2013) and the slab geometry of SLAB2 (Hayes et al. 2018).

• Additionally, we included seismic reflection and refraction profiles performed on the Chile margin (Araneda et al. 2003; Contreras-Reyes et al. 2008, 2014, 2015; Flueh et al. 1998; Krawzyk et al. 2006; Moscoso et al. 2011; Sick et al. 2006; Von Huene et al. 1997).

• Besides, we used sediment thickness maps from the intracontinental basin database ICONS (6 arc minute resolution, Heine 2007) and two oceanic sediment compilations: one along the southern trench axis (Völker et al. 2013) and another of global-scale (GlobSed; Straume et al. 2019).

To build the interfaces between the main lithospheric features, we compiled and interpolated these datasets on a regular grid with a surface resolution of 25 km. For that purpose, the convergent algorithm of the software Petrel was used. We assigned constant densities within each layer, except for the lithospheric mantle. In this case, we implemented a heterogeneous distribution by converting s-wave velocities from the SL2013sv seismic tomography (Schaeffer and Lebedev 2013) to densities. The python tool VelocityConversion was used for the conversion (Meeßen 2017).

To further constrain the crustal structure of the upper plate, a gravity forward modelling was carried out using IGMAS+ (Schmidt et al. 2010). The gravity anomaly from the model (calculated gravity) was compared to the free-air anomaly from the global gravity model EIGEN-6C4 (observed gravity; Förste et al 2014; Ince et al. 2019). Subsequently, the crystalline crust of the upper plate was split vertically into two layers of different densities. We inverted the residual between calculated and observed gravity to compute the depth to the interface between the two crustal layers. For the inverse modelling of the gravity residual, the Python package Fatiando a Terra was used (Uieda et al. 2013)

For each layer, the depth to the top surface, thickness and density can be found as separate files. All files contain identical columns:
- Northing as "X Coord (UTM zone 19S)";
- Easting as "Y Coord (UTM zone 19S)";
- depth to the top surface as "Top (m.a.s.l)" and
- thickness of each layer as "Thickness (m)".

The header ‘Density’ indicates the bulk density of each unit in kg/m3. For the oceanic and continental mantle units, a separate file is provided with a regular grid of the density distribution with a lateral resolution of 8 km x 9 km and a vertical resolution of 5 km. The containing columns are: Northing as "X Coord (UTM zone 19S)"; Easting as "Y Coord (UTM zone 19S)"; depth as "Depth (m.a.s.l)" and density as "Density (kg/m3)"