贛江流域氣候和土地利用變化對(duì)藍(lán)綠水的影響

李文婷，楊肖麗，任立良

(1.河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210098；2.河海大學(xué)水文水資源學(xué)院，江蘇 南京 210098)

氣候變化對(duì)區(qū)域水循環(huán)過(guò)程的影響已成為當(dāng)前研究熱點(diǎn)。氣候變化改變區(qū)域降水、氣溫以及蒸散發(fā)條件，進(jìn)而影響流域徑流過(guò)程，加劇水資源短缺問(wèn)題[1-2]。與此同時(shí)，隨著人口持續(xù)增長(zhǎng)，高強(qiáng)度的人類活動(dòng)改變了土地覆被條件，導(dǎo)致入滲、蒸發(fā)和地表徑流等一系列水文過(guò)程發(fā)生改變[3-4]。受氣候變化和土地利用變化的交互影響，地面上的水資源量被重新分配，進(jìn)一步加劇了水資源供需之間的矛盾[5]。因此，定量分析氣候變化和土地利用變化對(duì)水資源的影響，對(duì)于制定科學(xué)的水資源管理體系和提高水資源利用效率具有重要意義。

Falkenmark[6]最早提出了藍(lán)水和綠水的概念，為水資源研究與管理提供了新視角。藍(lán)水主要是指以地表徑流、土壤中流、地下徑流3種形式存在的水，綠水是指土壤水和實(shí)際蒸散量。綠水作為水分消耗的主體，其80%的水資源用于全球農(nóng)業(yè)生產(chǎn)[7]，對(duì)維持生態(tài)系統(tǒng)的穩(wěn)定起著不可替代的作用。由于藍(lán)水資源與人類發(fā)展需求(飲水、灌溉、航運(yùn)、發(fā)電、工業(yè)生產(chǎn)等)息息相關(guān)[8-9]，導(dǎo)致傳統(tǒng)的水資源研究大多集中于藍(lán)水資源，忽視了綠水資源在生態(tài)系統(tǒng)中的重要作用。研究表明，綠水資源作為植物生長(zhǎng)的基礎(chǔ)，對(duì)維持陸地生態(tài)系統(tǒng)的穩(wěn)定具有不可替代的作用[10]。藍(lán)綠水的變化主要受氣候變化和人類活動(dòng)兩方面因素的影響，在長(zhǎng)時(shí)間尺度上，氣候變化對(duì)流域水文循環(huán)過(guò)程的影響起主導(dǎo)作用；在短時(shí)間尺度上，土地利用變化則是影響流域水文要素變化的關(guān)鍵因素之一[11]。近年來(lái)，國(guó)內(nèi)外諸多學(xué)者陸續(xù)開展了流域內(nèi)藍(lán)綠水對(duì)氣候變化和土地利用變化的響應(yīng)研究。Zhang等[8]基于SWAT(soil and water assessment tool)模型，通過(guò)設(shè)置氣候變化與土地利用變化多種情景，定量研究了贛江上游多情景下藍(lán)綠水的變化情況；Huang[12]研究了氣候變化和土地利用變化交互影響下的全球農(nóng)業(yè)藍(lán)綠水資源的變化情況。Akbar等[13]研究了氣候和土地利用變化對(duì)巴基斯坦昆哈爾河流域徑流的影響。趙安周等[14]借助SWAT模型探討了氣候變化和人類活動(dòng)對(duì)渭河流域藍(lán)綠水資源的影響。

贛江是長(zhǎng)江的重要支流，地處亞熱帶濕潤(rùn)地區(qū)。贛江流域內(nèi)氣候變化劇烈，降水年內(nèi)、年際分配不均，導(dǎo)致了洪澇、干旱等氣象災(zāi)害頻發(fā)[15]。此外，由于流域內(nèi)經(jīng)濟(jì)快速發(fā)展、城市化進(jìn)程加快導(dǎo)致建筑用地?cái)U(kuò)張，下墊面條件發(fā)生改變。基于此，本文利用標(biāo)準(zhǔn)化降水指數(shù)(standardized precipitation index, SPI)分析典型年份贛江流域藍(lán)綠水資源的時(shí)空分布，通過(guò)設(shè)置不同情景來(lái)定量分析氣候變化和土地利用變化對(duì)贛江流域藍(lán)綠水的影響，以期為贛江流域水資源管理提供參考。

