一种考虑逐日数据波动特征的GCM校正方法技术

本发明专利技术公开了一种考虑逐日数据波动特征的GCM校正方法，以月为计算单元，提取某一种气象因子的历史时期逐日实测数据和同历史时期GCM逐日数据以及待矫正的未来时期GCM逐日数据，计算历史实测数据和历史GCM数据多年月均值，构建均值矫正因子，从而对未来气象数据进行逐日均值矫正，并构建方差矫正因子，对均值矫正后的未来气象数据进行方差矫正，以矫正后的方差为约束条件，得到一个新的未来GCM气象序列。本发明专利技术考虑了未来GCM数据的分布特征和实测数据的分布特征一致性，对未来GCM数据进行偏差矫正，从而使得采用GCM数据的模型模拟的结果更趋于准确。

【技术实现步骤摘要】
一种考虑逐日数据波动特征的GCM校正方法
本专利技术涉及一种考虑逐日数据波动特征的GCM校正方法，属于气候分析统计
技术背景全球气候模式是评估气候变化对人类社会及自然系统影响的重要数据基础。然而，大气环流模式(GCM)普遍分辨率较低，且存在着区域系统性误差。因此，需根据气象站点实测数据对气候模式数据进行降尺度，以实现大气环流模式数据与气象站点数据的尺度匹配。然而，传统的降尺度或偏差校正方法仅从月均值的角度来进行气候模式数据的校正处理，没有对逐日数据的准确性进行分析判断。事实上，在很多研究领域，如农业生长耗水模拟、水文循环模拟领域，逐日数据的准确性将直接关系到研究结果的有效性，因此需对气候模式数据进行逐日校正。偏差校正的基本思想是通过某种转换函数使得GCM历史时期的气象要素尽可能地接近于实测气象序列，然后将这种转换函数应用至未来气候模式数据，以获得更加可靠的未来气象序列。
技术实现思路
本专利技术所要解决的技术问题是克服现有技术的缺陷，提供一种考虑逐日数据波动特征的GCM校正方法，以月为计算单元，以多年月均值和多年月均方差为约束条件，重新构造气象序列，使得重新构造后的气象序列不仅在均值上得到校正，同时使得逐日数据在更加合理的范围内波动。为解决上述技术问题，本专利技术提供一种考虑逐日数据波动特征的GCM校正方法，包括以下步骤：1)提取某一种气象因子的历史时期逐日实测数据和同历史时期GCM逐日数据以及待矫正的未来时期GCM逐日数据；2)计算所述气象因子历史逐日实测数据和历史GCM逐日数据各月均值，并计算二者多年月均值；3)基于历史逐日实测数据多年月均值与历史GCM逐日数据多年月均值，构建均值校正因子；4)待矫正的未来GCM逐日数据分别乘以所述均值校正因子得到均值校正后的气象序列；5)基于所述步骤1)提取的气象因子的原始数据，分别计算该气象因子历史逐日实测数据和历史GCM逐日数据各月方差，并计算多年月方差均值；6)构建方差校正因子；7)构建逐日气象因子转换函数。前述的步骤2)中，气象因子历史逐日实测数据多年月均值的计算式如下：其中，Mo,i为历史逐日实测数据的第i月多年均值，i＝1,2,3…12；Xo,i,k,j历史时期第i月第k年第j天实测值，k＝1,2,3…s，j＝1,2,3…,n，n＝28/29/30/31，具体n的值以该月的实际天数为准，s为年数；所述历史GCM逐日数据多年月均值的计算式如下：其中，Mm,i为历史GCM逐日数据的第i月多年均值，i＝1,2,3…12；Xm,i,k,j历史时期第i月第k年第j天GCM值，k＝1,2,3…s；j＝1,2,3…,n，n＝28/29/30/31，具体n的值以该月的实际天数为准，s为年数。前述的步骤3)中，均值校正因子为：其中，θi为第i月均值校正因子。前述的步骤5)中，气象因子历史逐日实测数据多年月方差均值的计算式为：其中，σo,i为历史逐日实测数据的第i月的多年方差均值，σo,k,j为历史实测第k年第i月方差；所述历史GCM逐日数据的多年月方差均值计算式为：其中，σm,i为历史GCM逐日数据的第i月的多年方差均值，σm,k,i为历史GCM第k年第i月方差。前述的步骤6)中，方差校正因子为：其中，βi为第i月方差校正因子。前述的逐日气象因子转换函数：Xd,i,j＝Xi,j-(Xi,j-Mf,i)×(1-βi)其中，Xi,j为均值校正后的第i月第j天的未来GCM逐日数据，j＝1,2,3…,n，n＝28/29/30/31，具体n的值以该月的实际天数为准，Mf,i是均值矫正后的未来GCM逐日数据气象序列的逐月均值，Xd,i,j为方差校正后的第i月第j天的未来GCM逐日数据。有益效果：本专利技术基于历史实测数据和历史GCM数据，从气象因子分布特征出发，构建均值矫正因子和方差矫正因子矫正未来GCM数据，矫正后的未来的GCM数据序列的特征，包括均值和波动程度特征更接近于实测值，从而使得基于未来GCM数据的模型模拟的结果更加准确。附图说明图1为本专利技术的考虑逐日数据波动特征的GCM校正方法流程图；图2为历史实测辐射及偏差矫正前GCM辐射图；图3为校正后GCM辐射数据历史实测数据。具体实施方式下面对本专利技术作进一步描述。以下实施例仅用于更加清楚地说明本专利技术的技术方案，而不能以此来限制本专利技术的保护范围。如图1所示，本专利技术的考虑逐日数据波动特征的GCM校正方法，包括以下步骤：(1)以月为计算单元，提取某一种气象因子的历史时期逐日实测数据和同历史时期GCM逐日数据以及待矫正的未来时期GCM逐日数据，待矫正的未来时期GCM逐日数据由某些科研机构提供。基于此数据，分别采用算术平均法计算该气象因子历史逐日实测数据和历史GCM逐日数据各月均值，并计算二者多年月均值，计算方法如下：其中，Mo,i为历史逐日实测数据的第i月多年均值，i＝1,2,3…12；Xo,i,k,j历史时期第i月第k年第j天实测值，k＝1,2,3…s，j＝1,2,3…,n，n＝28/29/30/31，具体n的值以该月的实际天数为准；s为年数，一般取年限为30年以上，精确度高。历史GCM逐日数据的多年月均值(Mm,i)计算方法与历史逐日实测数据的计算方法相同，可表示为：其中，Mm,i为历史GCM逐日数据的第i月多年均值，i＝1,2,3…12；Xm,i,k,j历史时期第i月第k年第j天GCM值，k＝1,2,3…s；j＝1,2,3…,n，n＝28/29/30/31，具体n的值以该月的实际天数为准。(2)基于历史逐日实测数据多年月均值与历史GCM逐日数据多年月均值，构建均值校正因子：其中，θi为第i月均值校正因子。待矫正的未来GCM逐日数据分别乘以该均值校正因子得到均值校正后的气象序列。(3)基于步骤(1)中提取的某一种气象因子的原始数据，分别计算该气象因子历史逐日实测数据和历史GCM逐日数据各月方差，并计算多年月方差均值：其中，σo,i为历史逐日实测数据的第i月的多年方差均值，σo,k,i为历史实测第k年第i月方差。历史GCM逐日数据的多年月方差均值(σm,i)计算方法与历史逐日实测数据计算方法相同，计算式为：其中，σm,i为历史GCM逐日数据的第i月的多年方差均值，σm,k,i为历史GCM第k年第i月方差。(4)基于该气象因子实测数据的多年月方差均值与GCM逐日数据的多年月方差均值，构建方差校正因子：其中，βi为第i月方差校正因子。(5)计算经均值校正后的未来GCM逐日数据气象序列的逐月均值(Mf,i)，构建逐日气象因子转换函数：Xd,i,j＝Xi,j-(Xi,j-Mf,i)×(1-βi)其中，Xi,j为均值校正后的第i月第j天的未来GCM逐日数据，Mf,i是均值矫正后的未来GCM逐日数据气象序列的逐月均值，Xd,i,j为方差校正后的第i月第j天的未来GCM逐日数据。以矫正后的方差为约束条件，得到一个新的未来GCM气候序列。实施例对BCC-CSM1.1(m)2030年昆山地区辐射数据进行偏差校正。具体步骤如下：(1)提取昆山地区1961-2000年逐日实测和GCM辐射数据以及待进行偏差校正的2030年逐日辐射数据。基于此数据，分别计算历史实测和GCM历史时期辐射数据各月均值，并采用算术平均方法计算实测辐射数据多年月均值，计算方本文档来自技高网...

