基于CLDAS2.0驅(qū)動數(shù)據(jù)的中國區(qū)域土壤濕度模擬與評估

師春香， 姜立鵬， 朱 智， 姜志偉， 梁 曉， 韓 帥， 張 濤

(1.南京信息工程大學(xué)地理與遙感學(xué)院，江蘇南京 210044； 2.國家氣象信息中心，北京 100081)

土壤濕度是反映陸面表面狀況的一個重要物理量[1-2]，也是農(nóng)業(yè)干旱監(jiān)測中的重要指標(biāo)[3-5]。土壤濕度在時間尺度上具有周尺度甚至月尺度的記憶能力，可直接影響短期氣候變化和中尺度天氣過程[6]；同時土壤濕度可以通過影響地表反照率和蒸散發(fā)，進而對農(nóng)作物生長產(chǎn)生影響[7]。盡管土壤濕度在農(nóng)業(yè)氣象、氣候變化等研究中非常重要，但是觀測值的有限密度和數(shù)量制約了相關(guān)研究的開展[8-9]，因此研究者通過使用陸面模式對陸表狀況進行模擬[10]，獲取時空分布連續(xù)的土壤濕度模擬結(jié)果，進而開展相關(guān)研究。

盡管研究者們在陸面模型發(fā)展與改進方面開展了大量的工作，但是陸面模式的模擬結(jié)果仍然存在一定的誤差。改善陸面模式模擬結(jié)果的途徑主要有以下4種：(1)制作更為精確的大氣驅(qū)動數(shù)據(jù)[11]；(2)更換土壤質(zhì)地等地表參數(shù)[12]；(3)改進陸面模式中的物理過程[13]；(4)同化觀測資料[14]，其中高質(zhì)量、高時空分辨率的大氣驅(qū)動數(shù)據(jù)是準(zhǔn)確的陸面過程模擬中不可缺少的。

對于陸面模式研究者來說，缺乏長時間序列的、高分辨率的、接近真實的大氣驅(qū)動數(shù)據(jù)一直是困擾他們的一大問題[15]，并且大量關(guān)于大氣驅(qū)動數(shù)據(jù)的研究也證明了準(zhǔn)確的、有更多觀測信息存在的驅(qū)動數(shù)據(jù)對于提高陸面模式模擬結(jié)果的重要性，大氣驅(qū)動數(shù)據(jù)的質(zhì)量對陸面模式真實模擬地表狀況影響很大[16-19]。近年來，許多研究團隊致力于研究全球以及區(qū)域的大氣驅(qū)動數(shù)據(jù)作為陸面模式的驅(qū)動場[20-25]。盡管國際上已有幾套全球大氣強迫數(shù)據(jù)集可供陸面過程研究使用，但在全球各地的質(zhì)量有很大的差異。雖然中國研究者最近幾年也開展過關(guān)于大氣驅(qū)動數(shù)據(jù)集的研究[26-32]，但還缺少高質(zhì)量、高時空分辨率的大氣驅(qū)動數(shù)據(jù)集供相關(guān)研究使用，可以說高質(zhì)量、高時空分辨率的中國區(qū)域大氣驅(qū)動數(shù)據(jù)的缺乏已經(jīng)變成制約陸面過程模擬及相關(guān)研究的重要因素。

筆者基于高質(zhì)量、高時空分辨率的中國氣象局陸面數(shù)據(jù)同化系統(tǒng)(CMA Land Data Assimilation System，CLDAS)2.0驅(qū)動數(shù)據(jù)，運行Noah-MP陸面模式對中國區(qū)域土壤濕度進行模擬，并根據(jù)逐小時的土壤濕度觀測值對土壤濕度模擬結(jié)果進行評估，探討在使用高質(zhì)量、高時空分辨率驅(qū)動數(shù)據(jù)的情況下，陸面模式對土壤濕度的模擬精度，試圖為農(nóng)業(yè)干旱監(jiān)測提供更為準(zhǔn)確的土壤濕度格點數(shù)據(jù)。

1 數(shù)據(jù)與方法

1.1 土壤濕度觀測值及質(zhì)量控制方法

本研究所使用的土壤濕度觀測來源于國家氣象信息中心資料服務(wù)室，時間分辨率為逐小時，時間尺度是2013—2014年。

