As multichannel equipment is increasingly used, we developed two algorithms for multivariable wavelet analysis of GPR signals (multi-cross wavelet MCW and Fourier multi-cross function FMC) and applied them to analyze raw GPR traces of archeological subsurface strata. The traces were from the tunnel located beneath the Temple of the Feathered Serpent (The Citadel, Teotihuacan, Mexico). The MCW and FMC algorithms determined the periods of subsurface strata of the tunnel. GPR traces inside-and-outside the tunnel/chamber, outside the tunnel/chamber and inside the tunnel/chamber analyzed with the MCW and filtered FMC algorithms determined the periods of the tunnel and chamber fillings, clay and matrix (limestone-clay compound). The tunnel filling period obtained by MCW analysis (14.37 ns) reflects the mixed limestone-clay compound of this stratum since its value is close to that of the period of the matrix (15.22 ns); periods of the chamber filling (11.40 ns) and the matrix (11.40 ns) were almost identical. FMC analysis of the tunnel obtained a period (5.08 ns) close to that of the chamber (4.27 ns), suggesting the tunnel and chamber are filled with similar materials. The use of both algorithms allows a deeper analysis since the similarities of the tunnel and chamber filling periods could not have been determined with the MCW algorithm alone. The successful application of the new multi-cross algorithms to archeological GPR data suggests they may also be used to search water and other resources in celestial bodies.