请教关于综合气象干旱指数CI的滚动计算

午后阳光

       在计算CI的时候，M以及SPI的资料都是按月值存放，如果想实现CI的逐日滚动计算，那么将近30d资料处理成“月”，1a就多出5d，10a就有50d，月份之间就无法同期对比，错开了，这样计算的SPI值应该是有误的。请教各位高手，这种情况资料该怎么处理才能实现逐日的滚动计算？
      另外，国家气候中心后来发布的MCI，也就是改进后的CI，找了很多地方都没有相关的说明，在论坛里找到MCI的公式，但是SPIW60和MI30没有算法说明，如果有哪位大神有相关资料说明，可否分享一下，不甚感激！

sunrui_1126

   你好，对于新气象干旱综合指数构建思路和公式：前两项是降水近期影响，主要是气象干旱对农业的影响MCIA；后两项是降水长期影响，主要是气象干旱对水资源的影响MCIW。因此MCI公式可又分为两部分：MCI=MCIA+MCIW;MCIA=Ka *MCI,MCIW=KW*MCI，其中各地各月的Ka作物系数均有差异。
   你说的公式前面两项是针对干旱跳跃性发展问题，使用有效降水的概念，采用陆尔教授提出的权重平均降水WAP（Weighted Average of Precipitation）方法对前期降水进行衰减，你可以看下南信大外聘陆尔教授在逐日尺度上监测气象干旱方面的文章。
   MCI公式项SPIW60是基于60天滑动WAP降水SPI运算SPIW60 = SPI（WAP）：为了给定等级划分标准，采用SPI处理方式对WAP进行正态化处理，最后形成SPIW项。通过研究发现，一般情况下，近60天内的降水对农业影响最大，这里衰减系数取0.95,计算近60天内的权重平均降水。公式里为了便于计算，统一除去了WAP前面的0.05系数，Pn为距离当天前第n天降水量。
   至于你说的公式项MI30，MI是下发的全国相对湿润度指数，MI30就是全国近30天相对湿润度指数。

午后阳光

有研究过相关领域的前辈们可否交流指点一下？果然还是内容太冷了……

午后阳光

sunrui_1126 发表于 2016-9-19 18:34
你好，对于新气象干旱综合指数构建思路和公式：前两项是降水近期影响，主要是气象干旱对农业的影响MCIA ...

非常感谢！之前看到的MCI公式中提到MI30是累积相对湿润度指数，所以混淆了。我找找你推荐的文献再看看。另外，我目前已经计算出了CI的值，看了一些文献说新旧指数在趋势上以及统计意义上来说是基本一致，只是表征干旱程度可能有强弱区别，并且MCI较CI而言，平滑了曲线跳跃性突变性这种，那么我用CI来暂且表征干旱情况的话，会不会被认可，因为看到近几年的文献都是用的改进后的MCI，心里没底。

午后阳光

sunrui_1126 发表于 2016-9-19 18:34
你好，对于新气象干旱综合指数构建思路和公式：前两项是降水近期影响，主要是气象干旱对农业的影响MCIA ...

并且，在绘图时发现，如果是全国的160个站，干旱的区域性分布特点还是很明显，但是只画一个区域，比如一个地市，旱情的发生几乎是整个地区同步进行，没有分片区，看一些文献大多是针对省份来做分析，是不是地市的话，区域相对小了，所以画出来要么全部干旱，要么全部无旱？

午后阳光

sunrui_1126 发表于 2016-9-19 18:34
你好，对于新气象干旱综合指数构建思路和公式：前两项是降水近期影响，主要是气象干旱对农业的影响MCIA ...

最后还要请教你关于逐日滚动计算，请问你有做过吗？我目前的思路是将日资料按30d依次滚动处理成“月值”，因为我的SPI是用到了WMO的一个SPI可执行程序，对资料的形式有要求，所以只能这样。但如果资料这样处理，就会随着时间往后推移的同时，日期会错开，最后直至错开月份，这样SPI的值就不是在历史同期资料对比的情况下计算出来的，结果就有误。不知道该怎么处理这个点，谢谢！

午后阳光

sunrui_1126 发表于 2016-9-19 18:34
你好，对于新气象干旱综合指数构建思路和公式：前两项是降水近期影响，主要是气象干旱对农业的影响MCIA ...

你好，找到一篇陆尔教授关于逐日尺度上监测中国区域干旱的文献，其中给出的WAP公式如下：

其中文献中给出的a=0.95，N=60，但是和你提供的WAP公式差着分母一项，请问这个是怎么处理掉的？谢谢！

午后阳光

我的问题是不是太多了。。。。。

sunrui_1126

你好，你进行公式进一步演推一下，没有差异，然后再把文献中WAP指数进一步研究一下：

sunrui_1126

午后阳光 发表于 2016-9-20 11:04
我的问题是不是太多了。。。。。

还好，我慢慢给你交流