<p>This paper presents a study focused on near-surface geophysical exploration in a region characterized by significant karst features, specifically the Cenote Ring, located at the outer edge of the Chicxulub impact crater. The high density of near-surface sinkholes, locally known as ‘cenotes’, and cavities in the area introduces complex groundwater flow patterns and varying degrees of geological risk. To address these challenges, we propose a gravity inversion method to estimate the spatial distribution and geometry of near-surface karst features. Although gravity data acquisition is rapid and cost-effective, interpretation in karst terrains is complex due to the gravity response of such features typically being in the range of a few tenths of a milligal. To mitigate nonuniqueness and uncertainty in gravity inversion, we integrate interpretations from Electrical Resistivity Tomography (ERT), Ground-Penetrating Radar (GPR), and LiDAR data. The inversion process, using the Simulated Annealing (SA) method, incorporates these constraints to build a comprehensive model. The resulting model is analyzed to identify geological patterns and structures associated with karstification, providing information on the dynamics of the subsurface groundwater flow of the region and enhancing the assessment of geological risk.</p>