论南极气象学的现状及管理方案.doc

论南极气象学的现状及管理制度 　　摘要：水汽是一种在全球气候系统中起关键作用的温室气体，它不仅吸收和辐射太阳能量，对低层大气的化学作用也有影响，我们设法利用GPS定位中的误差源来反演大气中的水汽。 　　1地基GPS气象学技术的基本原理 　　如果地面上的GPS接收机的精确三维坐标通过其他方法确定，则可从GPS接收机测量的卫星信号传播的总延迟量中提取水汽造成的湿延迟，进而反演出信号路径上水汽的累积量。这就是地基GPS气象学的基本原理。GPS反演大气可降水量流程如图1所示。 　　2地基GPS反演大气可降水量处理流程 　　在应用GAMIT软件处理GPS数据反演可降水量时，GPS数据处理按照4个步骤进行。首先是数据准备，包括GPS观测数据(rinex格式的o文件)、广播星历(brdc文件)、精密星历文件(igs)、气象观测数据和相关表文件;接着用GAMIT软件处理各时段GPS的观测数据，得到一个基线解算文件;然后用GLOBK(卡尔曼滤波)进行多时段综合解算，以得到网平差的结果及各个测站的平均坐标和速度等参数;最后把精确的测站坐标带回GAMIT，对待求测站进行强约束，反演对流层的天顶延迟，得到可降水量。如果解算时发现各个时段的坐标相差较大，应该分析原因，重新进行解算。 　　3南极地基GPS数据解算实例 　　3.1站点选择和数据准备 　　本次实验所选取的站点采取引入长距离测站的方式获得绝对PWV。本次实验参与解算的站包括7个IGS站，分别是VESL，SYOG，MAW1，DAV1，CAS1，DUM1，MCM4;此外，还有中国的两个南极跟踪站长城站GRW1和中山站ZHN1(如图3所示)。为了在解算时较精确地得到可降水量，我们从SOPAC网站上查到2009年7个IGS站的精确坐标，长城站和中山站的坐标可以通过解算得到。这样在解算时对坐标点给予较强的约束。本次解算的是2009年全年9个站的GPS数据，准备好所有的数据文件后，便可以执行GAMIT的批处理命令让其自动解算一年的数据，就可以自动解算出9个站一年的数据。 　　3.2可降水量的提取及其精度分析 　　3.2.1提取可降水量 　　在GAMIT中提取可降水量的命令为sh_metutil。由于参与计算的天数有365d，要比对的站有5个，逐天逐条每个站下命令将很麻烦。这时，可以利用Linux系统的SHELL的功能，编写一个小程序，实现三重循环。用VI编辑器编写一个三重循环，以CAS1站为例，编写一个名为xunhuan.sh的文件，然后通过命令chmod使之成为可执行文件，具体命令为:Chmod+xxunhuan.sh。通过调用./xunhuan.sh来批量执行几百条命令，这样就比较方便，其他几个站以此类推。得到各个站点每天的可降水量，然后通过grep命令提取各个站点的可降水量。通过数据的整理分析，每个站点整理出一个可降水量文件。 　　3.2.2GPS/PWV与探空Radio/PWV的比较 　　目前，世界气象部门对大气中水汽的探测主要还是依赖于探空气球的探测，为了验证GPS/PWV的正确性和稳定性，将GPS/PWV和Radio/PWV进行比较。本次参与解算的9个站中有5个附近同时安装有探空站，分别是SYOG，DAV1，CAS1，DUM1，MCM4。下面将就这5个站解算出的GPS/PWV和探空PWV进行比较，其中探空PWV数据可以从美国怀俄明州大学的网站上下载。把2009年5个站点的GPS/PWV和Radio/PWV进行比较。由于两者的PWV在数值上比较接近，数据量又较大，为了能在图上清晰的显示，图4至图8画出GPS/PWV的时间序列图，以及GPS/PWV减去Radio/PWV的差值的时间序列图。5幅图中上面的曲线代表GPS/PWV，下面的曲线代表它和Radio/PWV的差值。从定性上去分析图4至图8可以得出以下特点:1)除个别点外，GPS/PWV比探空Radio/PWV系统性偏大，这说明用GPS探测水汽有偏高的迹象，并且当PWV值偏大时，GPS所测的数值会表现的更敏感一些，也就是说当绝对PWV较小时两者更为接近。2)图像显示PWV两端较大。差不多是年积日的0～50和300～365这两个范围PWV较大，这个时间段恰好是南极的夏天，天气较为湿润，出现降雨的可能性大些，这样和实际降水测量较为一致，解算的结果正确地反映了降水量的变化趋势，体现了解算的合理性和正确性。个别站点如CAS1站和SYOG中间也有一段时间(150～200d)PWV也偏大，这与南极部分地区季风性气候有关，造成这一时段降水天气的可能性增大。由上述5个图表显示及分析可知，GPS/PWV和Ra-dio/PWV在整体上是比较一致的，少部分天数两者的差距较大。总的来说，5个站点的PWV值常年都比较小，有少部分天数PWV急剧上升，这些天数发生降雨的可能