1 流域概況

贛江流域位于長(zhǎng)江中下游南岸，地處24 °31′N～28 °45′N、113 °34′E～116 °38′E之間，總面積8.35 萬(wàn)km2，是鄱陽(yáng)湖流域最大的子流域(圖1)。贛江上游流域主要為贛州市，中游流域主要為吉安市和撫州市，下游流域包括新余及宜春市的各市縣。贛江流域氣候?qū)儆趤啛釒Ъ撅L(fēng)氣候，氣候溫和，雨量充沛，多年平均降水量約為1 600 mm，多年平均氣溫約為18 ℃。贛江流域地勢(shì)大致呈南高北低的趨勢(shì)，西北部為九嶺山，東部與武夷山接壤，北部多為低丘崗地和平原。流域內(nèi)土地利用類型以林地和耕地為主，其中林地占比超過(guò)65%。流域土壤種類多樣，涵蓋紅壤、黃壤等20多種。

圖1 贛江流域概況

2 數(shù)據(jù)來(lái)源與研究方法

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

SWAT模型需要詳細(xì)的流域基礎(chǔ)數(shù)據(jù)，表1為構(gòu)建贛江流域SWAT模型所需的各類數(shù)據(jù)，主要包括DEM數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)、氣象數(shù)據(jù)、水文數(shù)據(jù)等。

表1 SWAT模型基礎(chǔ)數(shù)據(jù)

2.2 SPI

SPI用于表征某時(shí)間段內(nèi)降水量出現(xiàn)的概率，該指數(shù)未考慮具體的干旱機(jī)理，是一種計(jì)算簡(jiǎn)單、穩(wěn)定性較好的干旱指數(shù)，廣泛應(yīng)用于不同氣候區(qū)的干旱識(shí)別[16]。SPI具體計(jì)算過(guò)程見文獻(xiàn)[17]。干旱等級(jí)根據(jù)GB/T 20481—2017《氣象干旱等級(jí)》進(jìn)行劃分(表2)。本文采用Matlab程序方法，計(jì)算12月尺度的SPI值(SPI-12)，確定贛江流域1965—2016年的氣候典型年份。

表2 SPI的旱澇等級(jí)

2.3 SWAT模型

選用SWAT模型對(duì)區(qū)域藍(lán)綠水進(jìn)行模擬。SWAT模型是基于物理機(jī)制的大尺度半分布式水文模型，它可以結(jié)合流域內(nèi)多種地理要素來(lái)模擬復(fù)雜的流域水文過(guò)程[18]。SWAT模型能夠輸出藍(lán)綠水的各個(gè)組成分量，廣泛應(yīng)用于藍(lán)綠水模擬[19]。根據(jù)藍(lán)綠水的定義以及SWAT模型輸出結(jié)果可計(jì)算藍(lán)水量、綠水量以及綠水系數(shù)，計(jì)算過(guò)程如下：

WB=WYLD+WA_RCHG

(1)

WG=ET+WS

(2)

(3)

式中：WB為藍(lán)水量，mm；WYLD為子流域產(chǎn)水量，mm；WA_RCHG為深層含水層補(bǔ)給量，mm；WG為綠水量，mm；ET為實(shí)際蒸散發(fā)量，mm；WS為土壤含水量，mm；CGW為綠水系數(shù)，是指流域綠水量所占藍(lán)綠水總量的比例。

2.4 累積距平

累積距平是一種常用的判斷趨勢(shì)變化程度的方法，它是由曲線直觀判斷某一要素的變化趨勢(shì)[20]。對(duì)于序列xi，其某一時(shí)刻s的累積距平為

(4)

其中

2.5 土地利用動(dòng)態(tài)度

土地利用動(dòng)態(tài)度K用來(lái)表示某一時(shí)間段內(nèi)土地利用類型的變化趨勢(shì)和速度[21]，計(jì)算公式如下：

(5)

式中：Sa為研究初期的土地利用面積，km2；Sb為研究末期的土地利用面積，km2；T為時(shí)間長(zhǎng)度。

3 結(jié)果與分析

3.1 氣象要素與土地利用變化

3.1.1氣象要素變化

