趙建華，張 強(qiáng)，周廣勝

（1．甘肅?。袊?guó)氣象局干旱氣候變化與減災(zāi)重點(diǎn)實(shí)驗(yàn)室，中國(guó)氣象局蘭州干旱氣象研究所，甘肅 蘭州 730020；2．中國(guó)氣象科學(xué)研究院，北京 100081）

基于回路系統(tǒng)的干旱研究

趙建華1，張 強(qiáng)1，周廣勝2

（1．甘肅省／中國(guó)氣象局干旱氣候變化與減災(zāi)重點(diǎn)實(shí)驗(yàn)室，中國(guó)氣象局蘭州干旱氣象研究所，甘肅 蘭州 730020；2．中國(guó)氣象科學(xué)研究院，北京 100081）

干旱是一種損失最大、影響人口最多的自然災(zāi)害。陸面過(guò)程中，白天水分通常從土壤輸送到植被與大氣中，夜間相反，水分又部分回流到土壤中。這種水分在土壤、植被和大氣之間的流動(dòng)可以看作是一個(gè)回路，而構(gòu)成該回路的土壤、植被和大氣3個(gè)部件通過(guò)水分的紐帶作用構(gòu)成了一個(gè)表象上的完整開(kāi)放系統(tǒng)，稱(chēng)為回路系統(tǒng)。由此系統(tǒng)可以對(duì)氣象干旱、農(nóng)業(yè)干旱和水文干旱進(jìn)行統(tǒng)一研究，干旱是該系統(tǒng)內(nèi)在矛盾運(yùn)動(dòng)狀況的外在反映。對(duì)該回路系統(tǒng)的分析表明，該系統(tǒng)主要通過(guò)內(nèi)外2個(gè)過(guò)程與溫度控制過(guò)程維持。基于含水量、水分流量和熱量3個(gè)物理量，給出回路各部件與系統(tǒng)統(tǒng)一的水分熱動(dòng)力學(xué)方程組，并給出描述無(wú)旱澇過(guò)程的正態(tài)方程與描述旱澇過(guò)程的差量方程。最后，在含水量不變、流量不變、無(wú)植被的裸地、植被覆蓋密集的地方以及水分運(yùn)動(dòng)停止等特殊情況下對(duì)干旱進(jìn)行理論分析。結(jié)果表明：（1）對(duì)于干旱半干旱地區(qū)或無(wú)植被的裸地，干旱發(fā)生與否取決于水分源差量S′的符號(hào)；（2）對(duì)于濕潤(rùn)區(qū)或植被覆蓋度很大的地區(qū)，干旱發(fā)生與否取決于S′和?W′／?t的符號(hào)，其中?W′／?t是含水量差量的時(shí)間偏導(dǎo)數(shù)。隨后簡(jiǎn)單討論了干旱檢測(cè)的問(wèn)題。總之，以回路系統(tǒng)和水分熱動(dòng)力學(xué)方程組可以對(duì)干旱進(jìn)行系統(tǒng)和定量研究，這對(duì)于干旱的基礎(chǔ)理論研究、干旱模式以及干旱的預(yù)測(cè)具有重要意義。

干旱；回路系統(tǒng)；內(nèi)外過(guò)程；水分熱動(dòng)力學(xué)方程組

趙建華，張 強(qiáng)，周廣勝．基于回路系統(tǒng)的干旱研究［J］．干旱氣象，2016，34（6）：1087－1098，［ZHAO Jianhua，ZHANG Qiang，ZHOU Guangsheng．Drought Research Based on the Circulation System in Land Surface Process［J］．Journal of Arid Meteorology，2016，34（6）：1087－1098］，DOI：10．11755／j．issn．1006－7639（2016）－06－1087

引 言

干旱是損失最大、影響人口最多的自然災(zāi)害之一。全球干旱、半干旱地區(qū)總面積占陸地面積的34．9%，85%的耕地靠自然降水，42．9%的耕地位于干旱、半干旱區(qū)，干旱、半干旱現(xiàn)象日趨嚴(yán)重［1］。干旱每年造成的經(jīng)濟(jì)損失達(dá)數(shù)千億美元［2］，它對(duì)全球糧食、水資源、沙漠化、濕地、生物多樣性、社會(huì)經(jīng)濟(jì)、環(huán)境及人類(lèi)可持續(xù)發(fā)展具有嚴(yán)重威脅［2－3］。近幾十年來(lái)，我國(guó)干旱事件頻繁發(fā)生，如1998年北方發(fā)生的大范圍干旱［4］，2000年春夏北方又發(fā)生大面積的干旱［5］，2006年重慶遭遇百年一遇的特大伏旱［4］，2009年秋季到2010年春季西南出現(xiàn)的特大干旱［4，6］，2011年北方發(fā)生60 a一遇大旱、長(zhǎng)江中下游5省區(qū)出現(xiàn)50 a罕見(jiàn)的春夏連旱［7］等，這對(duì)我國(guó)人畜飲水、糧食安全和社會(huì)經(jīng)濟(jì)產(chǎn)生了重要影響。隨著氣候變暖的加劇、人口增長(zhǎng)以及工農(nóng)業(yè)的發(fā)展，水資源日趨短缺，加之極端天氣事件的周期縮短［8］，干旱事件呈現(xiàn)明顯上升趨勢(shì)［9－10］，其危害也日趨嚴(yán)重［3，10］。因此，研究干旱具有重要意義。

干旱很復(fù)雜，國(guó)內(nèi)外對(duì)于干旱的定義多達(dá)100多種［11］。世界氣象組織與美國(guó)氣象學(xué)會(huì)一般將干旱定義為“長(zhǎng)時(shí)間缺乏降水或降水明顯短缺”或“由于缺少降水，異常干的天氣時(shí)段持續(xù)到引起水文失衡”等［12］；王勁松等［13］定義干旱為“降水的短缺或當(dāng)水的供應(yīng)不能滿足對(duì)水的需求時(shí)的狀況，是指由水分收支或供求不平衡所形成的水分短缺現(xiàn)象”；張強(qiáng)等［14］認(rèn)為“干旱是指某一地域在某一具體時(shí)段內(nèi)的降水量比多年平均降水量顯著偏少，導(dǎo)致該地區(qū)的經(jīng)濟(jì)活動(dòng)（尤其是農(nóng)業(yè)生產(chǎn)）和人類(lèi)生活受到較大危害的現(xiàn)象。它是一種氣候?yàn)?zāi)害，也是一種持續(xù)性的氣象災(zāi)害”。干旱一般包括2種含義：一種是水分收入小于支出（Aridity），即干旱氣候，側(cè)重于空間；另一種是可利用水小于均值（Drought），即干旱災(zāi)害，側(cè)重于時(shí)間［11，15－16］。顯然，上述定義均有不同側(cè)重，其中王勁松等［13］給出的定義最貼近本質(zhì)。張強(qiáng)等［2］更進(jìn)一步認(rèn)為干旱一詞實(shí)際上包括了干旱氣候、干旱和干旱災(zāi)害3種含義。這種多學(xué)科交叉和含義的復(fù)合性，是干旱定義復(fù)雜的2個(gè)重要原因。

國(guó)內(nèi)外一般從成因、規(guī)律、影響、預(yù)報(bào)與監(jiān)測(cè)方法等方面對(duì)干旱展開(kāi)研究。如成因方面，下沉運(yùn)動(dòng)被認(rèn)為是重要原因［2，11，13－14，17］，而氣候變化［18－19］、北大西洋濤動(dòng) （NAO）［20－21］、ENSO［22－25］、北極濤動(dòng)［24］、海溫變化［24，26－27］、海陸氣之間的相互作用［28－29］、環(huán)流或風(fēng)場(chǎng)異常［9，30－32］以及人類(lèi)活動(dòng)［32－33］等也被認(rèn)為是干旱發(fā)生的原因；干旱規(guī)律方面，主要從干旱的空間分布、時(shí)間演變以及周期性等方面展開(kāi)研究［34－38］；干旱影響方面，主要研究了干旱對(duì)天氣氣候［39－40］以及植物［41］的影響；干旱預(yù)報(bào)方面，采用了包括水文模式［19］、氣候模式［28，42－43］、陸面模式［44－45］、旱澇持續(xù)性的尺度因子［46］、游程理論（Run Theory）［47］、自回歸模型［48］、馬爾科夫鏈模型［49］、神經(jīng) 網(wǎng) 絡(luò) 方 法［23，50］、條 件殘差重采樣技術(shù)（Conditional Residual Resample Technique）［51］、模糊邏輯模型［23］、低階離散自回歸滑動(dòng)平均（DARMA）模型［52］和泊松過(guò)程模式［53］等方法進(jìn)行研究；干旱監(jiān)測(cè)主要以干旱指數(shù)實(shí)現(xiàn)。據(jù)統(tǒng)計(jì)，全球提出了上百種干旱指數(shù)［2］，實(shí)際可能更多，常用的如帕爾默干旱指數(shù)、標(biāo)準(zhǔn)化降水指數(shù)、降水距平百分率、土壤相對(duì)濕度指數(shù)和CI指數(shù)等。其中，帕爾默干旱指數(shù)被認(rèn)為是干旱指數(shù)發(fā)展史上的轉(zhuǎn)折點(diǎn)。近年來(lái)，又產(chǎn)生了許多新的指數(shù)，比如快變指數(shù)（Rapid Change Index）［54］、農(nóng)業(yè)干旱參考指數(shù)［30］、歸一化差異水體指數(shù)［55］、K指數(shù)［56］、土壤濕度虧缺指數(shù)（Soil Moisture Deficit Index，SMDI）［57］、蒸散虧缺指數(shù)（Evapotranspiration Deficit Index，ETDI）［57］以及多尺度干旱指數(shù)標(biāo)準(zhǔn)化降水蒸散指數(shù)（Standardized Precipitation Evapotranspiration Index，SPEI）［58］等。

