Camacho, AG, A.G.Montesinos, F. G.Vieira, Ricardo2023-06-202023-06-202000Bosshard, E., and MacFarlane,D.J., 1970, Crustal structure of theWestern Canary Islands from seismic refraction and gravity data: J. Geophys. Res., 75, 4901–4918. Camacho, A. G., Montesinos, F. G., and Vieira, R., 1997, A threedimensional gravity inversion applied to Sao Miguel Island (Azores): J. Geophys. Res., 102, 7705–7715. Moritz, H., 1980, Advanced physical geodesy: HerbertWichmannVerlag. Pedersen, L. B., 1979, Constrained inversion of potential field data: Geophys. Prosp., 27, 726–748. Pérez-Torrado, F.J., Marti, J., Mangas, J., and Day, S., 1997, Ignimbrites of the Roque Nublo group, Gran Canaria, Canary Islands: Bull. Volcanol., 58, 647–654. Pick, M., Picha, J., and Vysk ˆ ocil,V., 1973, Theory of the Earth’s gravity field: Elsevier Science Publ. Co. René, R. M., 1986, Gravity inversion using open, reject, and “shapeof- anomaly” fill criteria: Geophysics, 51, 988–994. Rothman, D. H., 1985, Nonlinear inversion, statistical mechanics, and residual statics estimation: Geophysics, 50, 2784–2796. Silva, J.B.C. and Hohmann,G.W., 1983, Nonlinear magnetic inversion using a random search method: Geophysics, 46, 1645–1658. Tarantola, A., 1988, The inverse problem theory: Methods for data fitting and model parameter estimation: Elsevier Science Publ. Co. Zidarov, D., 1990, Inverse gravimetric problem in geoprospecting and geodesy: Elsevier Science Publ. Co.0016-803310.1190/1.1444729https://hdl.handle.net/20.500.14352/57460This paper presents a gravity inversion method for determining the volumes of bodies with pre-established density contrasts. The method works step-by-step on a prismatic partition of the subsurface volume, expanding the anomalous bodies to fit the observed gravity values in a systematic exploration of model possibilities. The process is treated in a 3-D context; at the same time, it can determine a simple regional trend. Moreover, positive and negative density contrasts are simultaneously accepted. The solution is obtained by a double condition: (1) the e(2)-fitness to the observed gravity data (model fitness) and (2) the minimization of the total (weighted) anomalous mass (model smoothness). A positive parameter is used to balance the two minimization terms. The method is applied to a simulated example and also to a real example: the volcanic island of Gran Canaria (Canary Islands, Spain). In both cases, the results obtained show the possibilities of the method.engjournal articlehttp://geophysics.geoscienceworld.org/content/65/1/95.full.pdf+htmlhttp://geophysics.geoscienceworld.org/open access550.3Canary-IslandsGran-CanariaGeochemistry & GeophysicsGeofísica2507 Geofísica