Long‐period zonal tidal Love numbers are more sensitive to lower‐mantle anelasticity (which is difficult to determine even from laboratory); however, the Love numbers of the zonal tides have not been adequately determined. In this study, we use the optimal sequence estimation method with the Global Positioning System (GPS) displacements to estimate the complex Love numbers of the long‐period zonal tides. For the first time, we retrieve the complex Love numbers h20 and l20 for the Msm, Msf, SN, Mstm, Mtm, Msqm, and Mqm tides from the Up‐ and North‐components of the global GPS displacements; we also re‐estimate the same complex Love numbers for the Mm (monthly) and Mf (fortnightly) tides with a higher precision. Based on the robust estimated h20 and l20 values for the nine zonal tides, we inverse their corresponding mantle anelastic parameters fr and fi, respectively. We finally obtain that the material‐dependent parameter α for the lower mantle is 0.189 ± 0.07, and find that the 18.6‐year nodal tide and the Chandler wobble, as well as the long‐period zonal tides, can be represented by a single absorption band. Our estimates for fr, fi, and α can be further used to construct attenuation and velocity models of the Earth's interior.