干旱的這種定義多、種類(lèi)多、成因復(fù)雜、規(guī)律多樣、影響面廣、危害大以及監(jiān)測(cè)手段眾多的特點(diǎn)，表明干旱并不是局部事件（如大氣、土壤、植被或社會(huì)經(jīng)濟(jì)等方面），它應(yīng)該站在系統(tǒng)或整體的角度上來(lái)理解，這樣才能對(duì)干旱進(jìn)行統(tǒng)一定義與研究；而從局部出發(fā)只能導(dǎo)致目前這種頭緒紛雜的局面。干旱主要是發(fā)生在土壤、植被和大氣邊界層中的一種現(xiàn)象，陸面過(guò)程對(duì)其具有重要意義。水文模式［19］和陸面模式［44－45］在干旱研究與預(yù)測(cè)中的應(yīng)用就體現(xiàn)了其重要性。許多研究指出，陸氣相互作用存在一個(gè)耦合態(tài)［59］。孫菽芬［60］指出，水熱在土壤、植被和大氣中的輸送存在一個(gè)物理回路過(guò)程。雖然也有許多文獻(xiàn)從理論方面對(duì)干旱展開(kāi)了研究［61－62］，但是從系統(tǒng)角度、物理回路過(guò)程出發(fā)進(jìn)行研究的很少見(jiàn)到。本文從回路過(guò)程出發(fā)，通過(guò)分析干旱的矛盾本質(zhì)和回路維持的機(jī)制，給出回路系統(tǒng)熱動(dòng)力學(xué)方程組，然后在一定條件下對(duì)之進(jìn)行理論分析。這對(duì)干旱的基礎(chǔ)研究具有重要價(jià)值。

1 刻畫(huà)干旱的回路過(guò)程

一般意義上的干旱主要發(fā)生在大氣邊界層下面。實(shí)際物理過(guò)程中水分在土壤、植被與空氣之間進(jìn)行遷移。以晴天為例，考慮一個(gè)土壤、植被和大氣邊界層構(gòu)成的理想封閉系統(tǒng)：白天水分自土壤中通過(guò)根系吸收到植被或者通過(guò)蒸發(fā)進(jìn)入到底層空氣中，然后植被中的水分通過(guò)蒸騰、底層空氣中的水分通過(guò)湍流輸送被帶至上層空氣中；夜間，地表變冷，空氣中的部分水分又通過(guò)地表或植被凝結(jié)／凝華（比如降露水）回流到土壤中，如圖1a所示。需要指出的是，降露水是陸面水分過(guò)程的重要分量，在濕潤(rùn)區(qū)與雨季，露水在促進(jìn)植物生長(zhǎng)中的作用幾乎可以忽略，而在干旱、半干旱地區(qū)則不然，露水是重要的水資源和濕度來(lái)源，它作為土壤的一種額外水分輸入，一定程度上緩解了土壤水分短缺，在某些干旱區(qū)降露水甚至超過(guò)了降水量，是干旱區(qū)植物生長(zhǎng)極為重要的液態(tài)水來(lái)源［63－65］。顯然，水分完成了一個(gè)表象上完整的日循環(huán)過(guò)程。該日循環(huán)為一個(gè)物理回路過(guò)程［60］，這個(gè)回路可看作是以水分為紐帶，且連通土壤、植被和大氣邊界層的回路系統(tǒng)。以此回路系統(tǒng)即可以對(duì)氣象干旱、農(nóng)業(yè)干旱與水文干旱進(jìn)行統(tǒng)一研究。當(dāng)然社會(huì)經(jīng)濟(jì)干旱也可以加入，此處為了突出重點(diǎn)，略。

當(dāng)然實(shí)際情況是回路系統(tǒng)很難封閉：水分會(huì)在輸送中被空氣和土壤中的水平運(yùn)動(dòng)攜帶至其它地方，從而不能回到原地；不過(guò)這并不影響回路的定義，因?yàn)槲覀儾皇菑腖agrange觀點(diǎn)研究問(wèn)題，而是從Euler觀點(diǎn)出發(fā)的，回路仍然是存在的，只不過(guò)流通的水分具有不同的背景。另外，回路中的水分也可能完成不了一個(gè)完整的循環(huán)，比如在極端干旱區(qū)，水分可能從空氣中回流不到土壤，這也不影響回路的定義，因?yàn)檫@種情況在數(shù)學(xué)上只不過(guò)表示流量為0，是回路中斷現(xiàn)象?；芈废到y(tǒng)是開(kāi)放系統(tǒng)。

陸面過(guò)程中常參考電路中計(jì)算阻抗的方法處理湍流通量在土壤、植被和空氣之間的輸送［60］，此處也采用之，這樣該回路系統(tǒng)就可以如圖1b所示。圖1b中，S代表土壤，其底部可取為研究時(shí)段內(nèi)土壤濕度沒(méi)有顯著變化的層次；V代表植被；A代表底層空氣，與植被同高度；T代表底層空氣以上至某一參考層的空氣，該參考層可以取到混合層頂或殘留層頂?shù)母叨?，這樣在一定程度上可保證水分守恒，因?yàn)樵跊](méi)有天氣過(guò)程并忽視部件水平方向、土壤底部與邊界層頂部的水分供給，以及忽視化學(xué)過(guò)程的條件下，水分主要分布在回路中。值得注意的是，圖1b是簡(jiǎn)略圖，并聯(lián)的V和A之間實(shí)際上是有水分交流的，為了簡(jiǎn)潔，圖中沒(méi)有繪出；而且T和S是有一定高度的，為了美觀，繪成了水平。孫菽芬認(rèn)為土壤是陸地表面水文循環(huán)的緩沖儲(chǔ)水池［60］，T顯然也具有這種功能，它既能蓄水，又能供水。

圖1 水分遷移回路系統(tǒng)示意圖（a）晝夜示意圖；（b）回路系統(tǒng)（a中虛線為晝夜分界線；b中帶箭頭的圓環(huán)表示輸送方向，實(shí)線代表水分自S輸送到T的過(guò)程，虛線反之）Fig．1 Sketches of water flows in circulation system （a）water flows at daytime and nighttime，（b）water flows in whole system （The dashed line in fig．a(chǎn) divides a day into daytime and nighttime．Circle with arrows represents the direction of transportation，solid lines in fig．b represent the processwater flows from S to T，while dashed lines represent the reverse）

2 干旱的矛盾本質(zhì)

基于上述認(rèn)識(shí)，以及對(duì)干旱反映水分收支或供求不平衡的普遍理解，將上述回路系統(tǒng)進(jìn)行條理化：（1）回路系統(tǒng)由土壤、植被和空氣等一系列部件構(gòu)成，干旱作為一種現(xiàn)象，是系統(tǒng)內(nèi)在運(yùn)動(dòng)狀況的外在反映；（2）作為干旱關(guān)鍵與決定因素的水分，是系統(tǒng)各部件連接的紐帶；（3）干旱是需水與供水之間的矛盾所致。這就是對(duì)干旱的認(rèn)識(shí)，對(duì)于澇漬的認(rèn)識(shí)類(lèi)似。

