So far, the contribution of ammonia-oxidizing archaea (AOA) to ammonia oxidation in wastewater treatment processes has not been well understood. In this study, two soil aquifer treatment (SATs) systems were built up to treat synthetic domestic wastewater (column 1) and secondary effluent (column 4), accomplishing an average of 95% ammonia removal during over 550 days of operation. Except at day 322, archaeal amoA genes always outnumbered bacterial amoA genes in both SATs as determined by using quantitative polymerase chain reaction (q-PCR). The ratios of archaeal amoA to 16S rRNA gene averaged at 0.70 ± 0.56 and 0.82 ± 0.62 in column 1 and column 4, respectively, indicating that all the archaea could be AOA carrying amoA gene in the SATs. The results of MiSeq-pyrosequencing targeting on archaeal and bacterial 16S rRNA genes with the primer pair of modified 515R/806R indicated that Nitrososphaera cluster affiliated with thaumarchaeal group I.1b was the dominant AOA species, while Nitrosospira cluster was the dominant ammonia-oxidizing bacteria (AOB). The statistical analysis showed significant relationship between AOA abundance (compared to AOB abundance) and inorganic and total nitrogen concentrations. Based on the mathematical model calculation for microbial growth, AOA had much greater capacity of ammonia oxidation as compared to the specific influent ammonia loading for AOA in the SATs, implying that a small fraction of the total AOA would actively work to oxidize ammonia chemoautotrophically whereas most of AOA would exhibit some level of functional redundancy. These results all pointed that AOA involved in microbial ammonia oxidation in the SATs.