利用GNSS垂向位移研究云南地區(qū)水文干旱特征

朱偉剛 馬 東 鄭 宇 楊興海

1 長(zhǎng)春工程學(xué)院勘查與測(cè)繪工程學(xué)院，長(zhǎng)春市同志街3066號(hào)，130021 2 西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院，成都市犀安路999號(hào)，611756

干旱指數(shù)作為監(jiān)測(cè)和評(píng)價(jià)旱情的重要工具，可以直觀、準(zhǔn)確、定量地表征干旱強(qiáng)度、時(shí)空分布等特征[1]。傳統(tǒng)干旱指數(shù)的計(jì)算依賴(lài)于地面單基站觀測(cè)與遙感反演等手段，如以降水與地表徑流為基礎(chǔ)的代表性指標(biāo)Z指數(shù)、標(biāo)準(zhǔn)化降雨指數(shù)SPI(standardized precipitation index)和標(biāo)準(zhǔn)化徑流指數(shù)SRI(standardized runoff index)等[2-3]。但基于單基站觀測(cè)數(shù)據(jù)計(jì)算的干旱指數(shù)受數(shù)據(jù)質(zhì)量和站點(diǎn)覆蓋面限制，難以監(jiān)測(cè)大尺度的干旱事件；基于傳統(tǒng)光學(xué)或微波遙感數(shù)據(jù)反演的干旱指數(shù)僅包含地表水、樹(shù)冠水與淺層地下水變化導(dǎo)致的氣象干旱事件，難以表征陸地水總儲(chǔ)量變化引起的水文干旱事件。

GRACE/GRACE-FO系列衛(wèi)星可觀測(cè)地球重力場(chǎng)變化，監(jiān)測(cè)地表水資源儲(chǔ)量變化情況，估計(jì)陸地水儲(chǔ)量盈虧。因此，可利用GRACE衛(wèi)星數(shù)據(jù)反演等效水高并進(jìn)一步計(jì)算干旱指數(shù)，識(shí)別并定量表征水文干旱特征[4-5]。然而，由于GRACE實(shí)際空間分辨率較低(300 km×300 km)，對(duì)小區(qū)域的陸地水儲(chǔ)量變化不夠敏感，不適合小區(qū)域水文干旱監(jiān)測(cè)；同時(shí)，GRACE任務(wù)結(jié)束和GRACE-FO任務(wù)開(kāi)始之間存在數(shù)據(jù)空白，會(huì)影響對(duì)水文干旱的評(píng)估效果。因此，需要尋找一種新的干旱監(jiān)測(cè)手段作為補(bǔ)充。

GNSS垂向位移時(shí)間序列包含環(huán)境負(fù)載(非潮汐大氣負(fù)載、非潮汐海洋負(fù)載、水文負(fù)載等)引起的地球彈性形變信息，并且對(duì)小區(qū)域的環(huán)境負(fù)載變化敏感[6]?？鄢浅毕髿庳?fù)載和非潮汐海洋負(fù)載等影響得到僅包含陸地水文負(fù)載信號(hào)的GNSS垂向位移時(shí)間序列，可以用來(lái)反演陸地水儲(chǔ)量變化及監(jiān)測(cè)因陸地水儲(chǔ)量變化引起的干旱事件[7]。國(guó)內(nèi)學(xué)者在GNSS垂向位移監(jiān)測(cè)干旱方面開(kāi)展大量研究[8-9]。但這些研究?jī)H將GNSS垂向位移與GRACE-DSI、CI(comprehensive meteorological drought index)、PDSI(Palmer drought severity index)等干旱指數(shù)進(jìn)行對(duì)比分析，未將其轉(zhuǎn)化為干旱指數(shù)并提供干旱嚴(yán)重程度的定量表征。

隨著中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)的建立，10多年來(lái)連續(xù)運(yùn)行的GNSS基準(zhǔn)站為水文大地測(cè)量提供了寶貴數(shù)據(jù)。本文從GNSS原始時(shí)間序列中提取僅包含水文負(fù)載信號(hào)的垂向位移時(shí)間序列，然后計(jì)算GNSS垂向位移的標(biāo)準(zhǔn)化異常數(shù)據(jù)GNSS-DSI，并將GNSS-DSI與GRACE-DSI進(jìn)行對(duì)比，最后利用GNSS-DSI分析云南地區(qū)的干旱特征，為干旱監(jiān)測(cè)與評(píng)估提供新的工具。

