Abstract High-resolution ultrasonic imaging for defects in anisotropic multilayer carbon fiber reinforced polymers (CFRPs) is challenging because of the severe ultrasonic attenuation and the low signal-to-noise ratio (SNR) of echoes. The existing delay-multiply-and-sum (DMAS) beamforming outperforms delay-and-sum (DAS) beamforming in resolution, but with high computational complexity and energy loss. This paper presents a novel delay-sum-and-square (DSAS) beamforming algorithm. It takes full advantage of spatial coherence of captured data in the receiving and transmitting apertures. The non-coherent components caused by background noise are suppressed during the imaging. The back-wall reflection method (BRM) is used to correct the direction-dependent velocity. Full-matrix data is experimentally captured and processed on three different CFRP samples. Compared with DAS and DMAS, DSAS has a significant improvement in resolution, SNR and contrast. It demonstrates excellent defect characterization and noise suppression capability with only 17.4% computation time of DMAS.