【技术保护点】
一种考虑逐日数据波动特征的GCM校正方法，其特征在于，包括以下步骤：1)提取某一种气象因子的历史时期逐日实测数据和同历史时期GCM逐日数据以及待矫正的未来时期GCM逐日数据；2)计算所述气象因子历史逐日实测数据和历史GCM逐日数据各月均值，并计算二者多年月均值；3)基于历史逐日实测数据多年月均值与历史GCM逐日数据多年月均值，构建均值校正因子；4)待矫正的未来GCM逐日数据分别乘以所述均值校正因子得到均值校正后的气象序列；5)基于所述步骤1)提取的气象因子的原始数据，分别计算该气象因子历史逐日实测数据和历史GCM逐日数据各月方差，并计算多年月方差均值；6)构建方差校正因子；7)构建逐日气象因子转换函数。

【技术特征摘要】
1.一种考虑逐日数据波动特征的GCM校正方法，其特征在于，包括以下步骤：1)提取某一种气象因子的历史时期逐日实测数据和同历史时期GCM逐日数据以及待矫正的未来时期GCM逐日数据；2)计算所述气象因子历史逐日实测数据和历史GCM逐日数据各月均值，并计算二者多年月均值；3)基于历史逐日实测数据多年月均值与历史GCM逐日数据多年月均值，构建均值校正因子；4)待矫正的未来GCM逐日数据分别乘以所述均值校正因子得到均值校正后的气象序列；5)基于所述步骤1)提取的气象因子的原始数据，分别计算该气象因子历史逐日实测数据和历史GCM逐日数据各月方差，并计算多年月方差均值；6)构建方差校正因子；7)构建逐日气象因子转换函数。2.根据权利要求1所述的一种考虑逐日数据波动特征的GCM校正方法，其特征在于，所述步骤2)中，气象因子历史逐日实测数据多年月均值的计算式如下：其中，Mo，i为历史逐日实测数据的第i月多年均值，i＝1，2，3...12；Xo，i，k，j历史时期第i月第k年第j天实测值，k＝1，2，3...s，j＝1，2，3...，n，n＝28/29/30/31，具体n的值以该月的实际天数为准，s为年数；所述历史GCM逐日数据多年月均值的计算式如下：其中，Mm，i为历史GCM逐日数据的第i月多年均值，i＝1，2，3...12；Xm，i，k，j历史时期第i月第k年第j天GCM值，k＝1，2，3...s；j＝1，2，3...，n，n＝28/29/30/31，具体n的值以该月的实际天数为准，s为年数。3....

【专利技术属性】
技术研发人员：王卫光，丁一民，鲍金丽，邢万秋，董青，傅健宇，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32