1 研究區(qū)域與數(shù)據(jù)來(lái)源

1.1 研究區(qū)域

云南省位于我國(guó)西南地區(qū)，地勢(shì)北高南低，海拔落差懸殊，氣候兼具低緯氣候、季風(fēng)氣候、山原、高原氣候等特點(diǎn)。全省包含長(zhǎng)江、珠江、元江、瀾滄江、怒江、大盈江6大水系，年降水量充沛。雖然云南水資源量很大，但在時(shí)空分布上極不均衡，并且受大氣環(huán)流異常、水汽輸送減少、青藏高原和社會(huì)環(huán)境等因素影響，云南省干旱頻發(fā)[10]。因此，開(kāi)展干旱監(jiān)測(cè)對(duì)主動(dòng)抗旱減災(zāi)、旱情風(fēng)險(xiǎn)管理與生態(tài)環(huán)境保護(hù)具有重要作用。

1.2 數(shù)據(jù)來(lái)源

1.2.1 GNSS坐標(biāo)時(shí)間序列

從國(guó)家地震科學(xué)數(shù)據(jù)中心獲取云南地區(qū)26個(gè)GNSS臺(tái)站2011～2019年的時(shí)間序列進(jìn)行分析，站點(diǎn)分布如圖1所示，其結(jié)果由GAMIT/GLOBK軟件解算得到，該產(chǎn)品已在數(shù)據(jù)解算過(guò)程中扣除固體潮、海潮、極潮等因素的影響[11]。GNSS時(shí)間序列中包含由構(gòu)造運(yùn)動(dòng)引起的趨勢(shì)性變化、由環(huán)境負(fù)載引起的季節(jié)性變化和地殼運(yùn)動(dòng)與儀器更換引起的階躍。其中，季節(jié)性變化主要體現(xiàn)在垂直方向，主要影響因素為大氣、非潮汐海洋和陸地水儲(chǔ)量負(fù)載[12]。因此，本文僅采用GNSS時(shí)間序列中的垂直分量進(jìn)行分析。

圖1 云南地區(qū)GNSS站垂向位移周年振幅及臺(tái)站分布Fig.1 Annual vertical displacement amplitude and distribution of GNSS stations in Yunnan

1.2.2 GRACE及降雨數(shù)據(jù)

從美國(guó)得克薩斯大學(xué)空間研究中心CSR(http:∥www2.csr.utexas.edu/grace/)下載GRACE/GRACE-FO RL06 V02 Mascon數(shù)據(jù)集。該數(shù)據(jù)產(chǎn)品已考慮橢球效應(yīng)，并應(yīng)用所有適當(dāng)?shù)男拚珻20、C30替換、地心一階項(xiàng)改正、冰川均衡改正(GIA)、加入大氣海洋模型(GAD)、扣除冰后回彈、非潮汐海洋負(fù)載和大氣負(fù)載影響等[13]。本文直接應(yīng)用該數(shù)據(jù)集，以研究區(qū)每個(gè)站點(diǎn)為圓心、300 km為半徑的格網(wǎng)值計(jì)算站點(diǎn)處等效水高(EWH)時(shí)間序列，并將其轉(zhuǎn)換為月尺度GRACE-DSI[14]。除2017-06～2018-06外，缺失的GRACE數(shù)據(jù)采用樣條插值方法進(jìn)行補(bǔ)充，該方法通過(guò)對(duì)鄰近點(diǎn)數(shù)值進(jìn)行三次插值得到。降雨數(shù)據(jù)從全球降水氣候?qū)W中心(http:∥gpcc.dwd.de)獲取，時(shí)間分辨率為1個(gè)月，空間分辨率為0.25°×0.25°。

2 研究方法

2.1 GNSS水文負(fù)載信號(hào)提取

為提取僅包含陸地水文負(fù)載變化引起的地表彈性形變信息，使用德國(guó)地學(xué)研究中心GFZ(https:∥www.gfz-potsdam.de/en/home/)提供的環(huán)境荷載產(chǎn)品消除非潮汐海洋負(fù)載(non-tidal ocean loading,NTOL)和非潮汐大氣負(fù)載(non-tidal atmospheric loading,NTAL)效應(yīng)的影響。從圖2可以看出，GNSS垂向位移時(shí)間序列受非潮汐大氣負(fù)載效應(yīng)影響較大，周年振幅在2～3 mm之間；受非潮汐海洋負(fù)載效應(yīng)影響較小，周年振幅最大為0.5 mm。最后對(duì)垂向位移時(shí)間序列進(jìn)行建模，基于最小二乘法估計(jì)所有參數(shù)并去除時(shí)間序列中的長(zhǎng)期線性趨勢(shì)及跳變，以減小構(gòu)造運(yùn)動(dòng)(如地震、火山噴發(fā)等)等其他因素的影響，得到以水文負(fù)載信號(hào)為主導(dǎo)的垂向位移時(shí)間序列[15]。模型公式為：

