化值按给定的目标函数度量方式达到最佳拟合，在各参数值的可行解范围内寻找各时段内各参数值的最优解；
[0054]S14、用所述陆面过程模式参数优化值时间序列代替原有的陆面过程模式参数，形成改进的陆面过程模式参数化方案。
[0055]步骤S14进一步包括:
[0056]S141、重复S11、S12、S13直至感兴趣的历史范围(如过去的某一年)结束，得到陆面过程模式参数化方案(例如土壤水参数化方案)中的各参数(包含饱和水力传导率、饱和土壤湿度、饱和土壤水势、Clapp和Hornberger经验常数等)优化值时间序列。
[0057]S142、利用前述步骤得到的陆面过程模式参数化方案(如土壤水参数化方案)中的各参数优化值时间序列改进参数化方案；
[0058]S15、使用所述改进参数化方案的陆面过程模式和拟研究时段的所述陆面过程输入数据，直接计算得出所述拟研究时段的陆面过程模式输出数据。
[0059]图3表示同化卫星遥感数据优化陆面过程模式的实施例，包含以下步骤:
[0060]S21、运行陆面过程模式，计算陆面变量的模拟值，作为陆面过程模式输出数据。
[0061]步骤S21进一步包括:
[0062]S211、根据研究问题的时空分辨率的要求，准备一定时空分辨率的大气强迫数据，包括近地面气温、近地面气压、近地面空气比湿、近地面全风速、地面向下短波辐射、地面向下长波辐射、地面降水率等；
[0063]S212、准备一定时空分辨率的地表参数，包括地表覆盖类型及其所占的比例、土壤质地(砂土和粘土)的比例、叶面积指数等，如模式自带有地表数据集，可以直接使用模式自带的地表数据集，则不需本步骤；
[0064]S213、陆面过程模式的起转过程(Spin-up):利用拟研究时段以前的长时间序列的大气强迫数据和地表参数运行陆面过程模式，以获得模式平衡状态；
[0065]S214、以模式的平衡状态为初始条件，利用拟研究时段的大气强迫数据和地表参数驱动陆面过程模式，计算陆面变量的模拟值。
[0066]S22、用数据同化算法校正所述陆面过程输出数据，得到输出数据同化值。
[0067]步骤S22进一步包括:
[0068]S221、用陆面过程模式输出数据(如表层土壤水分、地表温度、雪层深度和冠层温度以及植被类型、植被比例、土壤类型等参数)作为地表微波辐射传输模型的输入，运行地表微波辐射传输模型，生成模拟的卫星遥感数据；
[0069]S222、在有卫星遥感观测的时间步，用数据同化算法(如变分法、集合Kalman滤波、扩展Kalman滤波、粒子滤波等算法)同化卫星遥感观测(如地表微波辐射亮度温度值)，基于模拟的卫星遥感数据和观测的卫星遥感数据之间的偏差来校正陆面过程模式输出数据，得到拟同化的陆面变量的同化值；
[0070]S23、构建目标函数，将所述输出数据同化值引入陆面过程模式，通过数值优化方法计算，改变陆面过程模式参数值至所述目标函数最小。
[0071 ] 步骤S23进一步包括:
[0072]S231、构建目标函数(如均方根误差函数等)，以描述陆面状态变量的模拟值与同化值的拟合程度；
[0073]S232、对历史范围内(如过去的某一年)多个时段(如该年的各月份)，分别给定陆面过程模式某参数化方案(如土壤水参数化方案)中的各参数(如饱和水力传导率、饱和土壤湿度、饱和土壤水势、Clapp和Hornberger经验常数等)的可行解范围，用步骤S21的方法运行陆面过程模式，利用优化算法(如复合形混合演化算法SCE-UA等)极小化目标函数，通过反复迭代，使模式模拟值与同化值按给定的目标函数度量方式达到最佳拟合，在各参数值的可行解范围内寻找各时段内各参数值的最优解；
[0074]S24、用所述陆面过程模式参数优化值时间序列代替原有的陆面过程模式参数，形成改进的陆面过程模式参数化方案。
[0075]步骤S24进一步包括:
[0076]S241、重复S21、S22、S23直至感兴趣的历史范围(如过去的某一年)结束，可得到陆面过程模式某参数化方案(如土壤水参数化方案)中的各参数(如饱和水力传导率、饱和土壤湿度、饱和土壤水势、Clapp和Hornberger经验常数等)优化值时间序列。
[0077]S242、利用前述步骤得到的陆面过程模式参数化方案(如土壤水参数化方案)中的各参数优化值时间序列改进参数化方案；
[0078]S25、使用所述改进参数化方案的陆面过程模式和拟研究时段的所述陆面过程输入数据，直接计算得出拟研究时段所述陆面过程模式输出数据。
[0079]以上方法中所使用的卫星遥感观测为低频地表微波辐射亮度温度数据。所有能够提供此类数据的卫星传感器，例如搭载在NASA AQUA卫星上的高级微波扫描辐射计AMSR-E、美国雨云卫星系列Nimbus-5、6搭载的电子扫描微波辐射计ESMR、Nimbus-7的多通道微波辐射计SMMR、美国防卫气象卫星计划(DMSP)的特种微波成像仪SSM/1、美国和日本联合的热带降雨测量任务TRMM搭载的微波成像仪TM1、中国第二代极轨气象卫星风云3号(FY-3)、土壤湿度和海洋盐度卫星SM0S、土壤湿度主动被动探测卫星SMAP等，所提供的低频地表微波辐射亮度温度数据都可以用于该方法来优化陆面过程模式。本发明以优化非饱和土壤水模型为例，通过发展基于非饱和土壤水模型和扩展Kalman滤波算法并结合可变下渗能力模型VIC的土壤湿度同化方案，以月为同化窗，利用优化算法——复合形混合演化算法(SCE-UA)极小化目标函数，使模拟与同化的土壤湿度按给定的目标函数度量方式达到最佳拟合，得到同化时段(即上文所述感兴趣的历史范围)1986年非饱和土壤水模型各参数(饱和水力传导率、饱和土壤湿度、饱和土壤水势、Clapp和Hornberger经验常数)的优化值序列，然后利用优化的各参数优化值时间序列改进非饱和土壤水模型，最后利用改进的模型进行1986-1993年的数值模拟试验，结果表明:优化非饱和土壤水模型中的参数值，使得非饱和土壤水模型在物理机制上更为完善，提高了模型自身的模拟能力。本例为利用陆面数据同化理论及方法改进和完善陆面过程模式参数化方案研究提供了范例，具有较好的借鉴作用。
[0080]以上实施方式仅用于解释本发明的方法，而并非对本发明的限制，有关技术领域的普通技术人员在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变形，比如使用不同的陆面过程模式、不同的地表微波辐射传输模型、不同的同化算法以及不同的优化方法发展同化方案，同化不同的站点观测和卫星遥感观测。因此所有等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。
【主权项】
1.一种优化陆面过程模式的方法，其特征在于，包括以下步骤: 51、利用陆面过程模式输入数据，运行陆面过程模式，计算陆面变量的模拟值，生成陆面过程模式输出数据； 所述陆面过程模式输入数据包含大气强迫数据和地表参数； 52、用数据同化算法校正所述陆面过程模式输出数据，得到输出数据同化值； 53、构建目标函数，将所述输出数据同化值引入所述陆面过程模式，通过数值优化方法改变陆面过程模式参数值使所述目标函数值最小； 所述目标函数值最小时对应的陆面过程模式参数，作为陆面过程模式参数优化值； 54、对历史范围内多个时段重复步骤SI?S3，得到陆面过程模式参数优化值时间序列；用所述陆面过程模式参数优化值时间序列代替原有的陆面过程模式参数，形成优化参数化方案的陆面过程模式。
2.如权利要求1所述优化陆面过程模式的方法，其特征在于，还包含以下步骤 使用所述优化参数化方案的陆面过程模式和拟研究时段的所述陆面过程模式输入数据，直接计算得出所述拟研究时段的陆面过程模式输出数据。
3.如权利要求1或2所述优化陆面过程模式的方法，其特征在于， 所述大气强迫数据包含近地面气温、近地面气压、近地面空气比湿、近地面全风速、地面向下短波辐射、地面向下长波辐射、地面降水率。
4.如权利要求1或2所述优化陆面过程模式的方法，其特征在于， 所述地表参数包含地表覆盖类型及其所占的比例、土壤质地比例、叶面积指数；所述土壤质地指砂土或粘土。
5.如权利要求1或2所述优化陆面过程模式的方法，其特征在于， 所述步骤Si进一步包含:用拟研究时段以前的所述陆面过程模式输入数据运行所述陆面过程模式获得模式平衡状态； 以所述模式平衡状态为初始条件，用拟研究时段的所述陆面过程模式输入数据，运行所述陆面过程模式计算所述陆面过程模式输出数据。
6.如权利要求1或2所述优化陆面过程模式的方法，其特征在于， 所述数据同化算法包含变分法、集合Kalman滤波算法、扩展Kalman滤波算法、粒子滤波算法中至少一种。
7.如权利要求1或2所述优化陆面过程模式的方法，其特征在于， 所述步骤S2中，所述数据同化算法使用陆面变量模拟数据和陆面变量观测数据的偏差来校正陆面过程模式输出数据。
8.如权利要求1或2所述优化陆面过程模式的方法，其特征在于，所述步骤S2中， 所述数据同化算法使用模拟的卫星遥感数据与观测的卫星遥感数据的偏差来校正陆面过程模式输出数据； 所述模拟的卫星遥感数据是运用地表微波辐射传输模型对所述陆面过程模式输出数据进行计算得到。
9.如权利要求1或2所述优化陆面过程模式的方法，其特征在于，所述步骤S3中，所述目标函数为均方根误差函数。
10.如权利要求1或2所述优化陆面过程模式的方法，其特征在于，所述步骤S3中，所 述数值优化方法是复合形混合演化算法。
【专利摘要】本发明涉及一种优化陆面过程模式的方法，解决传统数据同化算法耗时长的技术问题，发明的方法包含：S1、运行陆面过程模式，计算陆面过程各物理量的模拟值；S2、用数据同化方法校正拟同化的物理量，得到陆面变量的同化值；S3、构建目标函数，将同化值引入陆面过程模式，通过数值优化方法，改变陆面过程模式参数值至目标函数最小，重复S1～S3得到模式参数的最优值时间序列；S4、用模式参数的最优值时间序列代替原模式参数，改进陆面过程模式。方案从物理机制上改进陆面过程模式参数化方案，提高了模式自身的模拟精度；直接运行改进的陆面过程模式得到高精度的模式输出数据，避免逐年逐日低效率同化过程。
【IPC分类】G06F19-00
【公开号】CN104573378
【申请号】CN201510033103
【发明人】张生雷, 陈良富, 苏林, 贾立
【申请人】中国科学院遥感与数字地球研究所
【公开日】2015年4月29日
【申请日】2015年1月22日