系統(tǒng)中的水分在系統(tǒng)內(nèi)部各個(gè)部件之間進(jìn)行輸送交流，滿足系統(tǒng)各部件的水分需求和系統(tǒng)的維持。正常狀況（即無(wú)旱澇）下，需水與水分供應(yīng)是既充足又平衡的。所謂的正常狀況，指的是一地長(zhǎng)期穩(wěn)定的自然狀況，是包括植物、動(dòng)物、人、土壤、空氣等在內(nèi)的多個(gè)方面平衡適應(yīng)的一種狀況，它不需要針對(duì)某個(gè)部件，而是系統(tǒng)的、整體的狀況。干旱就是針對(duì)此正常狀況產(chǎn)生的。特殊情況下，當(dāng)系統(tǒng)某些部件需水出現(xiàn)變化（比如植被生長(zhǎng)、氣溫變化、土壤水分下滲等引起需水增大或者減?。r(shí)，兩者之間的平衡會(huì)被打破；如果水分輸送不能適應(yīng)需水的變化，矛盾出現(xiàn)，旱澇就會(huì)在該部件上發(fā)生，直至兩者之間的矛盾被解決或調(diào)和為止；如果兩者不適應(yīng)時(shí)間較長(zhǎng)或者很難適應(yīng)，進(jìn)而還會(huì)引起該部件上下游（以水分輸送的方向決定）的其它部件出現(xiàn)“病變”，甚至可能導(dǎo)致系統(tǒng)整體出現(xiàn)病變，乃至瓦解、崩潰或重構(gòu)。從另一角度而言，當(dāng)水分供應(yīng)發(fā)生變化（加快、減慢或者停止），部件需水也會(huì)相應(yīng)出現(xiàn)變化；如果需水的調(diào)整不適應(yīng)水分供應(yīng)的變化，則矛盾出現(xiàn)，部件就會(huì)出現(xiàn)水分積累或虧損，旱澇發(fā)生，直至需水與水分供應(yīng)達(dá)到正常狀況且相互適應(yīng)為止。如果兩者不適應(yīng)時(shí)間較長(zhǎng)，部件病變會(huì)成為常態(tài)，整個(gè)系統(tǒng)都會(huì)異常。顯然，需水與供水之間的矛盾運(yùn)動(dòng)是旱與澇的內(nèi)因。需水與水分供應(yīng)之間的矛盾體現(xiàn)為水分在部件上出現(xiàn)了超出正常需水的大量積累為澇；反之為旱。

3 含水量和水分輸送的控制過(guò)程

3．1 含水量的控制

回路系統(tǒng)中，部件往往具有一定的持水能力，具體表現(xiàn)為部件含水量。含水量一般是有限度的，有最大值和最小值。部件含水量與部件自身的物理或生理特性密切相關(guān)。比如，空氣的氣溫、密度以及氣壓也是控制含水量的因素。風(fēng)速主要起控制水分流動(dòng)快慢（即水分流量，簡(jiǎn)稱(chēng)為流量）的作用。

外在因素也能起到控制部件含水量的作用，它可由回路系統(tǒng)中的其它部件實(shí)現(xiàn)，比如土壤含水量就可因植被與空氣的含水量或狀態(tài)的變化而變化，空氣含水量可因植被含水量的變化而變化。這本質(zhì)上反映了部件之間的相互聯(lián)系。

在土壤、植被和空氣3個(gè)部件中，空氣最具活性，其含水量的極值與空氣狀態(tài)有關(guān)。空氣的這種活性，使得它在回路系統(tǒng)中具有重要地位，在一定條件下，它可以通過(guò)自身狀態(tài)的變化調(diào)節(jié)回路系統(tǒng)中的水分運(yùn)動(dòng)，比如通過(guò)升溫，其所容納的水分會(huì)增多，植被與土壤中的水分就會(huì)補(bǔ)充到其中，從而引發(fā)或加劇后兩者的干旱；通過(guò)降溫，降低其水分需求，從而使水分保留在植被與土壤中，緩解乃至于解除干旱。

3．2 水分輸送的控制

水分輸送控制很重要。一般認(rèn)為水分是順梯度輸送的，此時(shí)系統(tǒng)中水分的輸送主要決定于部件兩端壓強(qiáng)或水勢(shì)。這是對(duì)部件而言。對(duì)回路系統(tǒng)，該系統(tǒng)是如何維持的？比如白天土壤S的水勢(shì)如果比空氣高，則夜間如何又從空氣倒流到土壤中，水分的供應(yīng)機(jī)制為何？對(duì)圖1所示的回路系統(tǒng)，主要有以下3種過(guò)程：

（1）內(nèi)部供水過(guò)程。通過(guò)部件內(nèi)在動(dòng)力系統(tǒng)與部件狀態(tài)變化驅(qū)動(dòng)，比如典型的蒸散—凝結(jié)／凝華過(guò)程：①土壤供水，通過(guò)臨近部件的植被（吸收水分）和空氣（通過(guò)湍流、升溫、減壓、體積減小以及粉塵氣溶膠數(shù)量增多或粒徑變小等過(guò)程）產(chǎn)生正水勢(shì)梯度，促使水分由S到A和V的輸送。即白天蒸散過(guò)程；②空氣供水，通過(guò)空氣降溫、升壓、水汽輻合、體積收縮以及粉塵氣溶膠數(shù)量減少與粒徑變大等過(guò)程促使水汽凝結(jié)，進(jìn)而回流到土壤，如夜間凝結(jié)。在這2種過(guò)程中，空氣的特性呈現(xiàn)出一種相反的變化，通過(guò)這種類(lèi)似于交流電的反相變化便可以實(shí)現(xiàn)回路系統(tǒng)的正常工作，簡(jiǎn)稱(chēng)此過(guò)程為內(nèi)過(guò)程。顯然，內(nèi)過(guò)程具有典型的日周期或年周期，可以在周期內(nèi)的各個(gè)位相或周期平均上考察干旱情況，即會(huì)得到相應(yīng)的小時(shí)、日、月、季或年尺度上的干旱情況。

（2）外部供水過(guò)程。理想的水分供應(yīng)應(yīng)該是除了能滿足回路系統(tǒng)中各個(gè)部件的正常需水之外，還能保證回路系統(tǒng)的正常運(yùn)行。但在實(shí)際中，由于內(nèi)過(guò)程的總水量守恒，各個(gè)部件的需水量常常變化，比如植被生長(zhǎng)需要水分；當(dāng)內(nèi)過(guò)程產(chǎn)生不了足夠的供應(yīng)或者不能維持回路正常工作時(shí)，干旱就可能發(fā)生。此時(shí)外部供水顯得非常重要，它的發(fā)生與否將是干旱緩解與加重的主要因素。這就是外部供水過(guò)程，簡(jiǎn)稱(chēng)外過(guò)程。外部供水的典型例子是降水與灌溉，當(dāng)然空氣抬升與下沉也是重要過(guò)程。外部供水發(fā)生的時(shí)機(jī)和程度是該過(guò)程非常重要的2個(gè)物理量，具體表現(xiàn)為降水或灌溉出現(xiàn)的時(shí)間與其水分補(bǔ)給量的大小。有效的外部供水應(yīng)該是在彌補(bǔ)前期內(nèi)部供水造成的水分缺失前提下，還能促進(jìn)回路系統(tǒng)的健康運(yùn)行。

（3）溫度控制過(guò)程。溫度對(duì)土壤中的水熱耦合運(yùn)動(dòng)［60］、植被生長(zhǎng)和大氣相對(duì)濕度具有重要的控制作用。通過(guò)溫度變化，控制土壤中熱量輸送引起的水分輸送過(guò)程，進(jìn)而對(duì)土壤中含水量產(chǎn)生影響；通過(guò)溫度變化，控制植被的生長(zhǎng)，從而加快或減緩植被對(duì)水分的吸收、輸送以及儲(chǔ)存；通過(guò)溫度變化，改變大氣飽和水汽壓與相對(duì)濕度，從而控制大氣中水分的含量與輸送。故溫度可以對(duì)部件含水量和水分輸送產(chǎn)生一定的控制作用，這種作用進(jìn)而對(duì)上下游部件也產(chǎn)生一定影響。

回路系統(tǒng)在內(nèi)過(guò)程維持下發(fā)生干旱時(shí)，如果沒(méi)有外過(guò)程參與，則：①系統(tǒng)中的水分需求旺盛部件會(huì)截留或榨取其它需求不旺盛部件所需水分，從而使回路系統(tǒng)病態(tài)發(fā)展，如高溫下植被的枯萎，空氣榨取了植被中的水分。或者局部部件雖然水分需求不旺盛，但是由于出現(xiàn)干旱，阻礙了水分的流通，也會(huì)造成回路系統(tǒng)病態(tài)發(fā)展，比如土壤干旱，這是局部干旱；②系統(tǒng)中各部件均一致缺水，造成回路系統(tǒng)中流量一致降低，各部件均干旱，這是系統(tǒng)干旱。另外，由局部干旱與系統(tǒng)干旱可以對(duì)干旱進(jìn)行等級(jí)劃分。對(duì)于系統(tǒng)干旱，以部分部件出現(xiàn)缺水為輕旱；以回路系統(tǒng)水分運(yùn)動(dòng)中斷、內(nèi)過(guò)程不能維持回路系統(tǒng)為中旱；以回路系統(tǒng)所有部件均無(wú)水分供應(yīng)為重旱。對(duì)于局部干旱，可稱(chēng)有水分供應(yīng)、但低于部件正常需要為輕旱；稱(chēng)無(wú)水分供應(yīng)、且部件含水量“滲出”為中旱；稱(chēng)既無(wú)水分供應(yīng)、又無(wú)部件含水量“滲出”為重旱。