圖2 GNSS站點(diǎn)NTAL與NTOL垂向位移周年振幅Fig.2 Annual vertical displacement amplitude of NTAL and NTOL at GNSS stations

y(ti)=a+bti+csin(2πti)+dcos(2πti)+

esin(4πti)+fcos(4πti)+

(1)

為保證計(jì)算干旱指數(shù)時(shí)數(shù)據(jù)的完整性，使用GMIS軟件補(bǔ)齊缺失的GNSS數(shù)據(jù)。該軟件可以很好地恢復(fù)GNSS時(shí)間序列中缺失的數(shù)據(jù)，并且減弱時(shí)間序列中白噪聲的影響[16]。以YNCX、YNHZ、YNMH站為例，使用GMIS軟件得到濾波與插值補(bǔ)齊后的時(shí)間序列(圖3)。圖中黑色為扣除長(zhǎng)期趨勢(shì)、階躍、非潮汐海洋負(fù)載和非潮汐大氣負(fù)載等影響的GNSS垂向位移時(shí)間序列，紅色為濾波與插值補(bǔ)齊后的垂向位移時(shí)間序列。

圖3 部分GNSS站濾波與插值補(bǔ)齊后的垂向位移時(shí)間序列Fig.3 Vertical displacement time series of some GNSS stations after filtering and interpolation

2.2 水文干旱指數(shù)計(jì)算

GRACE對(duì)大尺度的陸地水儲(chǔ)量變化十分敏感，該優(yōu)勢(shì)促進(jìn)了其在干旱監(jiān)測(cè)領(lǐng)域的發(fā)展。Zhao等[17]通過(guò)GRACE數(shù)據(jù)反演等效水高(EWH)，計(jì)算EWH的標(biāo)準(zhǔn)化異常數(shù)據(jù)GRACE-DSI，結(jié)果表明，GRACE-DSI可以獨(dú)立描述干旱事件。但GRACE存在長(zhǎng)期數(shù)據(jù)空白，且受空間分辨率限制，會(huì)影響評(píng)估精度。受上述干旱評(píng)估方法啟發(fā)，本文僅基于GNSS垂向位移數(shù)據(jù)計(jì)算干旱強(qiáng)度指數(shù)GNSS-DSI，計(jì)算方法與GRACE-DSI相同，且不存在數(shù)據(jù)缺失。

水文負(fù)載與其引起的地球彈性形變呈負(fù)相關(guān)關(guān)系，即水文負(fù)載增加，GNSS臺(tái)站垂向位置向下運(yùn)動(dòng)；水文負(fù)載減小，GNSS臺(tái)站垂向位置向上運(yùn)動(dòng)。為方便與GRACE-EWH、GRACE-DSI進(jìn)行對(duì)比分析，本文對(duì)所有垂向位移數(shù)據(jù)進(jìn)行反符號(hào)處理，利用處理后的垂向位移數(shù)據(jù)(vertical crustal displacement, VCD)進(jìn)一步計(jì)算GNSS-DSI。GNSS-DSI定義為GNSS-VCD的標(biāo)準(zhǔn)化異常，計(jì)算方法是通過(guò)月內(nèi)平均將每日VCD時(shí)間序列重新采樣得到每月VCD時(shí)間序列，減去每月(1～12月)VCD平均值(氣候平均值)得到VCD異常，再除以VCD的標(biāo)準(zhǔn)差。計(jì)算公式如下：

(2)

為評(píng)估GNSS-DSI的可靠性，結(jié)合GRACE數(shù)據(jù)進(jìn)行對(duì)比，按照上述方法計(jì)算GRACE-DSI，并將干旱指數(shù)等級(jí)劃分成11個(gè)類(lèi)別(表1)。

表1 干旱指數(shù)等級(jí)劃分Tab.1 Classification of drought index

3 結(jié)果與分析

