**Abstract**

The presence of spatial magma heterogeneities in volcanic monogenetic fields is a major observation discussed as well synthesized for worldwide volcanic fields. Magma heterogeneities still have not been visualized in the form of detailed spatial analyst tools, which could further help structuring works of geological mapping, volcanic hazard, and geoheritage evaluations. Here we synthetized 32 published datasets with a novel geochemical mapping model inspired by sub-disciplines of geomatic in one of the most documented monogenetic fields on earth: the Chichinautzin Volcanic Field (CVF) in Mexico. The volcanic units from CVF are covering the 2500 km2 area, and its neighbor stratovolcanoes are bordering the limit of most volcanic centers (Popocatepetl, Iztaccihuatl, and Nevado de Toluca). The results illustrate polygons and point map symbols from geochemical markers such as Alkalis vs SiO2, Sr/Y, and Ba/Nb. The geochemical heterogeneity of the CVF monogenetic bodies decreases as it approaches the Popocatepetl-Iztaccihuatl stratovolcanoes. This alignment is not observed in the occidental CVF portion near the flank of Nevado de Toluca, but geochemical anomalies associated to markers of continental crust interaction such as Sr/Y follow elongated patterns that are not strictly following structural lines and faults mapped on surface.

**Keywords:** monogenetic, spatial interpolation, trans-Mexican Volcanic Belt, geochemistry, Chichinautzin