由于目前土壤水分自動站點在處于布設(shè)-驗收-業(yè)務(wù)化的過程中，土壤濕度觀測值的質(zhì)量參差不齊，因此對于土壤濕度觀測值進行質(zhì)量控制是很有必要的。氣象數(shù)據(jù)質(zhì)量控制一般包括極限值檢查、時間一致性檢查、空間一致性檢查等步驟，本研究參考了韓帥提出的土壤濕度觀測值質(zhì)量控制方案[33]，并根據(jù)觀測值實際情況對質(zhì)量控制方案進行相應(yīng)的改動。本研究所使用的土壤濕度觀測值質(zhì)量控制方法包括以下3步：(1)由于觀測儀器的土壤體積含水量測量范圍是0～0.5 mm3/mm3，并且在土壤濕度觀測站點所在的場地不會出現(xiàn)土壤中沒有水分的情況，因此剔除了≤0或 >0.5 mm3/mm3的土壤濕度觀測值；(2)由于在0 ℃以下，土壤中的水分會以固態(tài)和液態(tài)2種形態(tài)存在，同時土壤濕度觀測儀器在0 ℃以下并不能正常工作，因此根據(jù)站點的10 cm土壤溫度觀測值對相應(yīng)時次的土壤濕度觀測值進行過濾，如果10 cm土壤溫度觀測值小于0 ℃，則剔除相應(yīng)時刻的土壤濕度觀測值；(3)由于各土壤濕度觀測站開始業(yè)務(wù)化的時間并不一致且各站點的土壤濕度觀測值質(zhì)量存在差異，須要選取有效觀測值較多的站點，因此本研究剔除了年觀測時間 <180 d 站點的所有土壤濕度觀測值。經(jīng)過質(zhì)量控制后，選取931個站點的0～10 cm土壤濕度觀測值，站點分布情況如圖1所示，由此可以看出經(jīng)過質(zhì)量控制后的土壤濕度觀測站呈現(xiàn)東南密集、西北稀疏的空間分布特征，并且由于剔除了10 cm土壤溫度<0 ℃時對應(yīng)的土壤濕度，在青藏高原地區(qū)只有零星站點在質(zhì)量控制后保留了下來。

由于在土壤濕度觀測值質(zhì)量控制過程中，北方地區(qū)的冬季土壤濕度觀測值大多數(shù)被剔除，為了在研究中與其他地區(qū)的土壤濕度情況進行對比，筆者選取了2013—2014年4—10月的0～10 cm土壤濕度作為研究對象。

1.2 CLDAS2.0驅(qū)動數(shù)據(jù)

《國家氣象科技創(chuàng)新工程(2014—2020年)實施方案》明確指出，研制中國區(qū)域千米分辨率地面氣溫、氣壓、濕度、風(fēng)、降水、輻射、土壤溫濕度、積雪等產(chǎn)品是國家氣象信息中心的任務(wù)之一。為了完成這一目標(biāo)，國家氣象信息中心進行了較長遠(yuǎn)的發(fā)展規(guī)劃，位于核心地位的是中國氣象局陸面數(shù)據(jù)同化系統(tǒng)(CLDAS)的建設(shè)，該系統(tǒng)計劃分為4個階段進行，即CLDAS1.0～CLDAS4.0。其中，CLDAS1.0系統(tǒng)目前已經(jīng)在國家氣象信息中心業(yè)務(wù)運行，其主要目標(biāo)是設(shè)計一個可擴展的陸面數(shù)據(jù)同化系統(tǒng)框架，并為系統(tǒng)升級預(yù)留接口。2015年，國家氣象信息中心氣象數(shù)據(jù)研究室開始了CLDAS2.0系統(tǒng)研發(fā)工作，其主要任務(wù)是建設(shè)一個較長時間序列(2008—2014年)的大氣驅(qū)動數(shù)據(jù)以及實現(xiàn)基于該大氣驅(qū)動數(shù)據(jù)的多個陸面模式運行，并且為建立長時間序列(1979年至今)的大氣驅(qū)動數(shù)據(jù)做準(zhǔn)備。本研究使用的CLDAS2.0驅(qū)動數(shù)據(jù)來源于國家氣象信息中心氣象數(shù)據(jù)研究室，包括近地面氣溫、氣壓、濕度、風(fēng)速、降水和短波輻射6個要素，空間分辨率為 0.062 5°，時間分辨率為1 h。

1.3 Noah-MP陸面模式