3.1 VCD、EWH與DSI及其相關(guān)性分析

為評(píng)價(jià)GNSS量化水文干旱程度的能力，本文對(duì)GNSS和GRACE觀測(cè)的均勻分布在云南不同地理位置(中部、東北部、東部、東南部、南部、西南部、西部、西北部)的8個(gè)臺(tái)站的VCD、EWH以及進(jìn)一步計(jì)算得到的DSI進(jìn)行相關(guān)性分析。 圖4為研究區(qū)26個(gè)GNSS站點(diǎn)處GNSS-VCD、GRACE-EWH和2種DSI的相關(guān)性系數(shù)以及位于不同相關(guān)性系數(shù)區(qū)間的臺(tái)站數(shù)量統(tǒng)計(jì)(無(wú)相關(guān)[0～0.2]、弱相關(guān)(0.2～0.4]、中等相關(guān)(0.4～0.6]、強(qiáng)相關(guān)(0.6～0.8]、極強(qiáng)相關(guān)(0.8～1])。從圖4(a)可以看出，整個(gè)云南地區(qū)的GNSS-VCD與GRACE-EWH表現(xiàn)出高度的一致性，兩者的相關(guān)性系數(shù)在0.58～0.94之間，位于中-極強(qiáng)相關(guān)區(qū)間，除東部地區(qū)存在2個(gè)臺(tái)站呈現(xiàn)中等相關(guān)外，其他地區(qū)的24個(gè)臺(tái)站均呈現(xiàn)強(qiáng)-極強(qiáng)相關(guān)。從圖4(b)、4(c)看出，云南地區(qū)整體的GNSS-DSI與GRACE-DSI具有較好的相關(guān)性，除呈無(wú)相關(guān)的3個(gè)臺(tái)站和呈弱相關(guān)的2個(gè)臺(tái)站外，其余21個(gè)臺(tái)站均呈中-強(qiáng)相關(guān)。結(jié)合圖4和圖1中水文負(fù)載位移周年振幅可以看出，云南東部、東南部、東北部地區(qū)陸地水儲(chǔ)量變化較弱，位于這3個(gè)區(qū)域的GNSS臺(tái)站季節(jié)性變化還包含其他地球物理因素[18]，導(dǎo)致這3個(gè)區(qū)域的GNSS-DSI與GRACE-DSI相關(guān)性較弱，并存在不穩(wěn)定性。

圖4 GNSS-VCD與GRACE-EWH、 GNSS-DSI與GRACE-DSI相關(guān)系數(shù)和不同相關(guān)系數(shù)區(qū)間的臺(tái)站數(shù)量統(tǒng)計(jì)Fig.4 Correlation coefficient between GNSS-VCD and GRACE-EWH, GNSS-DSI and GRACE-DSI and the number of stations in different correlation coefficient intervals

圖5為8個(gè)區(qū)域的GNSS-VCD、GRACE-EWH和2種干旱指數(shù)的時(shí)間序列以及相關(guān)性系數(shù)(藍(lán)色和黃色分別表示GRACE/GRACE-FO衛(wèi)星在2016-11～2017-06期間的數(shù)據(jù)質(zhì)量下降以及2017-07～2018-05期間的數(shù)據(jù)缺失)。從圖中可以看出，云南中部、西南部地區(qū)GNSS-VCD與GRACE-EWH、GNSS-DSI與GRACE-DSI具有很強(qiáng)的相關(guān)性，GNSS-VCD與GRACE-EWH的相關(guān)系數(shù)R值分別為0.89、0.92(圖5(a)、5(f))，GNSS-DSI與GRACE-DSI的R值分別為0.78、0.74(圖5(i)、5(n))；陸地水儲(chǔ)量年振幅較小的云南東北部、東部、東南部地區(qū)，GNSS-VCD與GRACE-EWH的R值分別為0.76、0.74、0.68(圖5(b)、5(c)、5(d))，GNSS-DSI與GRACE-DSI的R值分別為0.58、0.58、0.45(圖5(j)、5(k)、5(l))；其他3個(gè)地區(qū)的GNSS-VCD與GRACE-EWH均表現(xiàn)出良好的一致性，R值分別為0.83、0.93、0.91(圖5(e)、5(g)、5(h))，GNSS-DSI與GRACE-DSI也具有一定的相關(guān)性，R值分別為0.64、0.65、0.64(圖5(m)、5(o)、5(p))。從DSI變化趨勢(shì)來(lái)看，8個(gè)區(qū)域的2種DSI分別在2012、2013和2019年出現(xiàn)較低值，最低達(dá)到-2以下，顯示出3次較為嚴(yán)重的水文干旱事件。

