A better understanding of the spatiotemporal dynamics and influencing factors of sulfonamide antibiotic resistance genes (ARGs) distribution in subsurface flow constructed wetlands is essential to improve the ARGs removal efficiency. The spatiotemporal dynamics of sulfonamide ARGs were explored in the vertical upflow subsurface flow constructed wetland (VUSFCW). The results showed that the absolute abundance of ARGs presented a trend of bottom layer > middle layer > top layer. The relative abundance of ARGs decreased significantly from the bottom layer to the middle layer, but increased in the top layer. The bottom layer was the main stage to remove ARGs. The absolute abundance of ARGs at each point in summer was significantly higher than that in winter. Based on the spatiotemporal distribution of ARGs, the internal mechanism of ARGs dynamic change was explored by the partial least square path analysis model. The results showed that physical-chemical factors, microorganisms and antibiotics indirectly affected the spatiotemporal distribution of ARGs mainly through mobile genetic elements. The indirect influence coefficients of physical-chemical factors, microorganisms and antibiotics on the spatiotemporal distribution of ARGs were 0.505, 0.221 and 0.98 respectively. The direct influence coefficient of MGEs on the spatiotemporal distribution of ARGs was 0.895. The results of network analysis showed that the potential host species of ARGs in summer were more abundant than those in winter. The selection mode of sulfonamide ARGs to potential hosts was nonspecific. There is a risk of sulfonamide ARGs infecting pathogens in VUSFCW. Fortunately, VUSFCW has proven effective in reducing the absolute abundance of ARGs and the potential risk of pathogens carrying ARGs. These findings provide a model simulation and theoretical basis for effectively reducing the threat of ARGs.