Noah-MP(the Community Noah Land Surface Model with Muti-Parameterization Options)陸面模式是CLDAS2.0系統(tǒng)添加的主要陸面模式，也是美國新一代的陸面模式。目前Noah-MP陸面模式已經(jīng)廣泛應(yīng)用于陸面過程模擬研究。Yang等利用Noah-MP陸面模式對全球50個主要流域的水文狀況進行模擬，結(jié)果表明Noah-MP陸面模式能夠較好地模擬出地表溫度、土壤濕度等重要的陸表變量[34]；Cai等利用北美陸面數(shù)據(jù)同化系統(tǒng)(North-American Land Data Assimilation Systems，NLDAS)平臺測試了4個陸面模式(Noah、Noah-MP、CLM、VIC)在水文上的表現(xiàn)，結(jié)果表明，相對于其他3個陸面模式，Noah-MP陸面模式在模擬土壤濕度以及陸地水儲量上表現(xiàn)得最好[35]。

由于表層土壤濕度對于農(nóng)作物生長以及氣候變化研究是十分重要的，因此本研究選取了0～10 cm土壤濕度作為研究對象。同時為了在區(qū)域尺度上對土壤濕度的模擬結(jié)果進行評估，參考朱晨等的研究結(jié)果[8]，并結(jié)合實際觀測站點的空間分布，在中國區(qū)域選取6個研究區(qū)(圖2)，其中Ⅰ區(qū)為東北地區(qū)，Ⅱ區(qū)為華北地區(qū)，Ⅲ區(qū)為江淮地區(qū)，Ⅳ區(qū)為東南地區(qū)，Ⅴ區(qū)為西北東部地區(qū)，Ⅵ為西南地區(qū)。

2 結(jié)果與分析

2.1 土壤濕度模擬誤差的空間分布

圖3給出了由中國區(qū)域0～10 cm土壤濕度觀測值與Noah-MP模擬值的空間分布，由此可以看出，東南地區(qū)是土壤濕度的高值區(qū)，土壤體積含水量都在0.25～0.40 m3/m3之間，而河套地區(qū)是土壤濕度的低值區(qū)，土壤體積含水量在0.05～0.15 m3/m3之間，土壤濕度從西北地區(qū)向東南和西南地區(qū)遞增。從0～10 cm土壤濕度的空間分布來看，華北地區(qū)較為干燥，而東南和西南地區(qū)較為濕潤，這樣的土壤濕度分布特點與孫丞虎等的研究結(jié)果[36-37]較為一致。從圖3-b可以看出，Noah-MP模式對中國東部土壤濕度的模擬結(jié)果也呈現(xiàn)出從西北地區(qū)向東南和西南地區(qū)遞增的特點，但略有差異，主要體現(xiàn)在對于新疆南部地區(qū)土壤濕度出現(xiàn)了明顯的低估，這可能是由2個原因造成的：(1)觀測站點所在地與整體環(huán)境出現(xiàn)了差異；(2)剔除了10 cm土壤溫度<0 ℃時的0～10 cm土壤濕度觀測值，造成一些土壤濕度觀測低值沒有通過質(zhì)量控制流程。Noah-MP模式在河套地區(qū)出現(xiàn)了輕微的高估；而對于東南地區(qū)和西南地區(qū)的土壤濕度模擬得很好。從空間相關(guān)系數(shù)來看，Noah-MP土壤濕度模擬值與觀測值之間的相關(guān)系數(shù)為0.538，可以通過0.05水平的顯著性檢驗。

圖4給出了Noah-MP模式0～10 cm土壤濕度模擬值與觀測值偏差的空間分布，可以看出東南地區(qū)主要出現(xiàn)低估，呈負(fù)偏差，而在北方地區(qū)主要呈正偏差；Noah-MP在新疆地區(qū)存在一個負(fù)偏差高值中心，這可能與使用的觀測值有關(guān)；同時在四川省存在一個正偏差高值中心，這可能是由于觀測儀器的參數(shù)出現(xiàn)了問題，陳東東等研究表明，四川省的人工觀測土壤濕度和自動觀測土壤濕度之間存在較大的差異，而一般認(rèn)為通過人工烘干稱質(zhì)量觀測的土壤濕度觀測值較為準(zhǔn)確，因此四川省的土壤濕度自動觀測值是否準(zhǔn)確是值得商榷的[38]。從Noah-MP模式0～10 cm土壤濕度模擬值與觀測值偏差的統(tǒng)計情況(圖5)可以看出，偏差主要分布在 -0.09～0.09 m3/m3以內(nèi)，其中Noah-MP模式0～10 cm土壤濕度模擬值在69.1%的站點與觀測值的偏差在-0.06～0.06 m3/m3之間，在37.7%的站點與觀測值的偏差在 -0.03～0.03 m3/m3之間，在呈現(xiàn)負(fù)偏差的站點(占總站點數(shù)目的52.6%)略多于呈現(xiàn)正偏差的站點(占總站點數(shù)目的47.4%)。