氣候變化是影響流域內(nèi)藍(lán)綠水分布的重要因素之一，分析贛江流域1965—2016年降水量和氣溫的變化情況(圖2)，有助于理解藍(lán)綠水的變化特征。贛江流域多年平均降水量達(dá)1 636 mm，在2002年達(dá)到最高值(2 124 mm)；多年平均氣溫為17.96 ℃；多年平均最高氣溫為22.75 ℃；多年平均最低氣溫為14.52 ℃，三者均在2015年達(dá)到最大值(19.61 ℃、24.10 ℃和16.52 ℃)。利用非參數(shù)統(tǒng)計(jì)檢驗(yàn)對(duì)流域內(nèi)降水量、氣溫進(jìn)行趨勢(shì)性檢驗(yàn)，結(jié)果表明，多年平均降水量統(tǒng)計(jì)值Z為0.120 7，未通過(guò)95%的顯著性檢驗(yàn)，表現(xiàn)為不顯著的增大；多年平均氣溫、多年平均最高氣溫及多年平均最低氣溫的統(tǒng)計(jì)值Z分別為0.478 1、0.340 9和0.567 1，均通過(guò)顯著性檢驗(yàn)，呈現(xiàn)顯著的增大趨勢(shì)。

圖2 贛江流域1965—2016年降水量和氣溫變化

利用累積距平法對(duì)贛江流域1965—2016年的降水量、多年平均氣溫、平均最高氣溫和平均最低氣溫進(jìn)行突變點(diǎn)檢驗(yàn)分析(圖3)。流域內(nèi)降水量累積距平過(guò)程線波動(dòng)大，震蕩頻率較高，年際變化較大，1991年是降水量波動(dòng)上升的起始點(diǎn)；平均氣溫、平均最高氣溫和平均最低氣溫大致分為兩個(gè)變化階段，1965—1996年總體呈減小趨勢(shì)，1997—2016年總體呈增大趨勢(shì)，1997年是平均氣溫和平均最高氣溫的突變點(diǎn)，1996年是平均最低氣溫的一個(gè)突變點(diǎn)。

(a)降水量

基于贛江流域內(nèi)12個(gè)氣象站點(diǎn)的降水量、平均氣溫、平均最高氣溫、平均最低氣溫的年均值數(shù)據(jù)，采用反距離權(quán)重法進(jìn)行空間插值，得到各氣象要素的空間分布圖(圖4)。贛江流域內(nèi)降水空間分布不均，流域的東部、西部和西北部地區(qū)降水較為豐富，流域的南部和中部地區(qū)降水較少。流域內(nèi)平均氣溫、平均最高氣溫和平均最低氣溫三者的空間變化趨勢(shì)基本保持一致，均呈現(xiàn)東西分布不均，流域西部氣溫較低，東南部氣溫較高，溫差較大。

(a)降水量(單位：mm)

3.1.2土地利用變化

在贛江流域1990年和2015年兩期土地利用中(表3)，林地和耕地是主要土地利用類型，兩者占流域總面積90%以上，其中林地占比超過(guò)65%，草地、水域、建筑用地和裸地占比較少，均不超過(guò)6%。從變化趨勢(shì)來(lái)看，1990—2015年間贛江流域內(nèi)建筑用地單一動(dòng)態(tài)度為1.38%，其年增速最大，這與流域內(nèi)經(jīng)濟(jì)快速發(fā)展、城市化進(jìn)程加快密切相關(guān)。水域單一動(dòng)態(tài)度為0.11%，呈現(xiàn)增大的趨勢(shì)，這與流域內(nèi)水利設(shè)施的修建、退耕還湖相關(guān)。耕地、裸地、林地和草地單一動(dòng)態(tài)度分別為-0.06%、-0.39%、-0.01%和-0.27%，總體均呈現(xiàn)減小的趨勢(shì)。

表3 1990—2015年贛江流域土地利用類型變化

3.2 SWAT模型建立與評(píng)價(jià)

3.2.1SWAT模型建立

首先根據(jù)基礎(chǔ)數(shù)據(jù)建立土壤數(shù)據(jù)庫(kù)、氣象數(shù)據(jù)庫(kù)等，然后根據(jù)贛江流域DEM分辨率，通過(guò)設(shè)置集水區(qū)閾值和流域總出口點(diǎn)，經(jīng)模型運(yùn)算將贛江流域劃分為63個(gè)子流域。在生成的子流域基礎(chǔ)上，綜合考慮各子流域的土地利用、坡度、土壤類型分布情況，生成了941個(gè)水文響應(yīng)單元(HRU)。