不管是內(nèi)過(guò)程還是外過(guò)程，都離不開(kāi)系統(tǒng)與部件需水量。需水量對(duì)內(nèi)過(guò)程的維持、外過(guò)程的效果具有重要影響，也對(duì)干旱指數(shù)和干旱的預(yù)報(bào)、預(yù)警、緩解及消除具有重要影響。

總而言之，研究干旱需要研究以下問(wèn)題：①回路系統(tǒng)與部件需水量；②水分輸送；③需水量與水分輸送之間的矛盾；④回路系統(tǒng)維持機(jī)制與內(nèi)外供水過(guò)程；⑤溫度的控制作用等。

4 基于回路系統(tǒng)的干旱水分熱動(dòng)力學(xué)方程組

4．1 含水量、流量和熱量方程

含水量和流量是回路系統(tǒng)中的2個(gè)關(guān)鍵物理量，溫度也是重要的控制因素，因此研究干旱需要考慮這3個(gè)物理量。土壤、植被和大氣各自有其相關(guān)的含水量、流量與熱量的方程，為了統(tǒng)一研究，此處將它們合寫(xiě)為一個(gè)方程。不過(guò)注意的是，水分在流經(jīng)部件時(shí)其質(zhì)量可能是變化的，比如降水、灌溉、植被吸收、植被枯萎、植被受傷后流出傷流量、土壤中根系吸收、霧和枝葉上的露水、溫度變化產(chǎn)生的凍融過(guò)程或者吸濕水含量發(fā)生變化等，尤其是水分被植被吸收、植被氣孔關(guān)閉的情況下，水分近乎“消失”，不考慮質(zhì)量變化顯然不符合事實(shí)，也會(huì)引起誤差。因此含水量和流量的方程需要考慮這種變化。通常熱量平衡方程采用的是以溫度為預(yù)報(bào)變量的形式，但這樣難以統(tǒng)一各個(gè)部件，為此仍以熱量為預(yù)報(bào)變量給出其方程。另外需要指出的是，回路系統(tǒng)中，由于植被對(duì)水分的吸收、水分在空氣中的重新分布以及在土壤中的滲透，在回路系統(tǒng)中運(yùn)行的水分實(shí)際上很難保證是同一質(zhì)量，不過(guò)這并不影響基于Euler觀點(diǎn)對(duì)方程的推導(dǎo)。

設(shè)單位體積總含水量（包括液態(tài)、氣態(tài)與固態(tài)）為W（單位：kg·m－3），流量為I（單位：kg·s－1· m－2），熱量為Q（單位：J·m－3）。水分可以擴(kuò)散、被吸收或者消失，比如化學(xué)與生物過(guò)程。

對(duì)某一部件，其質(zhì)量變化等于流入量與其質(zhì)量變化率之和，其流量變化等于邊界流入量和外力強(qiáng)迫之和，其熱量變化等于流入量與外源項(xiàng)（包括轉(zhuǎn)化項(xiàng)）之和，容易得到土壤、植被和大氣3部件統(tǒng)一的含水量和流量方程如下：

其中，下標(biāo)i代表第i個(gè)部件，Wi、Ii與Qi分別是第i個(gè)部件的含水量、流量與熱量，土壤與大氣中存在水熱耦合效應(yīng)，Ii即是其綜合效應(yīng)，wi是水分垂直遷移速度。注意的是，這里定義的I和W都是復(fù)合項(xiàng)；Hi是由熱傳導(dǎo)、湍流以及液態(tài)水和氣態(tài)水運(yùn)動(dòng)引起的熱量在垂直方向的輸送；DhW，i、DhI，i與DhQ，i分別是水平方向的含水量散度、流量散度、熱量散度，DhW，i比如土壤與空氣中的水分輻合輻散，DhI，i比如空氣中水平方向的水汽通量散度，DhQ，i比如因空氣輻合輻散、土壤受到擠壓以及液態(tài)水與氣態(tài)水在橫向流動(dòng)時(shí)產(chǎn)生的熱量在水平方向的變化；δi是熱量的外源項(xiàng)，包括做功、耗散、動(dòng)能的轉(zhuǎn)化、輻射、水的相變以及光合作用等產(chǎn)生的熱量變化；顯然，（3）式是個(gè)極端綜合的熱量方程，其具體表達(dá)式隨部件的不同而有變化?！鞔韺?shí)際部件的空間尺度變化；pi是部件上下兩端水分的水勢(shì)；Si是其它形式的水分源匯項(xiàng)，比如部件自身耗水量或補(bǔ)水量；li是部件的線尺度，g是重力加速度，t是時(shí)間。Wthres，i是含水量的臨界值，包括飽和值與最低限度的取值，它是時(shí)間t、溫度Ti、含水量Wi、部件質(zhì)量Mi以及其它因素Oi（對(duì)生物，比如陽(yáng)光、營(yíng)養(yǎng)物、生長(zhǎng)速率與呼吸速率等）的函數(shù)。當(dāng)Wi＝Wthres，i時(shí)，部件含水量達(dá)到飽和或最低限度，（4）式的關(guān)系會(huì)對(duì)（1）式產(chǎn)生制約。

在陸面過(guò)程中，植被和空氣中的流量I和Q并不用（2）式和（3）式計(jì)算，而是通過(guò)阻抗或總體輸送系數(shù)計(jì)算。實(shí)際計(jì)算時(shí)也可以采用這2種方法。

由（1）式可知，回路系統(tǒng)中水平方向水分補(bǔ)給、內(nèi)過(guò)程（垂向水分補(bǔ)給）、外過(guò)程（水平方向水分補(bǔ)給、降水及灌溉等其它形式水源）和部件自身耗水量是干旱的決定因素。其物理意義是：內(nèi)外過(guò)程水分供應(yīng)不足或部件自身耗水量過(guò)大，均可導(dǎo)致干旱發(fā)生。由（2）式可知，水分輸送快慢既與流量散度有關(guān)，又與部件凈水勢(shì)梯度有關(guān)。其物理意義是：正的凈水勢(shì)梯度、流量輻合、垂向流量通量?jī)糨斎?，可加速水分流?dòng)，反之減慢。（3）式亦然。

將各部件的（1）、（2）、（3）式求和，即得系統(tǒng)整體的含水量、流量和熱量的方程：為系統(tǒng)整體的含水量、流量、熱量、含水量輻合輻散項(xiàng)、流量輻合輻散項(xiàng)、熱量輻合輻散項(xiàng)、質(zhì)量外源匯和熱量外源匯。顯然，系統(tǒng)的水熱方程組與部件的在形式上并無(wú)二致，區(qū)別僅在于對(duì)象不同；（1）—（3）式針對(duì)部件而言，可對(duì)部件各自的干旱情況進(jìn)行預(yù)測(cè)，（5）—（7）式則針對(duì)整個(gè)系統(tǒng)，預(yù)測(cè)的是整體情況。（1）—（7）式即是回路系統(tǒng)中干旱水分的熱動(dòng)力學(xué)方程組。

由于回路系統(tǒng)具有典型的日周期，故可對(duì)除（4）式之外的（1）—（7）式進(jìn)行24 h的時(shí)間積分，即得系統(tǒng)含水量、流量與熱量的日變化：

其中，下標(biāo)0表示初始時(shí)刻，1表示末時(shí)刻。（8）—（13）式即可預(yù)測(cè)各部件與系統(tǒng)整體在一日的干旱狀況。

當(dāng)然，由于該回路系統(tǒng)大致以白天和黑夜分界，還可以對(duì)白天和黑夜分別進(jìn)行積分，得到相關(guān)物理量的變化。

4．2 正態(tài)方程與差量方程

部件需水異常和水分輸送異常對(duì)干旱有重要影響，因此需要研究需水和水分輸送的正常狀態(tài)和對(duì)正常狀態(tài)的偏離?，F(xiàn)實(shí)中確實(shí)在許多情況下是不發(fā)生干旱的，這就是正常狀態(tài)。如前所述，正常狀態(tài)即包括植物、動(dòng)物、人、土壤、空氣等在內(nèi)的多個(gè)方面平衡適應(yīng)的一種自然狀況，它是系統(tǒng)的、整體的一種自適應(yīng)狀況，并不針對(duì)個(gè)別部件。干旱因此也可以看成是一種含水量、流量以及兩者之間關(guān)系對(duì)正常狀態(tài)的偏離，這種偏離就是矛盾，矛盾是事物發(fā)展的內(nèi)在動(dòng)力［66］。正常狀態(tài)可以試驗(yàn)確定、也可以統(tǒng)計(jì)確定，此處給出其方程（簡(jiǎn)稱(chēng)“正態(tài)方程”）。正態(tài)方程可以采用求平均的方法獲得，不過(guò)求平均時(shí)注意應(yīng)該將旱澇的時(shí)段去掉，因?yàn)殡y以保證旱、澇兩者的平均。對(duì)（1）—（7）式求時(shí)間平均后可得正態(tài)方程如下：