圖5 8個(gè)區(qū)域的垂向位移-等效水高和干旱指數(shù)相關(guān)性分析Fig.5 Correlation analysis of vertical displacement-equivalent water height and drought index in 8 regions

由于GNSS臺(tái)站在接收信號(hào)過(guò)程中會(huì)產(chǎn)生部分隨機(jī)噪聲和虛假的季節(jié)性信號(hào)，而GRACE存在長(zhǎng)期的數(shù)據(jù)質(zhì)量惡化及數(shù)據(jù)空白，并且GRACE空間分辨率較低(300 km)，對(duì)小區(qū)域水文負(fù)載變化不及GNSS-VCD與GNSS-DSI敏感，從而影響其對(duì)水文干旱事件的評(píng)估效果。因此，上述因素導(dǎo)致2種干旱指數(shù)存在一定差異，結(jié)合2種數(shù)據(jù)更有助于分析水文干旱事件。

3.2 GNSS-DSI表征水文干旱事件

為證明GNSS-DSI表征水文干旱事件的能力，本文調(diào)查云南地區(qū)的水文干旱特征。圖6為8個(gè)區(qū)域的GNSS-VCD、GNSS-VCD氣候平均值以及GNSS-VCD異常和GNSS-DSI與降水(PRECI)異常時(shí)間序列的比較。表2為云南地區(qū)2011～2019年通過(guò)GNSS-DSI探測(cè)到的水文干旱事件及其特征。結(jié)合圖6和表2可以看出，GNSS-DSI探測(cè)到了云南地區(qū)發(fā)生的持續(xù)時(shí)間和嚴(yán)重程度不同的水文干旱事件，其中，2011-08～2012-03、2012-11～2013-04、2019-04～2019-12這3個(gè)時(shí)期的干旱事件最為顯著，尤其是從2019-04開(kāi)始持續(xù)9個(gè)月的特大干旱(D4)，各區(qū)域的VCD異常與PRECI異常均出現(xiàn)大幅下降，VCD平均虧損值和虧損峰值分別為-4.5 mm、-16.1 mm，DSI值最小為-2.43，期間東北部、東部、東南部區(qū)域從2019-08左右開(kāi)始，西部、西北部從2019-10左右開(kāi)始，干旱嚴(yán)重程度逐漸減弱。2011-08～2012-03干旱事件嚴(yán)重程度相對(duì)較弱，除2012-02左右DSI出現(xiàn)波谷，VCD虧損峰值最低為-14.0 mm，DSI值最小為-1.75，等級(jí)為極端干旱(D3)，其他時(shí)期均為中度-嚴(yán)重干旱(D1-D2)。2011-08～2011-11干旱主要集中在東北部、東部、東南部區(qū)域，2011-12左右開(kāi)始擴(kuò)散到整個(gè)云南地區(qū)，于2012-03逐漸恢復(fù)。曹言等[19]基于TRMM降雨數(shù)據(jù)也探測(cè)到本次干旱事件。2012-11～2013-04干旱事件與2011年比較類(lèi)似，但持續(xù)時(shí)間較短，DSI于2013-03左右出現(xiàn)波谷，VCD虧損峰值最低為-15.1 mm，DSI值最小為-1.96，等級(jí)為極端干旱(D3)，同時(shí)期南部、東南部的干旱強(qiáng)度相對(duì)其他區(qū)域較弱。金燕等[20]也識(shí)別到本次干旱事件。

表2 云南地區(qū)干旱事件總結(jié)Tab.2 Summary of drought events in Yunnan