3.2.2參數(shù)敏感性分析

SWAT模型參數(shù)眾多，根據(jù)SWAT-CUP中的全局敏感性對(duì)參數(shù)進(jìn)行分析，其中t值表示參數(shù)敏感性程度，t的絕對(duì)值越大表示參數(shù)越敏感；p值表示參數(shù)敏感性的顯著性，p越接近0表示越顯著。根據(jù)敏感性分析結(jié)果，選取10個(gè)敏感性較高的參數(shù)作為贛江流域的率定參數(shù)(表4)。其中，對(duì)徑流產(chǎn)生較為敏感的參數(shù)有CH_K2(主河道水力傳導(dǎo)率)、ALPHA_BF(基流α因子)、CN2(SCS徑流曲線數(shù))。參數(shù)CH_K2反映河道匯流過(guò)程；參數(shù)ALPHA_BF主要影響地下水過(guò)程；參數(shù)CN2反映降水前流域下墊面特征，與土壤濕度有關(guān)，土壤濕度愈大敏感性愈強(qiáng)。在模型運(yùn)行過(guò)程中，利用贛江流域外洲水文站1965—2016年逐月徑流觀測(cè)數(shù)據(jù)，結(jié)合2015年土地利用數(shù)據(jù)，設(shè)置1960—1964年為模型預(yù)熱期，1965—1990年為率定期，1991—2016年為驗(yàn)證期。參數(shù)率定采用SWAT-CUP中的SUFI-2算法，外洲水文站敏感參數(shù)最優(yōu)值如表4所示。

表4 SWAT模型參數(shù)選擇及率定結(jié)果

3.2.3模型結(jié)果評(píng)價(jià)

本文選用決定系數(shù)R2、Nash-Suttcliffe系數(shù)(Ens)、相對(duì)誤差RE共3個(gè)指標(biāo)[22]對(duì)SWAT模型模擬結(jié)果進(jìn)行評(píng)價(jià)(圖5)。結(jié)果表明，率定期評(píng)價(jià)結(jié)果為：R2= 0.94，Ens= 0.93，RE=-4.1%；驗(yàn)證期評(píng)價(jià)結(jié)果為：R2= 0.93，Ens= 0.90，RE= 9.5%。總體模型模擬結(jié)果較好，說(shuō)明SWAT模型可用于研究贛江流域藍(lán)綠水對(duì)氣候變化和人類活動(dòng)的響應(yīng)模擬研究。

(a)率定期

3.3 典型干濕年份確定

選用SPI-12來(lái)確定贛江流域1965—2016年的典型年份。由圖6確定1972年、1998年和2014年分別為贛江流域的干旱年、濕潤(rùn)年和正常年。

圖6 贛江流域歷年SPI-12值

3.4 典型年份下的藍(lán)綠水時(shí)空變化特征

贛江流域在干旱年、正常年、濕潤(rùn)年的藍(lán)水量分別為638.99 mm、832.81 mm和1 283.14 mm。濕潤(rùn)年的藍(lán)水量是干旱年的2.01倍，是正常年的1.54倍，存在較大差異，這主要是由于濕潤(rùn)年的降水量明顯高于干旱年和正常年。與藍(lán)水量相比，贛江流域的綠水量在3個(gè)典型年份的變化比較穩(wěn)定，差異不大。干旱年、正常年、濕潤(rùn)年的綠水量分別為833.52 mm、830.95 mm、815.97 mm。此外，贛江流域干旱年、正常年、濕潤(rùn)年的綠水系數(shù)也存在一定差異，分別為57.26%、51.29%和39.64%。干旱年的綠水系數(shù)是濕潤(rùn)年的1.44倍，表明干旱年的綠水資源比例明顯高于濕潤(rùn)年，對(duì)于干旱年份，綠水資源對(duì)維持生態(tài)平衡尤為重要。