類(lèi)似于求湍流方程的方法，以（1）—（7）式減去（14）—（20）式可得干旱水分熱動(dòng)力學(xué)的差量方程：

量代表差量，即對(duì)正態(tài)的偏差。

差量的意義在于可以在規(guī)律的基礎(chǔ)上認(rèn)識(shí)變化。比如其中的W′i與W′反映了旱（含水量的偏差為負(fù)）與澇（含水量的偏差為正）的程度，偏差越大，越不正常，旱澇就會(huì)加重，風(fēng)險(xiǎn)性也會(huì)增大。由（21）—（27）式即可對(duì)旱澇的成因及其演變進(jìn)行預(yù)測(cè)。

從上述方程組可以看出，差量方程與正常狀態(tài)是有關(guān)的。只要預(yù)先知道正態(tài)和初邊值，就可以知道偏差狀態(tài)。再根據(jù)各學(xué)科對(duì)干旱不同的定義，即可以對(duì)干旱進(jìn)行日、月或年尺度的預(yù)測(cè)。另外，常說(shuō)的正能量與負(fù)能量也可以根據(jù)以上差量方程得到，此略。

4．3 干旱的環(huán)量描述

對(duì)于正態(tài)，流量沿回路是單向的。發(fā)生干旱時(shí)，流量變小、中斷乃至于消失。對(duì)于系統(tǒng)干旱，可以用流量沿回路的積分來(lái)表征，通過(guò)研究此積分的變化來(lái)研究干旱。（2）式乘以△li，略經(jīng)處理后可得：

其中，Γi＝Ii△li為局部部件的環(huán)量。

（28）式沿整個(gè)回路求和，即得整個(gè)回路系統(tǒng)的環(huán)量Γ方程如下：

5 干旱熱動(dòng)力學(xué)方程的理論分析

理想情況下，含水量與流量應(yīng)該遵從（14）—（20）式，此時(shí)沒(méi)有所謂的旱澇過(guò)程。但實(shí)際上由于氣候變化與人類(lèi)活動(dòng)等原因，這種正態(tài)過(guò)程往往被打破，從而出現(xiàn)了旱澇。以下對(duì)一些簡(jiǎn)單情況進(jìn)行討論。另外，由于溫度起控制作用，為簡(jiǎn)單計(jì)算，以下的討論將不考慮它。

5．1 含水量不變

考慮一極端理想的情況，對(duì)某一部件i，含水量Wi不變，沒(méi)有外源外匯（Si＝0），并不計(jì)側(cè)邊界流入（DhW，i＝0，DhI，i＝0）。由（1）和（2）式得到：

（30）式表明，流量在部件i的兩端是相等的，即流入多少，就流出多少，水分在流動(dòng)中不被部件截留，部件含水量達(dá)到定常的飽和狀態(tài)或環(huán)境溫壓恒常、植被不生長(zhǎng)時(shí)容易出現(xiàn)此種情況。但流量是時(shí)間的函數(shù)，它和外力以及部件兩端含水量的分布有關(guān)。若再設(shè)，則流量?jī)H是凈水勢(shì)的函數(shù)。若凈水勢(shì)梯度為常數(shù)，則流量是后者的線性函數(shù)。整個(gè)回路系統(tǒng)結(jié)論相似。

5．2 流量不變

再考慮另一極端理想情況，流量不隨時(shí)間和空間變化，且Si＝0，DhW，i＝0，DhI，i＝0，則由（1）和（2）式可得：

此時(shí)，Wi與水勢(shì)梯度以及垂直遷移速度梯度有關(guān)。若水勢(shì)梯度為常數(shù)，（33）易積得：

這是個(gè)超越方程。

5．3 內(nèi)過(guò)程與外過(guò)程

內(nèi)過(guò)程造成的水勢(shì)變化與其維持回路系統(tǒng)正常運(yùn)行的時(shí)長(zhǎng)，對(duì)預(yù)報(bào)干旱具有重要意義。外過(guò)程出現(xiàn)的時(shí)機(jī)與程度亦然。先以（11）、（12）式研究之。忽略水平側(cè)邊界的水分通量與流量通量的水平梯度，即DhW＝0，DhI＝0。由（11）式可知，此時(shí)系統(tǒng)中的水分由流量的垂直梯度和外源匯決定。如果，即系統(tǒng)在一日內(nèi)無(wú)凈的水分補(bǔ)給或消耗，則水分僅取決于流量梯度；若，則回路系統(tǒng)中流通的水分除了被部件攔截一部分之外，另一部分被部件所吸收；若，外部供水、植被枯萎、作物受傷或化學(xué)過(guò)程流出的水分，該水分可以部分維持部件需水與水分流通。

內(nèi)過(guò)程的停止意味著流通的水分不足以維持系統(tǒng)的正常需水，對(duì)一個(gè)一致變化的過(guò)程，其臨界狀態(tài)為。外過(guò)程亦然，如果其供給的水分不足以緩解干旱，也滿足該條件。以（25）、（26）式研究此兩過(guò)程。內(nèi)外過(guò)程一般需要進(jìn)行一段時(shí)間的積分，不過(guò)積分后方程內(nèi)核并沒(méi)有什么變化，此處為了簡(jiǎn)略，直接以不積分的方程進(jìn)行討論。

5．3．1 無(wú)植被的裸地

為簡(jiǎn)單計(jì)，考慮干旱半干旱區(qū)或植被非常稀疏的地區(qū)。此時(shí)植被可忽視，回路系統(tǒng)是由土壤與空氣2個(gè)部件以串聯(lián)形式構(gòu)成的，水分的輸送與吸釋僅發(fā)生在土壤和空氣之中：在晴天，白天水分自土壤中進(jìn)入到空氣中被存儲(chǔ)，夜間又通過(guò)凝結(jié)或凝華過(guò)程返回到土壤中。

設(shè)S′＝0，且△l′＝0，后者的意義是：忽略植被的異常生長(zhǎng)、大氣邊界層高度的日變化以及水分在土壤中輸送深度的變化等；S′＝0對(duì)內(nèi)過(guò)程而言是部件狀態(tài)為正態(tài)，對(duì)外過(guò)程而言是供給的水分處于正態(tài)之下。則（25）式化為：

其中，下標(biāo)A表示空氣部件，S表示土壤部件，見(jiàn)圖2。上式推導(dǎo)中已經(jīng)應(yīng)用等式I′S2＝I′A1和I′S3＝I′A4，因?yàn)樗鼈冋锰幱谙噜彶考慕唤缑?，并設(shè)I′A2＝I′A3和I′S1＝I′S4，當(dāng)然如果空氣部件T和土壤部件S取得足夠厚，這兩項(xiàng)可以認(rèn)為是0。一般而言，水分在空氣中輸送的距離大于在土壤中輸送的距離，即△＞△。

圖2 空氣與土壤構(gòu)成的回路系統(tǒng)Fig．2 Two－component circulation system composed by atmosphere and soil

（36）式表明，對(duì)于干旱、半干旱地區(qū)或無(wú)植被的裸地，在一定條件下的一段時(shí)間內(nèi)，只要土壤到空氣的流量偏差大于空氣到土壤中的流量偏差，即凈流入到空氣中的水分比正態(tài)下的多，就可能導(dǎo)致內(nèi)過(guò)程破壞與外過(guò)程失效，進(jìn)而可能引發(fā)干旱。該過(guò)程的一個(gè)典型例子是持續(xù)高溫，因?yàn)楦邷貢?huì)促進(jìn)蒸發(fā)，容易出現(xiàn)這種情況。

如果S′＜0，比如化學(xué)過(guò)程或氣粒轉(zhuǎn)化過(guò)程，或者外過(guò)程補(bǔ)給的水分低于正態(tài)，則（36）式的條件可以更寬松，即：

此時(shí)在一定條件下，I′S2可以小于I′S3，即凈流入空氣中的水分小于正態(tài)時(shí)，也可能使內(nèi)外過(guò)程破壞或失效。這是容易理解的，因?yàn)橥獠抗┧陀谡Ｋ?，自然發(fā)生干旱的條件會(huì)降低。

如果S′＞0，對(duì)內(nèi)過(guò)程而言，代表部件其它形式的水分被逸出過(guò)多或其它來(lái)源的水分供給過(guò)多，比如土壤吸濕水被蒸發(fā)過(guò)多；對(duì)外過(guò)程而言，供水超出了正態(tài)，也可能出現(xiàn)干旱，這與流量的變化以及部件的尺度有很大關(guān)系，不過(guò)這種情況出現(xiàn)干旱的概率較小。（36）式的約束更為苛刻，其結(jié)果形式同（37）式，此略。

