海南GPS网探测对流层水汽廓线的试验研究
Tropospheric Water Vapor Profiles Using GPS Network in Hainan
-
摘要: 介绍了2005年11月在海南地区进行的一次地基GPS小网观测试验。试验目的是研究利用组网的GPS倾斜路径观测进行对流层水汽层析反演, 给出站点上空水汽的垂直廓线结构信息。概述了试验中GPS原始数据处理方法以及层析反演的方法。将GPS层析得到的水汽垂直廓线与海口站探空观测的水汽廓线进行了对比, 结果表明:二者一致性较好, 均方根误差在0.5g·m-3左右, 层析结果较好地反映出试验期间水汽减少、大气变干的过程。另外, 采用3种不同的先验信息方案测试分析了GPS层析的结果, 表明GPS观测量对水汽先验信息有明显的调整作用。并对GPS水汽层析中可能存在的问题进行了讨论。试验证明高时间分辨率的GPS观测有能力层析出GPS测站上空水汽的廓线信息。Abstract: Water vapor plays a key role in atmospheric processes. Water vapor field is also one of the initial conditions needed by numerical weather prediction. Its distribution remains difficult to quantify due to water vapor's high variability in time and space and the deficiency of available measurements. The GPS has proved its capacity to measure the integ rated water vapor at zenith with the same accuracy as other methods, such as radiosonde and water vapor radiometers, which has been demonstrated by many experiments and GPS precipitable water vapor has been operated in many areas. Recent studies also show that it is possible to quantify the integ rated water vapor in the direction of signals of the GPS satellite. These observations can be used to study the water vapor tomography using GPS network. That is to use the water vapor amount along slant path of GPS satellites to form observation equations and to obtain the vertical profile information over GPS sites by solving these equations.The principle of tomography is introduced in which slant path observations from ground-based GPS is used and the method of calculating the observation equations is presented. During November in 2005, a small GPS network experiment is carried out in Hainan. The purpose of this experiment is to study the tomography tropospheric water vapor profiles based on GPS slant path observations. Without a priori information of water vapor, the absolute value of water vapor is difficult to resolve. After the method of GPS raw data processing is summarized, how to use a priori information is also analyzed. The period of tomography is from Nov 17 to 21 in 2005. Three different vertical tomographic resolutions of 600 m, 800 m and 1000 m are tested and analyzed. Results show that these resolutions have no large difference. Therefore, tomographic resolution from 500 m to 1000 m can be applied in practical use. In the experiment, tomographic vertical resolution of 1000 m is selected. Three different plans using a priori information are described.In the first plan, priori information is the earliest radiosonde at 00:00 on Nov 17, 2005. In the second plan, the latest radiosonde before the retrieved time is used as priori information. The average radiosondes are used in the third plan during the experiment. In the first and the third plan, priori information is hold but in the second plan priori information is changed and regenerated. Water vapor vertical profiles derived from GPS tomography are compared with that measured by radiosonde in Haikou site. The results show that they agree well with each other and the root-mean-square error (RMSE) is about 0.5 g·m -3. The process of water vapor decreasing and atmosphere drying is well reflected by GPS tomography. The tomographic results from different priori information used in tomographic method are also analyzed. It is demonstrated that water vapor priori information is obviously adjusted by GPS observation. And some problems in tomography are also discussed. Vertical profile distribution of tropospheric water vapor can be retrieved from GPS measurements with high time resolution through tomography inversion.
-
Key words:
- GPS;
- water vapor;
- tomography
-
图 4 2005年11月21日 (年积日325) 采用不同先验信息的层析结果与实测探空对比
(十字形:方案一, 即采用第321天00:00探空为先验值的层析结果; 圆形:方案二, 即采用时间最近的上次探空为先验值; 正方形:方案三, 即用平均探空为先验值; 菱形:第325天的反演时刻的单次探空)
Fig. 4 Tomographic results from different a priori information versus observations by radiosonde
(line with cross-shaped symbol is the first plan using radiosonde in 321 day of year as a priori values; line with circle symbol represents the second plan using the latest radiosonde before this time as a priori information; line with square symbol is the third plan using a priori information from average radionsonde; line with diamond symbol shows the single radiosonde in day 325)
图 5 试验期间3种先验值层析方案与探空对比的散点图
(a) 采用2005年11月17日00:00探空为先验值的层析结果, (b) 采用反演时刻前一次探空为先验值的层析结果, (c) 以平均探空廓线为先验值的结果
Fig. 5 Scatterplots of results from radiosonde and tomography during the experiment
(a) tomographic results using radiosonde at 00:00 as a priori information on Nov 17, 2005, (b) using latest radiosonde before the time of retrieval, (c) using average radiosonde
图 6 2005年11月21日GPS层析、先验信息分别与探空的对比散点图
(菱形:先验信息与探空的对比; 实线:通过原点的拟合直线; 三角形: GPS层析与探空的对比; 虚线:通过原点的拟合直线)
Fig. 6 Scatterplots of GPS tomography, a priori information and radiosonde on Nov 21, 2005
(diamond symbol denotes comparison between a priori information and radiosonde, solid line is fitting line passing origin; triangle symbol denotes the results from GPS and radiosonde; dashed line is fitting line passing origin)
表 1 不同分辨率下的层析结果统计分析
Table 1 Tomographic statistics results of different resolutions
-
[1] Bevis M, Businger S, Herring T A, et al. GPS meteorology; Remote sensing of atmospheric water vapor using the global positioning system. J Geophys Res, 1992, 97: 15787-15801. doi: 10.1029/92JD01517 [2] Rocken C, Ware R, Van Hove T, et al. Sensing atmospheric water vapor with the Global Positioning System. Geophys Res Lett, 1993, 20: 2631-2634. doi: 10.1029/93GL02935 [3] Ware R, Alber C, Rccken C, et al. Sensing integrated water vapor along GPS ray paths. Geophys Res Lett, 1997, 24:417-420. doi: 10.1029/97GL00080 [4] Alber C, Ware R, Rocken C, et al. Obtaining single path phase delays from GPS double differences. Geophys Res Lett, 2000, 27: 2661-2664. doi: 10.1029/2000GL011525 [5] Braun J, Rocken C, Ware R. Validation of line-of-sight water vapor measurements with GPS. Radio Sci, 2001, 36: 459-472. doi: 10.1029/2000RS002353 [6] Braun J, Rocken C, Liljegren James. Comparisons of line-of-sight water vapor observations using the Global Positioning System and a pointing microwave radiometer. J Atmos Ocean Technol, 2003, 20: 606-612. doi: 10.1175/1520-0426(2003)20<606:COLOSW>2.0.CO;2 [7] Flores A, Rius A, Ruffini G. 4D tropospheric tomography using GPS slant wet delays. Ann Geopkys, 2000, 18: 223-224. doi: 10.1007/s00585-000-0223-7 [8] Braun J, Rocken C. Water Vapor Tomography within the Planetary Boundary Layer Using GPS. International Workshop on GPS Meteorology, Tsukuba, Japan, 2003. https://www.researchgate.net/profile/John_Braun2/publication/228573377_Water_vapor_tomography_within_the_planetary_boundary_layer_using_GPS/links/00b7d5385ff8537f6e000000.pdf?origin=publication_detail [9] 毕研盟, 毛节泰, 刘晓阳, 等.应用地基GPS遥感倾斜路径方向大气水汽总量.地球物理学报, 2006, 49(2):335-342. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200602004.htm [10] 毛辉, 毛节泰, 毕研盟, 等.遥感GPS倾斜路径信号构筑水汽时空分布图.中国科学 (D辑), 2006, 36(12):1177-1186. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200612012.htm [11] Eugenia Kalnay.大气模式、资料同化和可预报性.蒲朝霞, 杨福全, 邓北胜, 等译.北京:气象出版社, 2005. [12] Ruffini G, Flores A, Rius A. GPS tomography of the ionospheic electron content with a correlation functional. IEEE Tram Geos Remote sensing, 1998, 36: 143-153. doi: 10.1109/36.655324 [13] Bi Yanmeng, Mao Jietai, Li Chengcai. Preliminary results of 4D water vapor tomography in the troposphere using GPS. Adv Atmos Sci, 2006, 23(4): 551-560. doi: 10.1007/s00376-006-0551-y [14] Hirahara K. Local GPS tropospheric tomography. Earth Planets Space, 2000, 52: 935-939. doi: 10.1186/BF03352308 [15] Champolliona C, Massona F, Bouin M, et al. GPS water vapour tomography: Preliminary results from the ESCOMPTE field experiment. Atmos Res, 2005, 74: 253-274. doi: 10.1016/j.atmosres.2004.04.003 [16] Bastin S, Champollion C, Bock O, et al. On the use of GPS tomography to investigate water vapor variability during a Mistral/sea breeze event in southeastern France. Geophy Res Lett, 2005, 32 : L05808, doi: 10. 1029/2004GL021907. [17] Noguchi W, Yoshihara T, Tsuda T, et al. Time-height distribution of water vapor derived by moving cell tomography during Tsukuba GPS campaigns. J Meteor Soc Japan, 2004, 82 (1B) : 561-568. doi: 10.2151/jmsj.2004.561