贛江流域典型年份的藍(lán)水量、綠水量和綠水系數(shù)空間分布見圖7。總體來(lái)看，3個(gè)典型年份的藍(lán)水量空間分布基本保持一致，均呈現(xiàn)自北向南減少的趨勢(shì)。在濕潤(rùn)年，藍(lán)水量最為豐富，流域各地藍(lán)水量均大于550 mm，北部和東南部大部分地區(qū)超過(guò)940 mm；在正常年，流域內(nèi)藍(lán)水量位于540～1 300 mm之間，南部地區(qū)藍(lán)水量不足600 mm；在干旱年，藍(lán)水量明顯減少，整體位于400～900 mm之間，南部地區(qū)藍(lán)水量不足500 mm。

(a)藍(lán)水量(單位：mm)

流域內(nèi)3個(gè)典型年份的綠水量空間分布總體一致，且差異較小，均呈現(xiàn)從東北向西南遞減的趨勢(shì)。在濕潤(rùn)年，大部分地區(qū)的綠水量在800 mm以下，東北部少數(shù)地區(qū)大于900 mm；在正常年，流域內(nèi)大部分地區(qū)綠水量為770～880 mm，北部綠水量較大，大于910 mm；在干旱年，綠水量較為豐富，流域大部分地區(qū)處于810 mm以上，北部地區(qū)綠水量為870～960 mm。

綠水系數(shù)的空間分布總體表現(xiàn)為由南部地區(qū)向北部地區(qū)遞減。在濕潤(rùn)年，流域內(nèi)大部分地區(qū)的綠水系數(shù)在39%以下，西南部少數(shù)地區(qū)在47%以上；在正常年，南部大多數(shù)地區(qū)綠水系數(shù)在49%以上，北部地區(qū)在43%以下；在干旱年，流域內(nèi)綠水系數(shù)與其他典型年份相比較高，大部分地區(qū)在54%以上，這主要是由于干旱年份降水量少、氣溫較高，流域內(nèi)蒸散發(fā)作用消耗的水資源量高于濕潤(rùn)年和正常年，從而導(dǎo)致干旱年份的綠水系數(shù)高于其他典型年份。

3.5 氣候與土地利用變化對(duì)藍(lán)綠水的時(shí)空影響

3.5.1情景設(shè)置

采用情景對(duì)比法來(lái)定量分析氣候變化和土地利用變化對(duì)贛江流域藍(lán)綠水的影響。贛江流域1965—2016年降水量突變年份為1991年，同時(shí)保證劃分的兩段氣象數(shù)據(jù)時(shí)間長(zhǎng)度保持一致，因此選擇1990年作為劃分氣象數(shù)據(jù)的時(shí)間點(diǎn)。具體設(shè)置見表5。

3.5.2藍(lán)綠水時(shí)空變化特征

藍(lán)綠水的變化主要受氣候變化和人類活動(dòng)兩方面因素的影響，土地利用變化是人類活動(dòng)作用的直接表現(xiàn)。因此，為定量分析氣候變化和土地利用變化對(duì)藍(lán)綠水的影響，本文將不同氣候情景、土地利用情景分別代入校準(zhǔn)好的SWAT模型，各情景下藍(lán)綠水的變化情況見表5。

表5 贛江流域氣候與土地利用情景設(shè)置

在土地利用變化情景下，藍(lán)水量、綠水量、綠水系數(shù)分別為843.22 mm、821.11 mm和49.75%，與基準(zhǔn)期相比，變化微弱，這主要是因?yàn)?990—2015年贛江流域土地利用類型未發(fā)生較大變化，僅有建筑用地的單一動(dòng)態(tài)度為1.38%，變化程度較大，其余土地利用類型的變化量均較小。因此，土地利用變化對(duì)流域內(nèi)水文過(guò)程影響作用較小。

在氣候變化情景下，藍(lán)水量、綠水量、綠水系數(shù)分別為919.23 mm、851.25 mm和48.55%，與基準(zhǔn)期相比，藍(lán)水量、綠水量分別增加了75.52 mm和30.65 mm，綠水系數(shù)減小了2.33%。這主要是由于流域內(nèi)1965—2016年降水量、氣溫均呈上升趨勢(shì)，其中降水量的增加會(huì)直接導(dǎo)致藍(lán)水量的增加，而氣溫的升高首先會(huì)增加流域內(nèi)的蒸散發(fā)量，其次蒸散發(fā)量作為綠水的主要組成部分，其增加將進(jìn)一步增加綠水量。綠水系數(shù)的減小可能是由于降水轉(zhuǎn)化的藍(lán)水量較多，使得綠水資源在藍(lán)綠水總量中所占比例有所下降。

