The potential energy surfaces of one, two, and three water molecule sequential adsorptions on the symmetrically chlorinated Si(100)-2 x 1 surface were theoretically explored with SIMOMM:MP2/6-31G(d). The first water molecule adsorption to the surface dimer requires a higher reaction barrier than the subsequent second water molecule adsorption. The lone pair electrons of the incoming water molecule nucleophilically attack the surface Si atom to which the leaving Cl group is bonded, yielding an S(N)2 type transition state. At the same time, the Cl abstracts the H atom of the incoming water molecule, forming a unique four-membered ring conformation. The second water molecule adsorption to the same surface dimer requires a much lower reaction barrier, which is attributed to the surface cooperative effect by the surface hydroxyl group that can form a hydrogen bond with the incoming second water molecule. The third water molecule adsorption exhibits a higher reaction barrier than the first and the second water molecule adsorption channels but yields a thermodynamically more stable product. In general, it is expected that the surface Si-Cl bonds can be subjected to the substitution reactions by water molecules, yielding surface Si-OH bonds, which can be a good initial template for subsequent surface chemical modifications. However, oversaturations can be a competing side reaction under severe conditions, suggesting that the precise control of surface kinetic environments is necessary to tailor the final surface characteristics.