输沙势计算公式为：\n
DPi=∑U 2×(U-Ut)/100×t\n
其中i表示16个方位，具体有N、NNE、NE、NEE、E、EES、ES、ESS、S、SSW、WS、WWS、W、WWN、NW和NNW；U表示的是在10m高度上大于或等于阈值风速的有效风的每个风速级中的平均风速；Ut表示的是参考粒径0.25和0.30mm对应的阈值风速5.3和5.8m/s，是10m高度上使沙子处于跃移状态的最小风速；而t是刮风的时间，风速统计表中以百分率表示。采用2m/s的间隔，对大于5.3m/s的有效风，采用的平均风速分别为6、8、10、12、14、16、18、20、22、24、26、28、30、32和33m/s；而对大于5.8m/s的有效风，采用的平均风速分别为7、9、11、13、15、17、19、21、23、25、27、29、31和32m/s。注意：为了确保Fryberger提出的风能环境的有效分类（低能<200VU，200VU≤中能<400VU，高能≥400VU），这些平均风速已经使用换算系数（1knot=0.5144m/s）将速度单位m/s表达成速度单位knot。\n
将16个风向上得到的输沙势投影到正东（RDP_E）或正北方位(RDP_N)上，目的是为了计算合成输沙势（RDP）、合成输沙方向（RDD）和风向变率（WDV）:\n
RDP=sqrt(RDP_E 2+RDP_N 2)\n
RDD=180+57.296atan2(RDP_E⁄RDP_N)\n
WDV=RDP/DP\n
注意计算有效风百分比时的分母使用的是平年8760，闰年8784（年总小时数）。\n
由10m高度U和V向风速得到的风速风向数据来自于1940-2025年0.25° ERA5小时风数据和1950-2025年0.1° ERA5-Land小时风数据。\n
DP为来自16个风向的输沙势总和；RDP为合成输沙势，表示各种方向风下的净输沙势；RDD为合成输沙方向，表示输沙净走向；WDV为风向变率，是RDP/DP的比率，表示风来自一个方向（接近1），风来自很多方向（接近0）。具体地说，DP_方位表示该风向上的输沙势；RDP_E或RDP_N是将16个风向上的输沙势投影到正东或正北方位上，目的是为了计算合成输沙势（RDP）、合成输沙方向（RDD）和最终的风向变率（WDV）；输沙势_风向_区间是该风向上2m/s区间的输沙势；方位_频数是该风向上2m/s区间内出现的起沙风次数；方位_频率是该风向上2m/s区间内的起沙风比例。所有数据以GeoTIFF格式存储。", "ds_source": "", "ds_process_way": "", "ds_quality": "", "ds_acq_start_time": "1940-01-01 00:00:00", "ds_acq_end_time": "2025-12-31 00:00:00", "ds_acq_place": "全球", "ds_acq_lon_east": 180.0, "ds_acq_lat_south": -90.0, "ds_acq_lon_west": -180.0, "ds_acq_lat_north": 90.0, "ds_acq_alt_low": null, "ds_acq_alt_high": null, "ds_share_type": "open-access", "ds_total_size": 152732675012, "ds_files_count": 0, "ds_format": "GeoTIFF", "ds_space_res": "0.1度, 0.25度", "ds_time_res": "年", "ds_coordinate": "无", "ds_projection": "", "ds_thumbnail": "deb6b5a1-8c9a-4ee9-969b-520ba4faa833.png", "ds_thumb_from": 0, "ds_ref_way": "", "paper_ref_way": "", "ds_ref_instruction": "None", "ds_from_station": null, "organization_id": "52b7b79b-860c-49a5-9083-9a70cf8bed5a", "ds_serv_man": null, "ds_serv_phone": null, "ds_serv_mail": null, "doi_value": "", "subject_codes": [ "170.1545", "170.4510" ], "quality_level": 0, "publish_time": "2026-07-13 16:13:22", "last_updated": "2026-07-13 16:13:22", "protected": false, "protected_to": null, "lang": "zh", "cstr": "11738.11.ncdc.desert.db7640.2026", "i18n": { "en": { "title": "A dataset of global sand drift", "ds_format": "GeoTIFF", "ds_source": "", "ds_quality": "", "ds_ref_way": "", "ds_abstract": "
The sand drift potential (DP, in vector units (VU)) is calculated by\r\n
DPi=∑U 2×(U-Ut)/100×t\r\n
where i represents the 16 wind directions (N, NNE, NE, NEE, E, EES, ES, ESS, S, SSW, WS, WWS, W, WWN, NW and NNW); U is the mean wind speed in each speed category with an effective wind speed greater than or equal to the threshold wind speed Ut at 10 m height; Ut = 5.3 and 5.8 m s-1 is the threshold wind speed corresponding to the reference diameter d = 0.25 and 0.30 mm, which is the minimum wind speed at the 10 m height to cause sand in saltation; t is time wind blew, expressed as a percentage in a wind summary. The 2 m s-1 bin is adopted in each direction of effective winds ≥5.3 m s-1 (≥5.8 m s-1), corresponding to mean wind speeds of 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32 and 33 m s-1 (7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31 and 32 m s-1), and the sum of these bins is the final DP in the wind direction. Note that these mean wind speeds have been expressed in knots by the conversion factor (1 knot = 0.5144 m s-1) to ensure the valid classification of wind energy environments (low energy, <200 VU; intermediate energy, ≥200 VU and <400 VU; high energy, ≥400 VU) developed by Fryberger. \r\n
The obtained DPs in the 16 wind directions are projected onto the due-east (RDP_E) and due-north (RDP_N) directions to calculate the RDP, RDD and WDV:\r\n
RDP = sqrt(RDP_E 2 + RDP_N 2)\r\n
RDD = 180 + 57.296atan2(RDP_E ⁄ RDP_N)\r\n
WDV = RDP / DP\r\n
Notably, the divisor used in calculating the percentage of effective winds from different directions is the total hour number of Julian years (8760 hours for common years or 8784 hours for leap years). \r\n
The wind speed and direction data calculated by 10 m u-component of wind and 10 m v-component wind were from the 0.25° ERA5 hourly wind data during 1940-2025 and the 0.1° ERA5-Land hourly wind data during 1950-2025. \r\n
The drift potential (DP) represents the sum of DPs from the 16 wind directions. The resultant drift potential (RDP) indicates net sand transport potential under different wind directions, while the resultant drift direction (RDD) reflects the net trend of sand drift. Wind directional variability (WDV) is defined as the ratio of RDP to DP, indicating whether winds predominately come from a single direction (approaching 1) or from multiple directions (approaching 0). Specifically, DP_Azimuthal refers to the drift potential from a given wind direction. RDP_E and RDP_N represent the projections of the DPs from the 16 wind directions onto the due-east and due-north directions, respectively, which are used to calculate the RDP, RDD and the final WDV. DP_azimuthal_bin is the DP in the 2 m s-1 bin of the wind direction. Azimuthal_absolute frequency represents the count of sand-moving winds in the 2m s-1 bin of the wind direction; azimuthal_relative frequency represents the frequency of sand-moving winds in the 2 m s-1 bin of the wind direction. All data are stored in GeoTIFF format.", "ds_time_res": "", "ds_acq_place": "global", "ds_space_res": "", "ds_projection": "", "ds_process_way": "", "ds_ref_instruction": "" } }, "submit_center_id": "ncdc", "data_level": 0, "recommendation_value": 0, "license_type": "https://creativecommons.org/licenses/by/4.0/", "doi_reg_from": "reg_local", "cstr_reg_from": "reg_local", "doi_not_reg_reason": null, "cstr_not_reg_reason": null, "is_paper_in_submitting": false, "belong_to_nieer": false, "ds_topic_tags": [ "输沙势", "合成输沙势", "合成输沙方向", "风向变率", "频数", "频率" ], "ds_subject_tags": [ "应用气象学", "自然地理学" ], "ds_class_tags": [], "ds_locus_tags": [ "全球" ], "ds_time_tags": [ 1940, 1941, 1942, 1943, 1944, 1945, 1946, 1947, 1948, 1949, 1950, 1951, 1952, 1953, 1954, 1955, 1956, 1957, 1958, 1959, 1960, 1961, 1962, 1963, 1964, 1965, 1966, 1967, 1968, 1969, 1970, 1971, 1972, 1973, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981, 1982, 1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019, 2020, 2021, 2022, 2023, 2024, 2025 ], "ds_contributors": [ { "true_name": "李国帅", "email": "liguoshuai@lzb.ac.cn", "work_for": "中国科学院西北生态环境与资源研究所", "country": "中国" }, { "true_name": "李新", "email": "xinli@itpcas.ac.cn", "work_for": "中国科学院青藏高原研究所", "country": "中国" } ], "ds_meta_authors": [ { "true_name": "李国帅", "email": "liguoshuai@lzb.ac.cn", "work_for": "中国科学院西北生态环境与资源研究所", "country": "中国" } ], "ds_managers": [ { "true_name": "李国帅", "email": "liguoshuai@lzb.ac.cn", "work_for": "中国科学院西北生态环境与资源研究所", "country": "中国" } ], "category": "沙漠与荒漠化" }