陆面模式公里小时大气驱动场及其扰动场的建立方案.pdfVIP

陆面模式 5 公里 3 小时大气驱动场及其扰动场的建立方案 大气驱动数据包括七个变量：近地面2米温度、地表压强、近地面2米相对 湿度、近地面 10 米风速、降水率、下行短波辐射和下行长波辐射。建立陆地模 式5公里3小时时间分辨率大气驱动数据的技术路线是：首先以薄板平滑样条函 数和已存在的格点驱动资料为自变量，以台站观测资料为应变量，建立驱动场的 趋势面；然后对趋势面的残差场进行订正；最后估计的订正场的误差协方差矩阵， 并根据它随机生成已建立的驱动场的扰动场集合。 1 温度场  1.1 观测资料的处理  1958-1989年的观测数据是日尺度的，每天台站观测个数大概在670个左右 (图 1)。在这个时间段里，用日尺度的观测资料建立模型得到日尺度的5公里驱 动场数据。然后用再分析资料进行3小时分配（见步骤3）。   图 1. 项目所用的中国大陆区域观测台站的位置及子区域的划分 在1990-2010年的温度观测数据中，1990-1997年是6小时分辨率的观测数 据，这时可利用6小时数据对3小时数据的插值估计出3小时的观测数据；虽然 1998-2010年的观测数据的时间分辨率是3小时，但是1998-2006年期间每天的 3、9、12、15 时的观测数据量比 0、6、12、18 四个时刻的观测数据量少很多， 所以我们首先要把数据量少的这些时刻的观测数据补全（见本节的附 1）。通过 以上方法，1990-2010年的3小时观测数据就准备好了。 1.2 建立方案  1.2.1 利用薄板平滑样条模型建立趋势面  t(x, y) f (x, y) =+α⋅z(x, y) +ε(x, y),                       (1)  y 其中 为经度， 为纬度，t(x, y) 为对应位置温度的真值，z (x ,y ) 为对应位 x 置的高程，f (x ,y ) 为薄板平滑样条函数， 为待估计的线性回归量系数，ε(x ,y ) α 2 为均值为 0，方差为σ （空间统计独立的）的误差变量。对于每一个时刻，利用 观测数据和经度纬度高程信息训练模型 Eq.(1)从而估计出f (x ,y ) 的各项系数和 2 α 值，以及ε(x ,y ) 的方差σ 。将5 公里格点的经纬度和高程信息输入模型 Eq.(1) 就得到趋势面，即ttrend (x, y) f (x, y) =+α⋅z(x, y) 。上述过程可以通过 R 软件“gam ” 程序包中的“mgcv ”和“predict ”两个函数实现。  为了加快计算速度，中国大陆区域被划分为四个边界重叠的子区域：W(西 部),Mw(中西部)，Ne （北东部），Se （南东部）（图  1）。先在每一个子区域上建 立趋势面，再将边界重叠的部分按到两边界距离加权平均，以形成整个中国大陆 区域的趋势面。  1.2.2 趋势面的订正  利用平滑样条方法建立的趋势面可能在有观测值的地方与观测的离差值远 远大于观测误差，所以我们用观测对趋势面进行订正。订正后的场就是我们估计 的驱动场的温度场。具体的订正方法如下：  tadj (x, y) ≡ttrend (x, y)   m( x, y ) m( x, y ) 2 + ∑w (x , y )×(t (x , y ) −t (x , y )) ∑w(x , y ),        (2)  i i obs i i trend i i i i