圖6給出了Noah-MP模式0～10 cm土壤濕度模擬值與觀測值相關(guān)系數(shù)的空間分布，可以看出土壤濕度模擬值與觀測值的相關(guān)性在東部地區(qū)較強，在大部分站點與觀測值的相關(guān)系數(shù)大于0.5；而相關(guān)性在西部地區(qū)較弱，特別是在新疆地區(qū)，模擬值在大部分站點與觀測值的相關(guān)系數(shù)均小于0.4。從Noah-MP模式0～10 cm土壤濕度模擬值與觀測值相關(guān)系數(shù)的分布情況(圖7)可以看出，相關(guān)系數(shù)在各個區(qū)間均有分布，但是主要分布在0.5及以上，其中在60%以上的站點，土壤濕度模擬值與觀測值的相關(guān)系數(shù)大于0.6，可以通過 0.05水平的顯著性檢驗。

2.2 土壤濕度模擬誤差的時間分布

在分析了站點尺度土壤濕度比較結(jié)果的基礎(chǔ)上，在區(qū)域尺度上對Noah-MP模式的土壤濕度模擬值與觀測值進行對比分析。由于在西北西部和青藏高原研究區(qū)的土壤濕度觀測站較為稀疏，代表性不足，本試驗并沒有對西北西部和青藏高原研究區(qū)的土壤濕度模擬誤差時間分布進行研究。從圖8-a、圖8-b可以看出，Noah-MP對于4—5月期間的東北研究區(qū)的0～10 cm土壤濕度存在明顯的低估，這表明Noah-MP模式在凍土融化時的土壤濕度模擬方面還須要進一步改進；從圖8-e、圖8-f可以看出，在江淮研究區(qū)，Noah-MP模式可以抓住土壤濕度的變化，但是存在輕微的低估。從全國尺度看，盡管土壤濕度模擬值與觀測值的變化趨勢基本一致，但是Noah-MP的土壤濕度模擬值在春季和秋季存在輕微的低估，而在夏季與觀測更為接近。

由表1可以看出，東北地區(qū)Noah-MP模式的0～10 cm土壤濕度模擬值與觀測值的相關(guān)性低于其他5個研究區(qū)，相關(guān)系數(shù)為0.665，而其他5個研究區(qū)的相關(guān)系數(shù)都大于0.9，從全國尺度來看，Noah-MP模式的0～10 cm土壤濕度模擬值與觀測值的相關(guān)性很強，相關(guān)系數(shù)達(dá)到0.939；從區(qū)域尺度和全國尺度來看，0～10 cm 土壤濕度模擬值與觀測值的偏差較小，都在 -0.03～0.03 m3/m3之間，其中江淮地區(qū)、西南地區(qū)的偏差較大，分別為-0.027、0.023 m3/m3；6個研究區(qū)0～10 cm土壤濕度模擬值與觀測值的均方根誤差都在-0.03～0.03 m3/m3以內(nèi)，其中江淮地區(qū)的均方根誤差最大，為 0.03 m3/m3，西北東部地區(qū)的均方根誤差最小，為 0.01 m3/m3。

3 結(jié)論與討論

基于國家氣象信息中心提供的CLDAS2.0驅(qū)動數(shù)據(jù)，使用Noah-MP陸面模式對中國區(qū)域的陸面過程進行模擬，提取土壤濕度模擬結(jié)果，并基于土壤水分自動站的逐小時土壤濕度觀測值對模擬結(jié)果進行評估與分析，結(jié)果表明：(1)Noah-MP模式能夠很好地模擬出中國區(qū)域的0～10 cm土壤濕度空間分布，偏差主要分布在-0.06～0.06 m3/m3以內(nèi)，相關(guān)系數(shù)主要分布在0.5以上；從區(qū)域尺度看，Noah-MP能夠很好地模擬出各區(qū)域的土壤濕度變化，各研究區(qū)的偏差都在-0.03～0.03 m3/m3以內(nèi)。(2)對于春季北方地區(qū)的0～10 cm 土壤濕度，Noah-MP模式的模擬結(jié)果存在明顯的低估現(xiàn)象，這可能是因為該模式在凍土融化時的土壤濕度模擬方面還存在一定的缺陷，須要進一步改進。