另外，對(duì)于外過(guò)程，其供給的水分何時(shí)能夠緩解干旱是一個(gè)重要問(wèn)題。若以?xún)?nèi)過(guò)程遭到破壞為干旱的臨界狀態(tài)，則其臨界狀態(tài)為＝0，且。顯然，只需將（36）和（37）式“≥”改為“＜”即可，具體分析類(lèi)似，此略。

5．3．2 植被全覆蓋的地表

此時(shí)對(duì)應(yīng)濕潤(rùn)區(qū)或植被覆蓋度很大的地區(qū)。其回路系統(tǒng)為土壤、植被和空氣3個(gè)部件組成的系統(tǒng)，水分的輸送發(fā)生在這3個(gè)部件中。同樣設(shè)△l′＝0，則（25）式化為：

其中，下標(biāo)V代表植被，已利用等式I′S2＝I′V1、I′V2＝I′A1、I′A4＝I′V3和I′V4＝I′S3。

（39）式的意義是回路系統(tǒng)中，上升支差量的總流量梯度不小于下沉支差量的總流量梯度時(shí)，這相當(dāng)于大氣中的次級(jí)環(huán)流，其向上的水分輸送大于向下的水分輸送，內(nèi)過(guò)程難以維持，外過(guò)程易于失效。特殊地，若設(shè)I′A2＝I′A3和I′S1＝I′S4，則相對(duì)于正態(tài)，土壤和空氣部件中額外流出的水分不低于植被額外截留的水分時(shí)，也會(huì)導(dǎo)致內(nèi)外過(guò)程破壞或失效，進(jìn)而因缺水而引發(fā)干旱。顯然，這是差量代表的類(lèi)似于次級(jí)環(huán)流的次級(jí)內(nèi)過(guò)程出現(xiàn)外耗性輸送所致。

當(dāng)S′≠0時(shí)，（38）式變?yōu)椋?/p>

該式相比于（39）式，多了S′，其意義類(lèi)似，只是受S′的影響。

5．4 水分運(yùn)動(dòng)停止

流量也是診斷干旱的重要物理量之一。在特殊情況下，系統(tǒng)中水分運(yùn)動(dòng)可以停止，如旱情較重時(shí)，作物氣孔會(huì)關(guān)閉。設(shè)某一部件初始時(shí)刻流量為正，某時(shí)段部件出口的流量遞減至0。用數(shù)學(xué)表示就是：，下標(biāo)0表示初始時(shí)刻，下標(biāo)1表示部件的入口端，2表示出口端。仍不計(jì)水平側(cè)邊界的通量，則（2）式變?yōu)椋?/p>

該式左端為流量通量，右端的意義是重力勢(shì)與水壓勢(shì)之差。整個(gè)（41）式的意義是流量通量過(guò)小而不足以克服凈水勢(shì)時(shí)，水分輸送得到抑制，部件得以?xún)?chǔ)存水分。

6 干旱檢測(cè)方法略議

干旱發(fā)生與否，完全可以用前面推導(dǎo)的干旱熱動(dòng)力學(xué)方程確定，但在實(shí)際應(yīng)用中存在一定的不便，故仍需要建立干旱指數(shù)來(lái)診斷干旱。

干旱可用含水量、流量與熱量確定。根據(jù)干旱的定義與要素，從局部和整體兩方面可以提出一系列干旱指數(shù)：

（a）基于含水量定義干旱指數(shù)。比如基于土壤濕度、相對(duì)濕度與葉片含水量等。可以部件實(shí)際含水量與正常需水之比或之差定義局部干旱指數(shù)；以系統(tǒng)需水量與內(nèi)過(guò)程所需水量的比值或差定義為系統(tǒng)干旱指數(shù)。也可以定義外部干旱指數(shù)為：降水或灌溉等外部供水量與前期內(nèi)過(guò)程缺失量之差再除以系統(tǒng)需水量或當(dāng)下內(nèi)過(guò)程維持所需水量；

（b）基于流量的指數(shù)，比如流量偏離正常值的程度，可用各種數(shù)學(xué)形式表現(xiàn)；

（c）基于水分源匯項(xiàng)的指數(shù)。水分源匯的變化，對(duì)部件需水和回路流量有重要影響，比如基于降水與蒸散定義指數(shù)；

（d）基于部件兩端水勢(shì)的指數(shù)，比如基于葉水勢(shì)或?qū)嶋H水壓與其正常值建立指數(shù)；

（e）基于部件阻抗的指數(shù)。阻抗增大超過(guò)限度，比如阻抗無(wú)窮大，此時(shí)水流消失，部件吸納水分，典型例子是不下雨；

（f）部件承載水分的容量有限，即它具有最大持水容量（比如田間持水量），故可以借助持水量定義干旱指數(shù)；

（g）基于內(nèi)過(guò)程持續(xù)時(shí)間的指數(shù)，比如內(nèi)過(guò)程遭破壞、而外過(guò)程又不發(fā)生時(shí)，持續(xù)的時(shí)間越長(zhǎng)，干旱越嚴(yán)重；

（h）基于部件狀態(tài)或性質(zhì)變化構(gòu)建指數(shù)，比如基于作物形態(tài)、作物生理、溫度變化、風(fēng)速變化、地表地下徑流變化等。

常見(jiàn)的干旱指數(shù)，如年降水量、降水距平百分率、PDSI（Palmer Drought Severity Index）、SPI（Standardized Precipitation Index）、降水與氣溫比值（干燥指數(shù)）、降水與蒸散比值（干燥指數(shù)）以及Thornthwaite水分指數(shù)法等即是這些方面的反映，這些指數(shù)各有側(cè)重。

另外，內(nèi)過(guò)程的維持與破壞對(duì)干旱有重要影響。依前面定義，可以?xún)?nèi)過(guò)程破壞為輕旱，它此時(shí)可用（36）、（37）、（39）和（40）等式確定，也可用（11）和（12）式、（25）和（26）式通過(guò)診斷的方法確定。

用（11）與（12）式研究時(shí)，只需判斷它們的左右兩端是否相等即可。如果相等，表明回路系統(tǒng)運(yùn)行正常，干旱出現(xiàn)概率較低；如果左端大于右端，表明回路系統(tǒng)運(yùn)行異常，干旱可能出現(xiàn)；如果左端小于右端，回路系統(tǒng)運(yùn)行依然異常，澇漬可能出現(xiàn)。

用（25）與（26）式研究時(shí)，需要注意：正態(tài)過(guò)程的W′、I′、S′等接近0，內(nèi)過(guò)程破壞意味著出現(xiàn)了W′、I′、S′。如果以日尺度研究干旱，則可以對(duì)（25）與（26）式進(jìn)行24 h積分，根據(jù)W′和I′的符號(hào)與其出現(xiàn)的日期很容易判斷干旱何日出現(xiàn)。（25）與（26）式的24 h積分形式如下：

若W′＜0和I′＜0，則傾向于出現(xiàn)旱情；若W′＞0和I′＞0，則傾向于出現(xiàn)澇漬；若W′I′＜0，則旱澇均有可能發(fā)生。另外，可以通過(guò)診斷這個(gè)方程，以確定旱澇發(fā)生的原因。

內(nèi)過(guò)程的破壞意味著外過(guò)程的發(fā)生具有一定的必要性。若內(nèi)過(guò)程破壞時(shí)，外過(guò)程沒(méi)有發(fā)生，比如干旱期間，若降水與灌溉沒(méi)有出現(xiàn)，則干旱會(huì)持續(xù)；對(duì)于澇漬，如果缺乏徑流或高壓系統(tǒng)等有效的排水因素，澇漬也將持續(xù)。反之，內(nèi)過(guò)程破壞時(shí)，即使有外過(guò)程發(fā)生，則旱澇可能就不再出現(xiàn)。

當(dāng)然，對(duì)于系統(tǒng)干旱，還可用（28）、（29）式的環(huán)量方程確定。當(dāng)環(huán)量符號(hào)相反或者絕對(duì)值減少時(shí)，即認(rèn)為系統(tǒng)干旱。另外只要系統(tǒng)回路中出現(xiàn)局部流量中斷，也即認(rèn)為出現(xiàn)了系統(tǒng)干旱。

7 結(jié) 論

（1）土壤、植被和空氣可構(gòu)成一個(gè)開(kāi)放的回路系統(tǒng)，該回路系統(tǒng)可以對(duì)氣象干旱、農(nóng)業(yè)干旱與水文干旱進(jìn)行統(tǒng)一研究。

（2）內(nèi)過(guò)程、外過(guò)程與溫度控制過(guò)程是回路系統(tǒng)維持的3個(gè)重要機(jī)制。

