Objectives Atmospheric delay is one of the main error sources in interferometric synthetic aperture radar (InSAR) surface deformation monitoring. Using the global meteorological model to correct the tropospheric delay has been successfully applied in recent years.

Methods An InSAR atmospheric correction method is presented, which introduces China's first generation global atmospheric and land reanalysis product (CRA). Taking into account the physical characteristics of tropospheric meteorological parameters, the method was tested and evaluated quantitatively in Los Angeles, Southern California. The results show the feasibility of CRA products in InSAR atmospheric correction. First, the meteorological parameters such as temperature, water vapor and air pressure are obtained by CRA, taking into account the spatial variation characteristics of atmospheric parameters, according to the vertical stratification characteristics of meteorological data, the piecewise function is used to interpolate in the vertical direction. By considering the spatial variability of the atmosphere, the interpolation is carried out in the horizontal plane. According to the results, the atmospheric delay is calculated along the satellite line of sight. The atmospheric correction result is obtained by difference with the original interferogram. Second, verify and analyze the results. First, the error analysis is carried out from three aspects: interferogram standard deviation, spatial correlation and phase-elevation correlation coefficient. It is compared with the results of GACOS (generic atmospheric correction online service for InSAR) and PyAPS (Python⁃based atmospheric phase screen mitigation library) methods. Finally, 78 interferograms obtained from 40 scenes Sentinel-1 data of Southern California based on short baseline principle are used to carry out sequential deformation measurement experiments, the proposed method is compared with classical small baseline subset InSAR (SBAS-InSAR) results, and verified by global positioning system （GPS) data.

Results The results show that: (1) The average standard deviation of the corrected interferogram is reduced by 34.7% and the average standard deviation of 65% of interferograms has been reduced by more than 20%, especially for the interferograms that are seriously affected by the atmosphere. It is better than GACOS and PYAPS with an average improvement of 30%. (2) The significant decrease of the expected square variance of the spatial structure function of the interferogram shows that the product can effectively suppress the long-wave atmosphere. (3) The change of the fitting coefficient of correlation shows that the product can effectively reduce the atmospheric vertical stratification component caused by the influence of elevation. Depending on the degree of atmospheric influence, CRA can improve the vertical stratification component of the atmosphere, ranging from 20% to 60% or more. (4) It is effective to take the spatial variation of convection into account when establishing a high-resolution InSAR tropospheric delay map using low-resolution global atmospheric model data.

Conclusions This paper confirms that CRA can correct and improve the overall accuracy of deformation monitoring. Through quantitative evaluation, researchers can fully understand the correction effect and performance of this pro-duct, which is helpful in promoting the application and development of this product. However, due to the differences of external data, follow-up research can combine global navigation satellite system(GNSS) data with CRA solution to further improve the accuracy of InSAR atmospheric correction.