[1]Li H，Robock A，Liu S，et al. Evaluation of reanalysis soil moisture simulations using updated Chinese soil moisture observations[J]. Journal of Hydrometeorology，2009，6(2)：180-193.

[2]Xia Y L，Sheffield J，Ek M B，et al. Evaluation of multi-model simulated soil moisture in NLDAS-2[J]. Journal of Hydrology，2014，512(4)：107-125.

[3]尹 楠，姜琦剛，孟治國，等. 基于RADARSAT-2全極化數(shù)據(jù)反演周期性地表土壤濕度[J]. 農(nóng)業(yè)工程學(xué)報，2013(17)：72-79.

[4]李明星，馬柱國. 基于模擬土壤濕度的中國干旱檢測及多時間尺度特征[J]. 中國科學(xué)(地球科學(xué))，2015(7)：994-1010.

[5]王素萍，張存杰，宋連春，等. 多尺度氣象干旱與土壤相對濕度的關(guān)系研究[J]. 冰川凍土，2013，35(4)：865-873.

[6]Yeh T，Wetherald R T，Manabe S. The effect of soil moisture on the short-term climate and hydrology change—a numerical experiment[J]. Monthly Weather Review，1972，112(3)：474.

[7]馬柱國，符淙斌，謝 力，等. 土壤濕度和氣候變化關(guān)系研究中的某些問題[J]. 地球科學(xué)進展，2001，16(4)：563-566.

[8]朱 晨，師春香，席 琳，等. 中國區(qū)域不同深度土壤濕度模擬和評估[J]. 氣象科技，2013，41(3)：529-536.

[9]朱 智，師春香. 中國氣象局陸面同化系統(tǒng)和全球陸面同化系統(tǒng)對中國區(qū)域土壤濕度的模擬與評估[J]. 科學(xué)技術(shù)與工程，2014，14(32)：138-144.

[10]林朝暉，劉輝志，謝正輝，等. 陸面水文過程研究進展[J]. 大氣科學(xué)，2008，32(4)：935-949.

[11]Wang A H，Zeng X B. Sensitivities of terrestrial water cycle simulations to the variations of precipitation and air temperature in China[J]. Journal of Geophysical Research，2011，116(6)：2166-2181.

[12]Shen Y，Yang K，Qin J，et al. Evaluation of AMSR-E retrievals and GLDAS simulations against observations of a soil moisture network on the central Tibetan Plateau[J]. Journal of Geophysical Research：Atmospheres，2013，118(10)：4466-4475.

[13]劉樹華，李新榮，劉立超，等. 陸面過程參數(shù)化模式的研究[J]. 中國沙漠，2001，21(3)：303-311.

[14]黃春林，李 新. 陸面數(shù)據(jù)同化系統(tǒng)的研究綜述[J]. 遙感技術(shù)與應(yīng)用，2004，19(5)：424-430.

[15]劉金婷. 多陸面模式、多驅(qū)動場對新疆地區(qū)陸面過程模擬研究及結(jié)果集成[D]. 青島：中國海洋大學(xué)，2010.

[16]Maurer E P，Wood A W，Adam J C，et al. A long-term hydrologically based dataset of land surface fluxes and states for the conterminous United States[J]. Journal of Climate，2002，15(22)：3237-3251.

[17]Berg A A，F(xiàn)amiglietti J S，Walker J P，et al. Impact of bias correction to reanalysis products on simulations of North American soil moisture and hydrological fluxes[J]. Journal of Geophysical Research Atmospheres，2003，1866(D16)：1211-1222.

[18]Fekete B M，V?r?smarty C J，Roads J O，et al. Uncertainties in precipitation and their impacts on runoff estimates[J]. Journal of Climate，2004，17(2)：294-304.

[19]Sheffield J，Ziegler A D，Wood E F，et al. Correction of the high-latitude rain day anomaly in the NCEP-NCAR reanalysis for land surface hydrological modeling[J]. Journal of Climate，2004，17(19)：3814-3828.

