本补充提供了用于原文中沿海暴露元分析方法论的重力模型的配置细节。它扩展了原始出版物中的描述，以确保我们结果的完全可重复性。具体配置细节如下所述，并在扩展数据表1-5中列出。重力模型信息来自国际全球地球模型中心的公开计算服务，GFZ亥姆霍兹地球科学中心（ICGEM）。全球范围内WGS84椭球面上重力高度异常的点数据（以米为单位），在EGM96、EGM2008和GOCO06s重力模型下，以0.085°的分辨率获得（数据下载日期分别为：EGM96和EGM2008为2023年8月22日，GOCO06s为2024年1月10日）。具体来说，ICGEM输出网格记录了最大纬度为90.000 °N和89.945 °S（2118个纬度平行线），最大经度为180 °W和179.975 °E（4236个经度平行线）。共下载了8971848个点用于每个重力模型。我们在重力模型配置中包含了ICGEM提供的零度项，并且对于每个重力模型，都参考了无潮汐系统配置。支撑这些重力模型的参考椭球WGS84由地心引力常数GM = 3.98600441800 × 10^14 m^3/s^2、半径a = 6378137.00 m、扁率f = 1/298.25722356300、角速度ω = 7.292115 × 10^–5 s^–1和正常势U0 = 6.263685171456948 × 10^7 m^2/s^2定义。我们在配置中未应用任何过滤器。因此，我们从ICGEM计算服务获得的重力模型数据包括EGM96（最大使用度360）特征值为加权均值−7.1416531 × 10^–3 m，最大值8.5891824 × 10^1 m，最小值−1.0651816 × 10^2 m以及加权均方根3.0566520 × 10^1 m；EGM2008（最大使用度2190）特征值为加权均值−8.2435599 × 10^–3 m，最大值8.6064395 × 10^1 m，最小值−1.0648962 × 10^2 m以及加权均方根3.0570532 × 10^1 m；以及GOCO06s（最大使用度300）特征值为加权均值−7.1516965 × 10^–3 m，最大值8.5419839 × 10^1 m，最小值−1.0651832 × 10^2 m以及加权均方根3.0569639 × 10^1 m。对于每个重力模型，高度异常点通过使用多重二次径向基函数进行插值，将其转化为全球栅格形式，从而提供最准确的插值结果，并且这也是引力场建模研究中常用的方法。插值是在ArcGIS Pro（版本3.1.3）的地质统计向导中进行的。每个插值的核参数都经过优化。标准邻域类型自动设置，邻居数最少10个，最多15个。使用了一个扇区类型为一扇区。对于我们的额外计算处理，包括GOCO2025s和EGM-DIR R4研究，它们在我们文章的后续处理中被添加。对于GOCO2025s和DIR R4重力模型，我们下载了与上述相同的WGS84椭球高度异常数据（两个重力模型的下载日期：2025年10月3日）。我们使用了1.000°的网格步幅和全球覆盖范围，从90.000 °N延伸到90.000 °S，从180.000 °W延伸到180.000 °E，因此获得了65341个网格点。与之前描述的重力模型相同，我们在重力模型配置中包含了ICGEM提供的零度项，并且对于每个重力模型都参考了无潮汐系统配置。因此，由ICGEM输出的配置包括GOCO2025s（最大使用度300）特征值为加权均值1.6370365 × 10^–2 m，最大值8.4607608 × 10^1 m，最小值−1.0634983 × 10^2 m，以及加权均方根3.0565538 × 10^1 m；和DIR R4（最大使用度260）特征值为加权均值1.6487188 × 10^–2 m，最大值8.4596036 × 10^1 m，最小值−1.0649776 × 10^2 m，以及加权均方根3.0565126 × 10^1 m。相应的重力模型点数据通过使用多重二次径向基函数进行插值到全球栅格，优化核参数，标准邻域类型最少10个，最多15个邻居，一个扇区类型。随后，通过减法计算得到了EGM96和GOCO06s重力模型的沿海偏差和平均动态地形。