在共同變化情景下，藍(lán)水量、綠水量和綠水系數(shù)分別為918.31 mm、852.06 mm和48.61%。共同變化情景下的水資源量變化趨勢(shì)與氣候變化情景下基本一致，表明氣候變化對(duì)贛江流域藍(lán)綠水量影響較大。

贛江流域各情景下的藍(lán)水量、綠水量和綠水系數(shù)空間分布見圖8。整體而言，4種情景下的藍(lán)水量基本上呈現(xiàn)從北向南減少的趨勢(shì)。土地利用變化對(duì)藍(lán)水量的影響不明顯；氣候變化情景下，流域內(nèi)藍(lán)水量整體呈增加趨勢(shì)，主要集中于流域西北部和東南部，這與流域內(nèi)降水的空間分布有一定關(guān)系，贛江流域降水分布由北向南呈現(xiàn)遞減的趨勢(shì)，降水量的分布直接影響藍(lán)水量的分布；共同變化情景下的藍(lán)水量空間變化與氣候情景下基本一致，總體呈增加趨勢(shì)。

(a)藍(lán)水量(單位：mm)

各情景下的綠水量的空間分布總體呈現(xiàn)由東向西遞減的趨勢(shì)。與藍(lán)水量變化趨勢(shì)基本相似，流域內(nèi)綠水量在土地利用變化情景下變化微弱，在氣候變化、共同變化情景下均呈增加趨勢(shì)。綠水量的空間分布一定程度上受流域內(nèi)氣溫分布的影響，贛江流域西部氣溫較低，東南部氣溫較高，進(jìn)而使綠水量具有相似的空間分布特征。4種情景下綠水系數(shù)的空間分布均呈現(xiàn)由西南向東北地區(qū)遞減的趨勢(shì)。土地利用情景下綠水系數(shù)變化不大，氣候變化情景下流域西南部的綠水系數(shù)呈減小趨勢(shì)，共同變化情景下綠水系數(shù)變化趨勢(shì)與氣候變化情景下基本一致。綠水系數(shù)的減小可能是由于在氣候變化情景下藍(lán)水量的增幅大于綠水量，綠水量占藍(lán)綠水總量的比例出現(xiàn)了一定程度的下降。

4 結(jié) 論

a.SWAT模型在贛江流域的月徑流模擬精度較好，表明模型可用于贛江流域藍(lán)綠水對(duì)氣候變化和人類活動(dòng)的響應(yīng)模擬。

b.根據(jù)SPI確定1972年、1998年、2014年分別為贛江流域的干旱年、濕潤(rùn)年和正常年。贛江流域濕潤(rùn)年的藍(lán)水量是干旱年的2.01倍，是正常年的1.54倍；贛江流域綠水量在3個(gè)典型年份的變化比較穩(wěn)定；贛江流域干旱年的綠水系數(shù)是濕潤(rùn)年的1.44倍。從空間分布來(lái)看，在3個(gè)典型年份下，藍(lán)水量呈現(xiàn)自北向南減少的趨勢(shì)；綠水量呈現(xiàn)從東北向西南遞減的趨勢(shì)；綠水系數(shù)總體表現(xiàn)為由南部地區(qū)向北部地區(qū)遞減。

c.土地利用變化情景下，藍(lán)水量、綠水量和綠水系數(shù)分別為843.22 mm、821.11 mm和49.75%；氣候變化情景下，藍(lán)水量、綠水量和綠水系數(shù)分別為919.23 mm、851.25 mm和48.55%；共同變化情景下，藍(lán)水量、綠水量和綠水系數(shù)分別為918.31 mm、852.06 mm和48.61%，表明氣候變化對(duì)贛江流域藍(lán)綠水量影響較大。從空間分布來(lái)看，4種情景下的藍(lán)水量整體呈現(xiàn)從北向南減少的趨勢(shì)；綠水量表現(xiàn)為由東向西遞減；綠水系數(shù)則呈由西南向東北地區(qū)遞減的趨勢(shì)。