（3）推導(dǎo)出了回路系統(tǒng)各部件與系統(tǒng)整體的干旱水分熱動(dòng)力學(xué)方程組，并給出了描述正常無(wú)旱澇狀態(tài)的正態(tài)方程組和描述旱澇過(guò)程的差量方程組。

（4）理論分析顯示，對(duì)于干旱、半干旱地區(qū)或無(wú)植被裸地，如果S′＝0，在一定條件下，只要凈流入到空氣中的水分比正態(tài)下的多，就可導(dǎo)致內(nèi)過(guò)程破壞與外過(guò)程失效，進(jìn)而可能引發(fā)干旱；如果S′＜0，很容易發(fā)生干旱，在一定條件下，凈流入空氣中的水分小于正態(tài)時(shí)，也可能使內(nèi)外過(guò)程破壞或失效；如果S′＞0，不易發(fā)生干旱。對(duì)于濕潤(rùn)區(qū)或植被蓋度很大的地區(qū)，當(dāng) S′＝0時(shí)，若回路系統(tǒng)中上升支的總流量梯度不小于下沉支的總流量梯度時(shí)，易引發(fā)干旱。

梯度小于下沉支的總流量梯度時(shí)，易引發(fā)干旱。

當(dāng)然，本文沒(méi)有對(duì)熱量過(guò)程進(jìn)行分析。而且，推導(dǎo)的熱動(dòng)力學(xué)方程組也沒(méi)有考慮風(fēng)的貢獻(xiàn)，這些需要以后進(jìn)一步展開(kāi)研究。

致謝：誠(chéng)摯感謝中國(guó)科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所錢(qián)正安教授對(duì)本文做出的細(xì)心修改與指導(dǎo)，誠(chéng)摯感謝中國(guó)氣象局蘭州干旱氣象研究所陽(yáng)伏伶博士、王勁松研究員和張良博士以及蘭州大學(xué)生命科學(xué)學(xué)院鄭紅俊的幫助。

［1］劉永忠，李齊霞，孫萬(wàn)榮，等．氣候干旱與作物干旱指標(biāo)體系［J］．山西農(nóng)業(yè)科學(xué)，2005，33（3）：50－53．

［2］張強(qiáng)，張良，崔顯成，等．干旱監(jiān)測(cè)與評(píng)價(jià)技術(shù)的發(fā)展及其科學(xué)挑戰(zhàn)［J］．地球科學(xué)進(jìn)展，2011，26（7）：763－778．

［3］陳亞寧，楊青，羅毅，等．西北干旱區(qū)水資源問(wèn)題研究思考［J］．干旱區(qū)地理，2012，35（1）：1－9．

［4］崔冬林．近50年中國(guó)區(qū)域持續(xù)性氣象干旱事件的客觀識(shí)別及變化研究［D］．蘭州：蘭州大學(xué)，2010．

［5］衛(wèi)捷，張慶云，陶詩(shī)言．1999及2000年夏季華北嚴(yán)重干旱的物理成因分析［J］．大氣科學(xué)，2004，28（1）：125－137．

［6］王素萍，段海霞，馮建英．2009／2010年冬季全國(guó)干旱狀況及其影響與成因［J］．干旱氣象，2010，28（1）：107－112．

［7］王素萍，段海霞，馮建英．2011年春季全國(guó)干旱狀況及其影響與成因［J］．干旱氣象，2011，29（2）：261－268．

［8］Mishra A K，Singh V P．A review of drought concepts［J］．Journal of Hydrology，2010，391：202－216．

［9］Zhang LX，Zhou T J．Droughtover east Asia：A review［J］．Journal of Climate，2015，28（8）：3375－3399．

［10］Sivakumar M V K，Wilhite D A，Pulwarty R S，etal．The High－Level Meeting on National Drought Policy［J］．Bulletin of the A-merican Meteorological Society，2014，95（4）：ES85－ES88．

［11］張書(shū)余．干旱氣象學(xué)［M］．北京：氣象出版社，2008：1－55．

［12］Richard R，Heim J．A review of Twentieth－Century drought indices used in the United States［J］．Bulletin of the American Meteorological Society，2002，83（8）：1149－1165．

［13］王勁松，李耀輝，王潤(rùn)元，等．我國(guó)氣象干旱研究進(jìn)展評(píng)述［J］．干旱氣象，2012，30（4）：497－508．

［14］張強(qiáng)，王潤(rùn)元，鄧振鏞．中國(guó)西北干旱氣候變化對(duì)農(nóng)業(yè)與生態(tài)影響及對(duì)策［M］．北京：氣象出版社，2012：1－62．

［15］Mainguet M．Aridity：Droughts and Human Development［M］．Berlin：Springer－Verlag，1999：27．

［16］Kassas M．Aridity，drought and desertification［M］．Lebanon：Arab environment future challenges．Technical Publications and Environment＆Development Magazine，2008：96－110．

［17］Namias J．Some causes of United States drought［J］．Journal of Applied Meteorology，1983，22（1）：30－39．

［18］Dai X G，Liu Y，Wang P．Warm－dry collocation of the recent drought in southwestern China tied to moisture transport and climate warming［J］．Chinese Physics B，2015，24（4）：546－554．

［19］Lehner B，D?ll P，Alcamo J，etal．Estimating the impactof global change on flood and drought risks in Europe：a continental，integrated analysis［J］．Climatic Change，2006，75（3）：273－299．

［20］Wong W K，Beldring S，Engenskaugen T，et al．Climate Change Effects on Spatiotemporal Patterns of Hydroclimatological Summer Droughts in Norway［J］．Journal of Hydrometeorology，2011，12 （6）：1205－1220．

［21］Xin X G，Yu R C，Zhou T J，etal．Drought in late spring of South China in recent decades［J］．Journal of Climate，2006，19（13）：3197－3206．

［22］Cook E R，Seager R，Cane M A，et al．North American drought：Reconstructions，causes，and consequences［J］．Earth Science Reviews，2007，8（1）：93－134．

［23］Ummenhofer C C，Gupta A S，Briggs PR．Indian and Pacific O-cean Influences on Southeast Australian Drought and Soil Moisture ［J］．Journal of Climate，2011，24（5）：1313－1336．

［24］Barriopedro D，Gouveia C M，Trigo R M，et al．The 2009／10 Drought in China：Possible Causes and Impacts on Vegetation［J］．Journal of Hydrometeorology，2012，13（4）：1251－1267．

［25］?zger M，Mishra A K，Singh V P．Long Lead Time Drought Forecasting Using a Wavelet and Fuzzy Logic Combination Model：A Case Study in Texas［J］．Journal of Hydrometeorology，2012，13 （1）：284－297．

［26］Hoerling M，Eischeid J，Perlwitz J，et al．On the Increased Frequency of Mediterranean Drought［J］．Journal of Climate，2012，25（6）：2146－2161．

［27］Skinner C B，Ashfaq M，Diffenbaugh N S．Influence of Twenty－First－Century Atmospheric and Sea Surface Temperature Forcing on West African Climate［J］．Journal of Climate，2012，25（2）：527－542．

［28］Schubert SD，Suarez M J，Pegion P J，et al．Causes of longterm drought in the U．S．Great Plains［J］．Journal of Climate，2004，17（3）：485－503．

［29］Kingston D G，F(xiàn)leig A K，Tallaksen L M，et al．Ocean－Atmosphere Forcing of Summer Streamflow Drought in Great Britain［J］．Journal of Hydrometeorology，2011，14（1）：331－344．

［30］Woli P，Jones JW，Ingram K T，et al．Agricultural reference index for drought（ARID）［J］．Agronomy Journal，2009，104（2）：287－300．

［31］Seager R，Goddard L，Nakamura J，etal．DynamicalCausesof the 2010／11 Texas－Northern Mexico Drought［J］．Journal of Hydrometeorology，2014，15（1）：39－68．

［32］Delworth TL，Zeng F，RosatiA，etal．A Link between the Hiatus in Global Warming and North American Drought［J］．Journal of Climate，2015，28（9）：3834－3845．

［33］Bagley JE，Desai A R，Harding K J，etal．Droughtand Deforestation：Has Land Cover Change Influenced Recent Precipitation Extremes in the Amazon？［J］．Journal of Climate，2014，27（1）：345－361．

［34］馬柱國(guó)．華北干旱化趨勢(shì)及轉(zhuǎn)折性變化與太平洋年代際振蕩的關(guān)系［J］．科學(xué)通報(bào)，2007，52（10）：1199－1206．

［35］馬柱國(guó)，符淙斌．20世紀(jì)下半葉全球干旱化的事實(shí)及其與大尺度背景的聯(lián)系［J］．中國(guó)科學(xué)D輯，2007，37（2）：222－233．

［36］侯威，楊萍，封國(guó)林．中國(guó)極端干旱事件的年代際變化及其成因［J］．物理學(xué)報(bào)，2008，57（6）：3932－3940．

