The variational adjoint method was applied to retrieving the turbulivity of the atmospheric Ekman boundary layer along with the regularization principle. The validity of the method was verified by using the idealized data, and then the turbulivity profile and the geostrophic wind profile were retrieved through it for real observational wind filed data.