Implementation of the Gauss-Kronrod Quadrature Method (G7, K15) on 2D Gravity Anomaly Modeling in Basins with a Polynomial Variation of Density Distribution with Depth

Forward modeling of 2D gravity anomalies, considering density contrasts that vary polynomially with depth, was performed to examine basin structures. This process involved two main stages: deriving analytical formulas and executing numerical integration. The Gauss-Kronrod Quadrature Method, utilizing 7 Gauss points and 15 Kronrod points, was employed to precisely compute these integrals. Initial modeling applied to theoretical basement scenarios with fixed density contrasts showed gravity anomalies that accurately reflected the curvature of the basement. To validate the approach, it was then applied to real-world cases including the Sebastian Vizcaino Basin, San Jacinto Graben, and Sayula Basin. By incorporating suitable density contrasts, modeling lengths, and basement curvature shapes, the results revealed that both fixed-density and depth-variable density models produced gravity anomalies with patterns consistent with the actual basement curvature. These findings validate the modeling technique’s effectiveness in representing real geological features accurately. The study confirms that the Gauss-Kronrod Quadrature Method (G7, K15) is robust for analyzing 2D gravity anomalies, providing a reliable tool for understanding the influence of varying density contrasts on gravity responses.