［37］王春林，鄒菊香，麥北堅(jiān)，等．近50年華南氣象干旱時(shí)空特征及其變化趨勢(shì)［J］．生態(tài)學(xué)報(bào)，2015，35（3）：595－602．

［38］李韻婕，任福民，李憶平，等．1960—2010年中國(guó)西南地區(qū)區(qū)域性氣象干旱事件的特征分析［J］．氣象學(xué)報(bào)，2014，72（2）：266 －276．

［39］Wang Y，Xie Y，Cai L，etal．Impact of the 2011 Southern U．S．Drought on Ground－Level Fine Aerosol Concentration in Summertime［J］．Journal of the Atmospheric Sciences，2015，72（3）：1075 －1093．

［40］Winguth A M，Kelp B．The Urban Heat Island of the North－Central Texas Region and Its Relation to the 2011 Severe Texas Drought［J］．Journal of Applied Meteorology＆ Climatology，2013，52（11）：2418－2433．

［41］吳芹，張光燦，裴斌，等．3個(gè)樹(shù)種對(duì)不同程度土壤干旱的生理生化響應(yīng)［J］．生態(tài)學(xué)報(bào)，2013，33（12）：3648－3656．

［42］Poitras V，Sushama L，Seglenieks F，et al．Projected Changes to Streamflow Characteristics over Western Canada as Simulated by the Canadian RCM［J］．Journal of Hydrometeorology，2011，12 （6）：1395－1413．

［43］Mo K C，Long L N，Schemm JE．Characteristics of drought and persistentwet spells over the United States in the atmosphere－land －ocean coupled model experiments［J］．Earth Interactions，2012，16（16）：1－26．

［44］Zhang L，Zhang H Q，Zhang Q，et al．Comparison of a land surfacemodel and Palmer drought severity index inmonitoring drought in China［A］．IGARSS［C］．Canada，July 24－29，2011，58 （11）：1985－1988．

［45］Crow W T，Kumar SV，Bolten JD．On the utility of land surface models for agricultural droughtmonitoring［J］．Hydrology＆Earth System Sciences Discussions，2012，9（4）：5167－5193．

［46］龔志強(qiáng)，封國(guó)林．中國(guó)近1000年旱澇的持續(xù)性特征研究［J］．物理學(xué)報(bào)，2008，57（6）：3920－3931．

［47］MoyéL A，Kapadia A S，Cech IM，etal．The theory of runswith applications to drought prediction［J］．Journal of Hydrology，1988，103（1－2）：127－137．

［48］Kendall D R，Dracup JA．On the generation of drought events using an alternating renewal－reward model［J］．Stochastic Envi-ronmental Research and Risk Assessment，1992，6（1）：55－68．

［49］Akyuz D E，BayazitM，Onoz B．Markov Chain Models for Hydrological Drought Characteristics［J］．Journal of Hydrometeorology，2012，13（1）：298－309．

［50］Mishra A K，Desai V R．Drought forecasting using feed－forward recursive neutral network［J］．Ecological Modeling，2006，198（1 －2）：127－138．

［51］Hwang Y，Carbone G J．Ensemble forecasts of drought indices using a conditional residual resampling technique［J］．Journal of Applied Meteorology＆Climatology，2009，48（7）：1289－1301．

［52］Chung C H，Salas JD．Drought occurrence probabilities and risks of dependent hydrological process［J］．Journal of Hydrologic Engineering，2000，5（3）：259－268．

［53］Yoon JH，Mo K，Wood E F．Dynamic－Model－Based Seasonal Prediction of Meteorological Drought over the Contiguous United States［J］．Journal of Hydrometeorology，2012，13（2）：463－482．

［54］Otkin JA，Anderson M C，Hain C，etal．Examining the Relationship between Drought Development and Rapid Changes in the E-vaporative Stress Index［J］．Journal of Hydrometeorology，2014，15 （3）：938－956．

［55］Ji L，Zhang L，Wylie B．Analysis of dynamic thresholds for the normalized difference water index［J］．Photogrammetric Engineering＆Remote Sensing，2009，75（11）：1307－1317．

［56］王勁松，郭江勇，傾繼祖．一種K干旱指數(shù)在西北地區(qū)春旱分析中的應(yīng)用［J］．自然資源學(xué)報(bào)，2007，22（5）：710－717．

［57］Narasimhan B，Srinivasan R．Development and evaluation of Soil Moisture Deficit Index（SMDI）and Evapotranspiration Deficit Index（ETDI）for agricultural droughtmonitoring［J］．Agriculture＆Forest Meteorology，2005，133（1－4）：69－88．

［58］Vicente－Serrano SM，Beguería S，López－Moreno J I．A multi －scalar drought index sensitive to globalwarming：The standardized precipitation evapotranspiration index－SPEI［J］．Journal of Climate，2010，23（7）：1696－1718．

［59］Roundy JK，Wood E F．The Attribution of Land－Atmosphere Interactions on the Seasonal Predictability of Drought［J］．Journal of Hydrometeoroly，2014，16（2）：1095－1099．

［60］孫淑芬．陸面過(guò)程的物理、生化機(jī)理和參數(shù)化模型［M］．北京：氣象出版社，2005：166－177．

［61］Shiau JT，Modarres R．Copula－based drought severity－duration －frequency analysis in Iran［J］．Meteorological Applications，2009，16（4）：481－489．

［62］Lu C，Singh V P，Guo S．Drought analysis based on copulas［J］．Journal of Hydrologic Engineering，2013（18）：797－808．

［63］問(wèn)曉梅，張強(qiáng)，王勝，等．陸面露水特征及生態(tài)氣候效應(yīng)的研究進(jìn)展［J］．干旱氣象，2008，26（4）：5－11．

［64］張強(qiáng)，問(wèn)曉梅，王勝，等．關(guān)于陸面降露水測(cè)量方法及其開(kāi)發(fā)利用研究［J］．高原氣象，2010，29（4）：1085－1092．

［65］葉有華，彭少麟．露水對(duì)植物的作用效應(yīng)研究進(jìn)展［J］．生態(tài)學(xué)報(bào)，2011，31（11）：3190－3196．

［66］毛澤東．毛澤東選集［M］．北京：人民出版社，2008：299－340．

Drought Research Based on the Circulation System in Land Surface Process

ZHAO Jianhua1，ZHANG Qiang1，ZHOU Guangsheng2

（1．Key Laboratory of Arid Climatic Changing and Reducing Disaster of Gansu Province／Chinese Meteorological Administration，Lanzhou Institute of Arid Meteorology of Chinese Meteorological Administration，Lanzhou 730020，China；2．Chinese Academy of Meteorological Sciences，Beijing 10081，China）

Drought is one of the natural disasters with the biggest losses and affectingmost persons．In the land surface process，water usually transports from soil into plantand atmosphere in the daytime and then it partly reflows into soil atnight．Thiswater transporting process in soil，vegetable and atmosphere can be considered as a circuit and the three components of soil，vegetable and atmosphere linked with water flows comprise an open system in representation called circulation system．Themeteorological drought，hydrological drought and agricultural drought except socioeconomic drought can be studied in a whole by using it，and drought is an external reflection of its inner contradictionmovement ofwater．After the detailed analysis of the system，it shows that inner process，outer process and temperature－controlling process is three importantmechanisms to support the circulation system．Next，the thermaland dynamical water equations for the system and its each component，comprising three physical variableswhich are water content，water flux and heat quantity，are deduced in an uniform form．The equations for normal－state without drought or waterlogging and for abnormal－state with drought or waterlogging are deduced too．Finally，the theoretical analysis are also carried out to analyze the drought under the specified conditions such as constantwater content，constantwater flux，bare land withoutvegetable，land covered with dense vegetable and ceased water flows．The results are as follows：（1）For arid and semi－arid areas or bare land without vegetable，the trends to emerge drought depend on the signs ofwhichmeanswater extraneous source differing from its normal－state value．（2）Forwetarea or dense－vegetable covering area，the trends to appear drought depend on signs of both and meaning the time partial derivative ofwater content difference from the normal state．At last，the building of drought index is discussed simply．In short，the drought can be systematically and quantitatively studied with the circulation system and thermal and dynamical water equations，which has an important significance for the basic theoretical search on drought，drought numericalmodel and drought's forecast．

drought；circulation system；inner and outer process；thermal and dynamical water equations of drought

1006－7639（2016）－06－1087－12

10．11755／j．issn．1006－7639（2016）－06－1087

P423，P426．61，P425．6

A

2016－07－08；改回日期：2016－10－26

國(guó)家自然科學(xué)基金（41375019）與國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃（2013CB430200，2013CB430206）共同贊助

趙建華（1976－），男，寧夏中衛(wèi)人，主要從事大氣湍流、干旱與陸面過(guò)程研究．E－mail：53277357＠qq．com