[20]Sheffield J，Goteti G，Wood E F. Development of a 50-year high-resolution global dataset of meteorological forcings for land surface modeling[J]. Journal of Climate，2006，19(13)：3088-3111.

[21]Qian T，Dai A，Trenberth K E，et al. Simulation of global land surface conditions from 1948 to 2004.Part Ⅰ：forcing data and evaluations[J]. Journal of Hydrometeorology，2006，7(5)：953-975.

[22]Rodell M，Houser P R，Jambor U，et al. The global land data assimilation system[J]. Bulletin of the American Meteorological Society，2004，85(3)：381-394.

[23]Robock A，Luo L，Wood E F，et al. Evaluation of the North American land data assimilation system over the southern Great Plains during the warm season[J]. Journal of Geophysical Research：Atmospheres，2003，108(D22)：239-244.

[24]Xia Y L，Mitchell K，Ek M，et al. Continental-scale water and energy flux analysis and validation for North American Land Data Assimilation System project phase 2 (NLDAS-2)：2. Validation of model-simulated streamflow[J]. Journal of Geophysical Research，2012，117(D3)：812-819.

表1 6個研究區(qū)及全國0～10 cm土壤濕度Noah-MP模擬值與實測的相關(guān)系數(shù)、偏差和均方根誤差

[25] Xia Y L，Ek M，Wei H L，et al. Comparative analysis of relationships between NLDAS-2 forcings and model outputs[J]. Hydrological Processes，2012，26(3)：467-474.

[26]師春香，謝正輝，錢 輝，等. 基于衛(wèi)星遙感資料的中國區(qū)域土壤濕度EnKF數(shù)據(jù)同化[J]. 中國科學(xué)(地球科學(xué))，2011(3)：375-385.

[27]師春香. 基于EnKF算法的衛(wèi)星遙感土壤濕度同化研究[D]. 北京：中國科學(xué)院研究生院，2008.

[28]何 杰. 中國區(qū)域高時空分辨率地面氣象要素數(shù)據(jù)集的建立[D]. 北京：中國科學(xué)院研究生院，2010.

[29]劉 波，馬柱國，馮錦明. 1960—2004年新疆地區(qū)地表水熱過程的數(shù)值模擬研究Ⅰ.以觀測資料為基礎(chǔ)的陸面過程模型大氣驅(qū)動場的發(fā)展[J]. 中國沙漠，2012，32(2)：491-502.

[30]左志燕，張人禾. 中國東部春季土壤濕度的時空變化特征[J]. 中國科學(xué)(地球科學(xué))，2008(11)：1428-1437.

[31]張志富. 自動站土壤水分資料質(zhì)量控制方案的研制[J]. 干旱區(qū)地理，2013，36(1)：101-108.

[32]黃飛龍，李昕娣，黃宏智，等. 基于FDR的土壤水分探測系統(tǒng)與應(yīng)用[J]. 氣象，2012(6)：764-768.

[33]韓 帥. 基于CLDAS驅(qū)動數(shù)據(jù)的CLM3.5和SSIB2陸面模式模擬評估及干旱監(jiān)測應(yīng)用[D]. 南京：南京信息工程大學(xué)，2015.

[34]Yang Z L，Niu G Y，Mitchell K E，et al. The community Noah land surface model with multiparameterization options (Noah-MP)：2. Evaluation over global river basins[J]. Journal of Geophysical Research-Atmospheres，2011，116(D12)：1248-1256.

[35]Cai X T，Yang Z L，Xia Y L，et al. Assessment of simulated water balance from Noah，Noah-MP，CLM，and VIC over CONUS using the NLDAS test bed[J]. Journal of Geophysical Research，2014，119(24)：13751-13770.

[36]孫丞虎，李維京，張祖強，等. 淮河流域土壤濕度異常的時空分布特征及其與氣候異常關(guān)系的初步研究[J]. 應(yīng)用氣象學(xué)報，2005，16(2)：129-138.

[37]張文君，周天軍，宇如聰. 中國土壤濕度的分布與變化Ⅰ.多種資料間的比較[J]. 大氣科學(xué)，2008，32(3)：581-597.

[38]陳東東，王明田，張玉芳，等. 四川省土壤濕度自動站和人工觀測數(shù)據(jù)對比分析[J]. 安徽農(nóng)業(yè)科學(xué)，2011，39(29)：18066-18068，18147.