可降水量(PWV)对中小尺度灾害性天气预报具有指示意义,是降雨强度预报的重要参数之一。自20世纪80、90年代起,国外已经对GPS探测大气可降水量进行了大量研究。国内最早由毛节泰引进并详细介绍了地基GPS反演大气可降水量的方法,近年来发展迅速并取得了一定的成果。随着中国北斗卫星导航系统(BDS)提供服务并不断地建设完善,其已经能够应用到气象学领域,可以获得大气可降水量。文献[11]利用河北省6个测站的数据分别估算GPS、BDS天顶总延迟量(ZTD),两者之间的标准差为5~6mm。文献[12]利用亚太地区10个测站BDS数据采用精密单点定位方法进行大气可降水量反演,反演结果与GPS结果相比均方根误差为2mm左右。文献[13]采用精密单点定位方法,分别采用3种方案(BDS、GPS和BDS+GPS)来估算PWV,BDS精密单点定位可作为新的数据源,可独立提供高精度的ZTD与PWV产品且精度较高,同时BDS+GPS组合系统反演得到的可降水量精度更高。文献[14]利用上海市气象局建立的北斗气象站的观测数据及PANDA软件实现了基于北斗数据的大气可降水量反演,与探空数据计算结果及GPS数据反演结果之间的均方根误差约为3.5mm,相关系数均在0.95以上。
此外,GAMIT软件自10.5版本开始,就逐步加入对其他全球卫星导航系统(GNSS)的支持。随着10.61版本的发布,软件已经可以支持北斗观测数据的解算,并且附带的共用表文件也包含了多个GNSS系统的信息。截至目前,GAMIT软件已更新到10.70版本。因此,本文首次利用GAMIT软件对利用北斗系统数据反演大气可降水量进行研究,并结合探空数据计算得到的PWV与GPS反演得到的PWV进行精度验证,以期增强北斗卫星导航系统在大气水汽监测中的研究和应用。
1 大气可降水量计算方法1.1 无线电探空资料大气可降水量计算方法
基于探空气球的大气可降水量是利用探空气球观测气象对流层中各高度上的比湿q(单位g/kg),然后对大气压力p(单位hPa)从地面到对流层上界进行垂直积分而得,即
(1)
式中,PWV为可降水量,单位为mm;g为地球重力加速度。实际计算时,可以用各个标准等压面上的比湿值进行差分计算以代替积分,而且由于大气水汽几乎全部集中在气象学定义的对流层(0~12km),计算时采用分层叠加的办法进行。计算步骤如下:
(1) 根据各高度上的露点温度测值,算出各高度上的水汽压e(单位hPa),计算式为
(2)
式中,T