高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于高分辨率再分析风场的高原涡三维识别技术及应用

汤欢 傅慎明 孙建华 周象贤

汤欢, 傅慎明, 孙建华, 等. 2023. 基于高分辨率再分析风场的高原涡三维识别技术及应用[J]. 大气科学, 47(3): 698−712 doi: 10.3878/j.issn.1006-9895.2112.21127
引用本文: 汤欢, 傅慎明, 孙建华, 等. 2023. 基于高分辨率再分析风场的高原涡三维识别技术及应用[J]. 大气科学, 47(3): 698−712 doi: 10.3878/j.issn.1006-9895.2112.21127
TANG Huan, FU Shenming, SUN Jianhua, et al. 2023. Three-Dimensional Objective Identification of the Tibetan Plateau Vortex Based on a Reanalysis Wind Field with a High Spatial and Temporal Resolution [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 47(3): 698−712 doi: 10.3878/j.issn.1006-9895.2112.21127
Citation: TANG Huan, FU Shenming, SUN Jianhua, et al. 2023. Three-Dimensional Objective Identification of the Tibetan Plateau Vortex Based on a Reanalysis Wind Field with a High Spatial and Temporal Resolution [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 47(3): 698−712 doi: 10.3878/j.issn.1006-9895.2112.21127

基于高分辨率再分析风场的高原涡三维识别技术及应用

doi: 10.3878/j.issn.1006-9895.2112.21127
基金项目: 国网浙江省电力公司科技项目5211DS19001W,国家自然科学基金项目41775046、42075002,中国科学院青年创新促进会项目,国家重大科技基础设施项目“地球系统数值模拟装置”,高原与盆地暴雨旱涝灾害四川省重点实验室开放研究基金项目
详细信息
    作者简介:

    汤欢,女,1995 年出生,博士研究生,主要从事中尺度气象学研究。 E-mail: tanghuan@mail.iap.ac.cn

    通讯作者:

    傅慎明,E-mail: fusm@mail.iap.ac.cn

  • 中图分类号: P447

Three-Dimensional Objective Identification of the Tibetan Plateau Vortex Based on a Reanalysis Wind Field with a High Spatial and Temporal Resolution

Funds: State Grid Zhejiang Electric Power Company Science and Technology Project (Grant 5211DS19001W), National Natural Science Foundation of China (Grants 41775046, 42075002), Youth Innovation Promotion Association, Chinese Academy of Sciences, the National Key Scientific, Technological Infrastructure Project "Earth System Science Numerical Simulator Facility", Heavy Rain and Drought-Flood Disasters in Plateau and Basin Key Laboratory of Sichuan Province
  • 摘要: 高原涡(TPV)是生成于青藏高原主体的一类浅薄中尺度涡旋系统,其发生频繁、影响范围广、造成灾害强,是我国最重要的致灾中尺度系统之一。全面揭示高原涡的统计特征是本领域研究的重要基础。其中,高原涡的精准识别是认识其统计特征的关键。随着高时空分辨率再分析资料的出现,高原涡的研究有了更好的数据基础,然而,无论是人工识别方法还是基于较粗分辨率的客观识别算法都难以高效地适用于当前的新再分析资料。因此,亟需发展一种高精度的、适用于高时空分辨率再分析资料的高原涡客观识别方法。本文提出了一种适用于高分辨率再分析资料、基于风场的限制涡度高原涡客观识别算法(Restricted-vorticity based Tibetan-Plateau-vortex Identifying Algorithm,简称RTIA)。该方法首先判断高原涡候选点,然后以候选点为中心,划分多个象限,通过象限平均风场限定条件和象限组逆时针旋转(北半球)条件确定高原涡中心,无需复杂计算及对各气压层分别设定阈值,即可快速实现高原涡的水平和垂直追踪。基于1979~2020年共42个暖季(5~9月)、15466个高原涡(共计99090时次)大样本的评估表明,RTIA方法识别高原涡的平均命中率超过95%,平均空报率低于9%,平均漏报率少于5%,可以十分准确地对高原涡进行识别。此外,评估还表明RTIA方法应用于不同空间分辨率的再分析资料(如0.5°或0.25°)时,仍能保持高原涡识别的高准确率,其识别结果主要受涡旋自身强度的影响,对弱涡旋的识别精度比强涡旋偏低。该方法对其他中尺度涡旋识别也具有一定的借鉴意义。
  • 图  1  高原涡客观识别(RTIA)方法流程图

    Figure  1.  Flowchart of the RTIA (Restricted-vorticity-based Tibetan-Plateau-vortex Identifying Algorithm) method

    图  2  (a–d)2009年7月3日14:00(协调世界时,下同)和(e–h)2010年6月1日23:00的500 hPa流场(灰色箭头)分布:(a、e)涡度(阴影,单位:10−5 s−1);(b、f)限制涡度(阴影,单位:10−5 s−1);(c、g)有效涡度(阴影,单位:10−5 s−1);(d、h)位势高度场(阴影区,单位:gpm)。蓝色线表示高原边界;符号“×”表示高原涡可能的中心,(a、e)由涡度的极大值确定,(b、f)由限制涡度的极大值确定,(c、g)由有效涡度的极大值确定,(d、h)由位势高度的极小值确定

    Figure  2.  Distributions of 500 hPa streamline field (gray line with arrows) at 1400 UTC July 3, 2009 (left column) and 2300 UTC June 1, 2010 (right column): Panels (a, e) Relative vorticity (shading, units: 10−5 s−1); (b, f) restricted vorticity (shading, units: 10−5 s−1); (c, g) effective vorticity (shading, units: 10−5 s−1); (d, h) indicate the geopotential height (shading, units: gpm). Blue lines show the Tibetan Plateau boundaries; the black crosses show the candidate TPV centers, in (a, e) the centers are determined using the maxima of relative vorticity, in (b, f) the centers are determined using the maxima of restricted vorticity, in (c, g) the centers are determined using the maxima of effective vorticity, in (d, h) the centers are determined using the minima of geopotential height

    图  3  (a)四象限和(b)八象限限制涡度涡旋逆时针环流识别方案示意图[红色的点${O}'$表示候选点中心,识别半径Rf确定以${O}'$为起始点用于计算的风场数据,红色箭头表示各象限平均风场(u表示纬向风,v表示经向风)满足气旋性旋转,黑色虚线箭头表示各象限组(如第一象限和第二象限,第二象限和第三象限等)平均风向呈现逆时针旋转]

    Figure  3.  Schematic of the four-quadrant (a) and eight-octant (b) vortex counterclockwise identification approach, where the large red dot at the center is the center of the candidate vortex (${O}'$), and the identification radius (Rf) determined the wind data from the center of the candidate vortex (${O}'$) used for the calculation. The red arrow indicates that the octant-averaged wind should satisfy the cyclonic pattern, as shown by u (zonal wind) and v (meridional wind). The thick black dashed circle with arrows indicates the average wind direction of the octant group (e.g., the first and second octants; the second and third octants) is counterclockwise

    图  4  高原涡源地的识别范围(阴影表示地形,单位:m;黑色线表示高原边界)

    Figure  4.  Identification area (red rectangle) of the TPV source (the shaded area indicates terrain, units: m; the black solid line indicates the Tibetan Plateau’ s boundary)

    图  5  高原涡垂直追踪示意图。2018年5月14日06:00 450~550 hPa流场,阴影为相对涡度,单位:10−5 s−1,紫色线为高原边界线,R表示垂直搜索半径,$ {C}_{k} $表示该层(第k层)追踪的高原涡中心,$ {C}_{k+1} $表示向下一层追踪的高原涡中心,$ {C}_{k-1} $表示向上一层追踪的高原涡中心

    Figure  5.  Schematic of TPV Vertical Tracking. Streamlines field at 0600 UTC May 14, 2018, from 450 to 550 hPa; the shaded area indicates relative vorticity, units: 10−5 s−1; the purple solid line shows the boundary of the Tibetan Plateau; R is the vertical search radius; $ {C}_{k} $ represents the center of TPV tracked by the pressure layer (pressure layer k);$ {C}_{k+1} $ represents the center of TPV tracked by the next lower pressure layer;$ {C}_{k-1} $ represents the center of TPV tracked by the next upper-pressure layer

    图  6  2002年7月(a)26日15:00、(b)27日15:00、(c)29日10:00 500 hPa风场(黑色风羽,一根长羽为4 m s−1)、高度场(蓝色等值线,单位:gpm)、涡度场(阴影区,单位:10−5 s−1)及RTIA方法识别的高原涡中心(绿点)分布,褐色线为高原边界

    Figure  6.  Distributions of 500 hPa wind field (black wind barb, a full bar indicates 4 m$ {\mathrm{s}}^{-1} $), geopotential height field (blue contours, units: gpm), relative vorticity field (shading, units: 10−5 s−1) and the center of the TPVs (green dots) by RTIA at (a) 1500 UTC July 26, (b) 1500 UTC July 27, and (c) 1000 UTC July 29, 2002. The brown lines indicate the Tibetan Plateau’ s boundaries

    图  7  2002年7月27日15:00 200~700 hPa风场(黑色风羽,一根长羽为4 m s−1)、高度场(蓝色等值线,单位:gpm)、涡度场(阴影区,单位:10−5 s−1)及RTIA方法识别的高原涡中心(绿点)分布

    Figure  7.  Distributions of wind field (black wind barbs, a full bar indicates 4 m s−1), geopotential height field (blue contours, units: gpm), relative vorticity field (shading, units: 10−5 s−1), and the TPVs’ center (green dots) by RTIA with 50 hPa intervals from 200 to 700 hPa at 1500 UTC July 27, 2002

    图  8  2000年5月25日(a)21:00、(b)22:00、(c)23:00 500 hPa流场分布。阴影为相对涡度,单位:10−5 s−1;红点代表高原涡中心

    Figure  8.  Distributions of 500 hPa stream fields (shading: relative vorticity, units: 10−5 s−1; red dots: the centers of the TPVs) at (a) 2100 UTC, 2200 UTC, and (c) 2300 UTC May 25, 2000

    图  9  1979~2020年5~9月高原涡生命史分布(柱子上方数字表示生命史在该区间内的高原涡数量占全部高原涡的百分比,橙色表示RTIA方法,绿色表示实况)

    Figure  9.  The lifespan distribution of TPVs from May to September during 1979–2020, the number above the column represents the percentage of the number of TPVs within this lifespan period in the total TPVs from the RTIA (orange) and the actual situation (green)

    图  10  1979~2020年5~9月高原涡各时次2°×2°风场相关系数(绿线表示该年风场相关系数从小到大排第5百分位对应的值,灰色柱子是第10百分位对应的值,红线是该年所有时次风场相关系数的平均值)

    Figure  10.  Temporal correlation coefficient of 2°×2° wind field of TPVs from May to September during 1979–2020. The green line represents the value corresponding to the 5th percentile of wind field correlation coefficient from small to large in this year, the gray column represents the value corresponding to the 10th percentile, and the red line is the average value of all temporal wind correlation coefficients in this year

    图  11  1979~2020年5~9月高原涡RTIA效果评估(a)箱形图(线段最高点为最大值,线段最低点为最小值,箱型上部框线为第75%分位值,箱形下部框线为第25%分位值,箱内线为中位数,箱内点为平均值)以及(b)年变化(红线:命中率,蓝线:空报率,黑线:漏报率)

    Figure  11.  (a) Box-and-whisker plot (the lower and upper whiskers cover the minima and maxima value, the boxes cover the 25th–75th percentiles, the horizontal lines in the boxes mark the median value, and the dots represent the mean value) and (b) annual variations (the red line is the hit ratio, the blue line is the false alarm ratio, and the black line is the missing ratio) of TPVs from the RTIA from May to September during 1979–2020

    表  1  基于ERA5风场的不同涡旋中心识别的效果验证

    Table  1.   Validation of the effect of identifying different vortex centers based on the ERA5 wind field

    样本时段涡旋类型检出数命中率空报率漏报率
    1999年5~9月高原涡(500 hPa)309995.98%5.29%4.02%
    2008年5~6月西南涡(700 hPa) 41395.72%2.42%4.28%
    2015年7~8月大别山涡(850 hPa) 47996.06%3.34%3.94%
    下载: 导出CSV
  • [1] Blender R, Fraedrich K, Lunkeit F. 1997. Identification of cyclone-track regimes in the North Atlantic [J]. Quart. J. Roy. Meteor. Soc., 123(539): 727−741. doi: 10.1002/qj.49712353910
    [2] Curio J, Chen Y R, Schiemann R, et al. 2018. Comparison of a manual and an automated tracking method for Tibetan Plateau vortices [J]. Adv. Atmos. Sci., 35(8): 965−980. doi: 10.1007/s00376-018-7278-4
    [3] Curio J, Schiemannm R, Hodges K I, et al. 2019. Climatology of Tibetan Plateau vortices in reanalysis data and a high-resolution global climate model [J]. J. Climate, 32(6): 1933−1950. doi: 10.1175/JCLI-D-18-0021.1
    [4] Feng X Y, Liu C H, Rasmussen R, et al. 2014. A 10-yr climatology of Tibetan Plateau vortices with NCEP climate forecast system reanalysis [J]. J. Appl. Meteorol. Climatol., 53(1): 34−46. doi: 10.1175/JAMC-D-13-014.1
    [5] Fu S M, Li W L, Sun J H, et al. 2015. Universal evolution mechanisms and energy conversion characteristics of long-lived mesoscale vortices over the Sichuan Basin [J]. Atmos. Sci. Lett., 16(2): 127−134. doi: 10.1002/asl2.533
    [6] Fu S M, Zhang J P, Sun J H, et al. 2016. Composite analysis of long-lived mesoscale vortices over the middle reaches of the Yangtze River valley: Octant features and evolution mechanisms [J]. J. Climate, 29(2): 761−781. doi: 10.1175/JCLI-D-15-0175.1
    [7] Fu S M, Mai Z, Sun J H, et al. 2019. Impacts of convective activity over the Tibetan Plateau on plateau vortex, southwest vortex, and downstream precipitation [J]. J. Atmos. Sci., 76(12): 3803−3830. doi: 10.1175/JAS-D-18-0331.1
    [8] Fu S M, Zhang J P, Tang H, et al. 2020. A new mesoscale-vortex identification metric: Restricted vorticity and its application [J]. Environ. Res. Lett., 15(12): 124053. doi: 10.1088/1748-9326/abcac6
    [9] 关良, 李栋梁. 2019. 青藏高原低涡的客观识别及其活动特征 [J]. 高原气象, 38(1): 55−65. doi: 10.7522/j.issn.1000-0534.2018.00067

    Guan L, Li D L. 2019. Objective identifying and activity characteristics of Qinghai-Tibetan Plateau vortex [J]. Plateau Meteorology (in Chinese), 38(1): 55−65. doi: 10.7522/j.issn.1000-0534.2018.00067
    [10] Han Y Y J, Yang Q K, Liu N N, et al. 2021. Analysis of wind-speed profiles and optical turbulence above Gaomeigu and the Tibetan Plateau using ERA5 data [J]. Mon. Not. R. Astron. Soc., 501(4): 4692−4702. doi: 10.1093/mnras/staa2960
    [11] Hersbach H, Bell B, Berrisford P, et al. 2020. The ERA5 global reanalysis [J]. Quart. J. Roy. Meteor. Soc., 146(730): 1999−2049. doi: 10.1002/qj.3803
    [12] Holton J R. 2004. An Introduction to Dynamic Meteorology [M]. 4th ed. Amsterdam: Academic Press.
    [13] Hou J, Wang P, Zhuang S. 2017. A new method of characterizing flow patterns of vortices and detecting the centers of vortices in a numerical wind field [J]. J. Atmos. Oceanic Technol., 34(1): 101−115. doi: 10.1175/JTECH-D-15-0197.1
    [14] 黄一航, 陈权亮, 汪正林, 等. 2018. 1979-2016年青藏高原低涡的变化特征研究 [J]. 暴雨灾害, 37(2): 97−104. doi: 10.3969/j.issn.1004-9045.2018.02.001

    Huang Y H, Chen Q L, Wang Z L, et al. 2018. Variation characteristics of Tibetan Plateau vortex during 1979-2016 [J]. Torrential Rain and Disasters (in Chinese), 37(2): 97−104. doi: 10.3969/j.issn.1004-9045.2018.02.001
    [15] Huang J P, Chen W, Wen Z P, et al. 2019. Review of Chinese atmospheric science research over the past 70 years: Climate and climate change [J]. Sci. China Earth Sci., 62(10): 1514−1550. doi: 10.1007/s11430-019-9483-5
    [16] 李国平. 2021. 青藏高原动力气象学[M]. 3版. 北京: 气象出版社, 81–88.

    Li G P. 2021. Dynamic Meteorology of the Tibetan Plateau (in Chinese) [M]. 3rd ed. Beijing: China Meteorological Press, 81–88.
    [17] Li L, Zhang R H, Wu P L. 2020a. Evaluation of NCEP-FNL and ERA-interim data sets in detecting Tibetan Plateau vortices in May–August of 2000-2015 [J]. Earth Space Sci., 7(3): e2019EA000907. doi: 10.1029/2019EA000907
    [18] Li L, Zhang R H, Wu P L, et al. 2020b. Roles of Tibetan Plateau vortices in the heavy rainfall over southwestern China in early July 2018 [J]. Atmos. Res., 245: 105059. doi: 10.1016/j.atmosres.2020.105059
    [19] 林志强, 周振波, 假拉. 2013. 高原低涡客观识别方法及其初步应用 [J]. 高原气象, 32(6): 1580−1588. doi: 10.7522/j.issn.1000-0534.2012.00153

    Lin Z Q, Zhou Z B, Jia L. 2013. Objective identifying method of Qinghai–Xizang Plateau vortex using NCEP/NCAR reanalysis dataset [J]. Plateau Meteorology (in Chinese), 32(6): 1580−1588. doi: 10.7522/j.issn.1000-0534.2012.00153
    [20] Lin Z Q, Guo W D, Jia L, et al. 2020. Climatology of Tibetan Plateau vortices derived from multiple reanalysis datasets [J]. Climate Dyn., 55(7): 2237−2252. doi: 10.1007/s00382-020-05380-6
    [21] 卢敬华. 1986. 西南低涡概论[M]. 北京: 气象出版社, 270pp

    Lu J H. 1986. The introduction of the Southwest Vortices (in Chinese) [M]. Beijing: China Meteorological Press, 270pp.
    [22] Lu C H, Kong Y, Guan Z Y. 2020. A mask R-CNN model for reidentifying extratropical cyclones based on quasi-supervised thought [J]. Sci. Rep., 10(1): 15011. doi: 10.1038/s41598-020-71831-z
    [23] Orlanski I. 1975. A rational subdivision of scales for atmospheric processes [J]. Bull. Amer. Meteor. Soc., 56(5): 527−530. doi: 10.1175/1520-0477-56.5.527
    [24] 彭广, 李跃清, 郁淑华, 等. 2011. 青藏高原低涡切变线年鉴(2009) [M]. 北京: 科学出版社.

    Peng G, Li Y Q, Yu S H, et al. 2011. The Yearbook of Tibetan Plateau Vortex and Shear Line (2009) (in Chinese) [M]. Beijing: Science Press.
    [25] Tao S Y, Ding Y H. 1981. Observational evidence of the influence of the Qinghai–Xizang (Tibet) Plateau on the occurrence of heavy rain and severe convective storms in China [J]. Bull. Amer. Meteor. Soc., 62(1): 23−30. doi: 10.1175/1520-0477(1981)062<0023:OEOTIO>2.0.CO;2
    [26] Xie M, Li Y, Cao K. 2020. Global cyclone and anticyclone detection model based on remotely sensed wind field and deep learning [J]. Remote Sens., 12(19): 3111. doi: 10.3390/rs12193111
    [27] 叶笃正, 李麦村. 1964. 中小尺度运动中风场和气压场的适应 [J]. 气象学报, 34(4): 409−423. doi: 10.11676/qxxb1964.041

    Ye D Z, Li M C. 1964. The adaptation between the pressure and the wind field in the meso-and small-scale motion [J]. Acta Meteorologica Sinica (in Chinese), 34(4): 409−423. doi: 10.11676/qxxb1964.041
    [28] 叶笃正, 高由禧. 1979. 青藏高原气象学[M]. 北京: 科学出版社, 122–126

    Ye D Z, Gao Y X. 1979. The Tibetan plateau meteorology (in Chinese) [M]. Beijing: Science Press, 122–126.
    [29] 郁淑华, 高文良. 2017. 高原低涡与西南涡结伴而行的不同活动形式个例的环境场和位涡分析 [J]. 大气科学, 41(4): 831−856. doi: 10.3878/j.issn.1006-9895.1612.16213

    Yu S H, Gao W L. 2017. Analysis of environmental background and potential vorticity of different accompanied moving cases of Tibetan Plateau vortex and Southwest China vortex [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 41(4): 831−856. doi: 10.3878/j.issn.1006-9895.1612.16213
    [30] 张敬萍, 傅慎明, 孙建华, 等. 2015. 夏季长江流域两类中尺度涡旋的统计与合成研究 [J]. 气候与环境研究, 20(3): 319−336. doi: 10.3878/j.issn.1006-9585.2015.14164

    Zhang J P, Fu S M, Sun J H, et al. 2015. A statistical and compositional study on the two types of mesoscale vortices over the Yangtze River basin [J]. Climatic and Environmental Research (in Chinese), 20(3): 319−336. doi: 10.3878/j.issn.1006-9585.2015.14164
    [31] 张博, 李国平, 段炼, 等. 2018. 基于客观识别技术的高原低涡近30a气候特征 [J]. 兰州大学学报(自然科学版), 54(1): 104−111,119. doi: 10.13885/j.issn.0455-2059.2018.01.016

    Zhang B, Li G P, Duan L, et al. 2018. Climatic characteristics of Tibetan Plateau vortex based on the objective identification in the recent 30 years [J]. Journal of Lanzhou University (Natural Sciences) (in Chinese), 54(1): 104−111,119. doi: 10.13885/j.issn.0455-2059.2018.01.016
    [32] Zhao P, Xu X D, Chen F, et al. 2018. The third atmospheric scientific experiment for understanding the Earth–Atmosphere Coupled System over the Tibetan Plateau and its effects [J]. Bull. Amer. Meteor. Soc., 99(4): 757−776. doi: 10.1175/BAMS-D-16-0050.1
    [33] Zhao Y, Chen D L, Deng Y, et al. 2021. How were the eastward-moving heavy rainfall events from the Tibetan Plateau to the lower reaches of the Yangtze River enhanced? [J]. J. Climate, 34(2): 607−620. doi: 10.1175/JCLI-D-20-0226.1
  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  611
  • HTML全文浏览量:  107
  • PDF下载量:  148
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-07-19
  • 录用日期:  2021-12-07
  • 网络出版日期:  2022-01-05
  • 刊出日期:  2023-05-15

目录

    /

    返回文章
    返回