The codes for simulating surface load Love numbers to higher-degree using VSMS method and calculating load Green's functions, and the results for five Earth models

Here, we provide a zipped file ‘codes&datasets’, containing two directories, 'codes' and ‘datasets’. Firstly, in the 'codes' directory, we provide the codes for simulating surface load Love numbers (LLNs) up to higher-degree using the VSMS method and calculating load Green’s functions (LGFs), with the PREM Earth model as an example. Corresponding input/output datasets are included. Secondly, in order to support load deformation studies for global, oceanic, and continental regions‌, the new reference model REM1D and the crustal model ECM1 are combined to construct ‌three spherically layered elastic Earth models‌—‌’globalc‌’, ‘‌oceanc’‌, and ‌’landc‌’—each incorporating the ‌average crustal physical structure‌ of its respective region. Using the VSMS method together with ‌fourth-order Runge-Kutta integration, we computed LLNs up to degree 1000000 for the three regional models and two reference models, REM1D and PREM (both with ocean-removed) and derived the corresponding LGFs for ten geophysical quantities: surface radial displacement, tangential displacement, geoid change, gravity perturbation, gravity change, vertical deflection, tilt change, radial strain, co-latitudinal strain, and longitudinal strain. Higher-degree LLNs pass asymptotic analytical solution tests. In the ‘datasets’ directory, we release comprehensive datasets for the three regionally adapted models (globalc, oceanc, landc) and two reference models (REM1D and PREM), including surface LLNs from degree 0 to 1000000 and LGFs for ten physical quantities under CE, CM, CF, CL, and CH reference frames. These datasets for the globalc, oceanc and landc models provide ‌new, complete, and high-quality fundamental parameters‌ for surface load deformation research across global, oceanic, and continental scales. The LLNs and LGFs dataset from the REM1D model can provide a valuable reference for the three aforementioned regional models. The results for two reference models (REM1D and PREM) serve as benchmark data for some loading simulation studies.