This study estimates the Colombian territory's Curie Point Depth using magnetic anomaly data and two inversion methods, providing crucial insight into the crust's thermal structure. New national magnetic information acquired allows the data complementation with EMAG2; thus, we elaborate updates to the Curie isotherm distribution for northwestern South America. The CPD is calculated using the centroid and fractal nonlinear adjustment curve methods by moving square windows of 295 km and then compared with main tectonic and geologic features, geothermal gradient data from the National Hydrocarbon Agency, and element concentrations in active sediment. CPD average value is 43.3 ± 140.6 km for the centroid method and 20.1 ± 9.1 km for the fractal nonlinear adjustment method. Furthermore, shows shallower CPD values in the Southern Andes, coinciding with areas of volcanic density and geothermal anomaly areas like Las Animas - Chiles, Cerro Bravo - Machin or Huila - Sucubun. Unlikely, deeper values coincide with the Northern Andean sector, mainly in the Magdalena Basin and the Caldas Tear zone, the Paipa - Iza geothermal regime, and aging geologic bodies. Those deep values around the Caldas Tear have spatial correlations with the convergence of inferred Nazca, South American, and Caribbean plates due to the supposed subduction angle change. Furthermore, the CPD shows a main trend NE-SW probably controlled by the orogenic system of the three Andean Cordilleras. Lastly, the correlation between the CPD values and the element concentrations in active sediments shows moderate concordance for some elements, such as Ca or K.