除3次比較嚴(yán)重的干旱外，GNSS-DSI還探測(cè)到4次持續(xù)時(shí)間短、強(qiáng)度較弱的干旱。第1次為2015-06～09左右，主要分布在云南東部、西部、西北部(圖6(c)、6(g)、6(h)、6(k)、6(o)、6(p))，VCD異常均為負(fù)值，DSI在2015-08左右出現(xiàn)波谷，3個(gè)區(qū)域的DSI最低值分別為-1.29、-1.4、-1.29，等級(jí)為中度-嚴(yán)重干旱(D1-D2)。第2次為2016-07～10左右，主要分布在云南東部和東南部(圖6(c)、6(d)、6(k)、6(l))，本次干旱與2015年相似，DSI在2016-09左右出現(xiàn)波谷，DSI最低值分別為-1.42、-1.11，等級(jí)為中度-嚴(yán)重干旱(D1-D2)，通過(guò)查閱中國(guó)氣象局國(guó)家氣候中心干旱事件監(jiān)測(cè)信息，證實(shí)了本次干旱事件。第3次為2018-03～07左右，干旱范圍較大，除云南西南部以外其他7個(gè)區(qū)域均有干旱，DSI在2018-04出現(xiàn)波谷，其中東南部強(qiáng)度最大，VCD平均虧損值和虧損峰值分別為-6.1 mm、-12.7 mm，DSI最低值為-1.81，強(qiáng)度為極端干旱(D3)，其他地區(qū)顯示為中度-嚴(yán)重干旱(D1-D2)，于維等[21]也發(fā)現(xiàn)了本次干旱事件。第4次為2018-10～12左右，主要分布在云南東北部，2018-11 DSI達(dá)到低值-1.41，等級(jí)為嚴(yán)重干旱(D2)，VCD平均虧損值與虧損峰值分別為-1.1 mm、7.4 mm，同時(shí)期東部、東南、西北也出現(xiàn)中度干旱(D1)。Jiang等[22]利用GNSS垂向位移監(jiān)測(cè)云南地區(qū)陸地水儲(chǔ)量變化，同樣發(fā)現(xiàn)了本文除2016年以外的6次干旱事件，并指出期間陸地水儲(chǔ)量均呈虧損狀態(tài)，說(shuō)明該時(shí)期個(gè)別月份的PRECI異常雖為正值，但少量降水并未起到緩解干旱的作用，導(dǎo)致GNSS-DSI仍呈下降趨勢(shì)。

上述研究結(jié)果表明，基于GNSS-VCD的標(biāo)準(zhǔn)化異常GNSS-DSI可成功識(shí)別和量化因陸地水儲(chǔ)量虧損引起的水文干旱事件，說(shuō)明GNSS-DSI可作為水文干旱程度的評(píng)價(jià)指標(biāo)應(yīng)用于實(shí)際水文干旱監(jiān)測(cè)。

4 結(jié) 語(yǔ)

本文將GNSS-VCD與GRACE反演的GRACE-EWH進(jìn)行對(duì)比分析，發(fā)現(xiàn)GNSS-VCD可以很好地體現(xiàn)陸地水的季節(jié)性變化。基于此，本文利用GNSS-VCD進(jìn)一步計(jì)算GNSS-DSI，然后對(duì)比分析GNSS-DSI與GRACE-DSI，并利用GNSS-DSI研究云南地區(qū)水文干旱事件，得到以下結(jié)論：

1)云南省整體的GNSS-DSI與GRACE-DSI具有較好的相關(guān)性，共21個(gè)臺(tái)站呈中強(qiáng)相關(guān)，中部、南部、西南部、西部、西北部5個(gè)區(qū)域的GNSS臺(tái)站GNSS-DSI與GRACE-DSI相關(guān)性最好。東部、東南部、東北部個(gè)別GNSS臺(tái)站季節(jié)性變化包含其他地球物理因素，導(dǎo)致這些區(qū)域的GNSS-DSI與GRACE-DSI相關(guān)性較弱。

2)GNSS-DSI探測(cè)到云南地區(qū)7次水文干旱事件并可提供定量表征，期間陸地水儲(chǔ)量均呈虧損狀態(tài)，最為嚴(yán)重的是2019年干旱事件，強(qiáng)度大且持續(xù)時(shí)間長(zhǎng)。7次干旱期間降雨異常多為負(fù)值或呈下降趨勢(shì)，個(gè)別月份降水異常雖為正值，但由于降水量不足、氣候變化和人類(lèi)活動(dòng)等因素影響，陸地水儲(chǔ)量仍處于虧損狀態(tài)，導(dǎo)致GNSS-DSI為負(fù)值。

3)GNSS-DSI可以作為GRACE-DSI的補(bǔ)充，將GNSS-DSI與GRACE-DSI相結(jié)合，可以提高水文干旱監(jiān)測(cè)的時(shí)空分辨率，更全面地了解干旱特征，改善對(duì)水文干旱的預(yù)測(cè)、評(píng)估與管理。