氣候變化影響下瀾湄流域下游水稻生產用水量模擬與分析

謝詩猛 劉登峰 劉慧 胡宏昌 董志強 王天賜 明廣輝

摘要：農業用水是社會經濟用水的主要行業，氣候變化對農業用水的影響可能改變區域的水資源供需情勢。以瀾滄江-湄公河流域下游為研究區，基于ERA5.Land數據集和最新的CMIP6氣候預測數據，選用SSP1.2.6、SSP2.4.5、SSP5.8.5三種排放情景，AquaCrop模型將非生產性土壤蒸發和生產性作物蒸騰分離，將水稻生長期內蒸騰總量作為水稻生產用水量，模擬了瀾滄江-湄公河流域下游歷史和未來時期水稻生產用水量，分析了水稻生產用水量的變化情況及其與溫度、降水和CO2濃度之間的相關關系。結果表明：瀾滄江-湄公河流域下游水稻生產用水量在空間上表現為北多南少，總體呈現逐年減少趨勢，且在SSP5.8.5情景下趨勢更加明顯。SSP5.8.5情景下的未來遠期，瀾滄江-湄公河流域下游水稻生產用水量將減少29.7%。與溫度和降水相比，水稻生產用水量和CO2濃度之間的相關性最強。SSP5.8.5情景下的未來遠期，在泰國的相關系數為-0.875，其余國家各季度下的相關系數均小于-0.9。

關鍵詞：氣候變化；作物模型；用水量；相關性；瀾湄流域下游

中圖分類號：TV93? 文獻標識碼：A? 文章編號：1001.9235（2024）01.0013.21

Simulation and Analysis of Water Consumption for Rice Production in Downstream Lancanng.Mekong River Basin under the Influence of Climate Change

XIE Shimeng1，LIU Dengfeng1*，LIU Hui2，HU Hongchang3，DONG Zhiqiang2，WANG Tianci1，MING Guanghui4

（1.State Key Laboratory of Eco.hydraulics in Northwest Arid Region of China，School of Water Resources and Hydropower，

Xi'an University of Technology，Xi'an 710048，China; 2.China Institute of Water Resources and Hydropower Research，

Beijing 100038，China; 3.Department of Hydraulic Engineering，Tsinghua University，Beijing 100084，China;

4.Key Laboratory of Water Management and Water Security for Yellow River Basin （Ministry of Water Resources），Yellow River Engineering Consulting Co.，Ltd.，Zhengzhou 450003，China）

Abstract： Agriculture water consumption is the main sector of socio.economic water consumption.The impact of climate change on agriculture water consumption may change the situation of water supply and demand in a region.The downstream Lancang.Mekong River Basin is selected as the study area.Based on the ERA5.Land dataset and the latest CMIP6 climate projection data，three emission scenarios （SSP1.2.6，SSP2.4.5，and SSP5.8.5） are selected in this study，and the AquaCrop model is used to separate non.productive soil evaporation from productive crop transpiration.The total transpiration during the rice growth period is taken as water consumption for rice production，and the historical and subsequent water consumption for rice production in the downstream Lancang.Mekong River Basin is simulated.The changes in water consumption for rice production and its correlation with temperature，precipitation，and CO2 concentration are analyzed.The results show that the water consumption for rice production in the downstream Lancang.Mekong River Basin is manifested as more water consumption in the north and less water consumption in the south and has an overall decreasing trend annually.The trend is more obvious under the SSP5.8.5 scenario.In the SSP5.8.5 scenario，water consumption for rice production will reduce by 29.7% in the downstream Lancang.Mekong Basin in the far future.Compared with temperature and precipitation，the strongest correlation is found between water consumption for rice production and CO2 concentration.In the far future of the SSP5.8.5 scenario，the correlation coefficient is -0.875 in Thailand and less than -0.9 in other countries in all seasons.

Keywords：climate change;crop model;water consumption;correlation;downstream Lancang-Mekong River Basin

農業是對氣候變化最敏感，受影響最大的產業［1-2］。隨著氣候變化的加劇，農業氣候條件正在惡化，干旱壓力增加，水資源的時空分布也發生了顯著變化，農業水資源供需矛盾愈發突出［3］。近百年來，氣候變化的主要特征是全球變暖，在過去的50 a里，溫度的上升速度已達到了0.13℃/10a，且氣溫有著持續上升的趨勢［4］。研究發現：溫度和CO2濃度升高會干擾植物的正常光合作用，全球多個研究區域的主要作物（如小麥、玉米、水稻）物候期已經開始發生改變，未來農作物生產具有更大的不確定性［5-6］。氣候變化在影響作物生長發育的同時，也導致作物的耗水過程發生了變化［7］。有的研究指出氣候變暖使水稻生產用水量增加，也有研究認為降雨增加、日照時數減少和風速降低等各方面的綜合原因，水稻生產用水量將減少［8-9］。

瀾滄江-湄公河流域（簡稱瀾湄流域）是亞洲最重要的跨境河流域，農業是該流域內國家重要的經濟命脈［10］。在氣候變化背景下，瀾湄流域農業生產面臨著巨大挑戰，特別是在人口密集的下游地區，正面臨著與日俱增的干旱與洪澇事件的威脅［11-12］。眾多學者針對流域內主要作物對氣候變化的響應和干旱等氣象災害的時空演變特征展開了廣泛深入的研究，研究表明瀾湄流域的干旱有著向下游擴散的特征，盡管水稻產量隨降水和CO2濃度的升高而增加，但在未來的降水模式下仍需擴大灌溉面積［13-14］。在未來，盡管流域內降水量呈增加趨勢，但瀾湄流域的干旱季節特征愈發明顯，雨季將更加濕潤，而旱季干旱將更為嚴重［15-16］。Yamauchi等［17］研究也表明，雨季干旱期延長和前期降雨不足將增加雨養水稻的干旱風險。

為了揭示未來氣候變化情景下瀾滄江-湄公河流域的農業灌溉用水對氣候變化的響應，本研究基于ERA5.Land數據集和最新的CMIP6氣候預測數據，選用SSP1.2.6、SSP2.4.5、SSP5.8.5三種排放情景，基于AquaCrop模型模擬了瀾湄流域下游水稻生產用水量，并利用Mann.Kendall趨勢檢驗、Pearson相關系數分析了水稻生產用水量的時空變化趨勢及其與氣象要素的相關關系。研究成果為流域內應對氣候變化、優化灌溉制度、提高灌溉效率、合理規劃水資源利用提供參考。

1 研究區概況

瀾滄江-湄公河流域的農田主要分布在流域下游，其主要農作物為水稻，水稻生產幾乎占柬埔寨、老撾、泰國和越南4個國家80%以上的農業產出［18］。本研究以瀾湄流域下游為研究區，主要包括柬埔寨、老撾、泰國和越南4個國家共76個省，其中柬埔寨20個省，老撾17個省，泰國23個省，越南16個省，總面積約64萬km2，約為瀾湄流域總面積的3/4［19-20］，因缺乏緬甸的相關數據，本研究中不考慮緬甸，研究區域見圖1。

瀾湄流域下游屬熱帶季風氣候，有明顯的雨季和旱季，雨季約為每年的5—10月，旱季約為11月至次年4月，全年降水量分布不均，雨季降水量約占全年降水量80%［21-22］。受地形、氣候的影響，降水在空間分布上也非常不均勻，降水時空分布不均勻導致瀾湄流域下游洪澇與干旱的頻繁發生［23-24］。流域內氣溫變化均勻，年平均氣溫為26℃左右，平均相對濕度在60%～80%［25］。

2 數據來源和研究方法

2.1 數據來源

收集了瀾湄流域下游4個國家的共76個省份的氣候數據、作物數據、土壤數據和灌溉數據。選用的氣候數據主要包括日最高溫度Tmax、日最低溫度Tmin、日降雨量P，日參考蒸散發Et0和年平均CO2濃度。日參考蒸散發為基于其他數據變量采用FAO推薦Penman.Monteith公式計算得出，CO2數據來自美國國家海洋和大氣管理局的全球監測實驗室（NOAA）。

歷史時期（1982—2020年）氣候數據提取自ECMWF（European Entre for Medium.Range Weather Forecasts）發布的再分析數據集ERA5.Land。

未來時期（2021—2100年）氣候數據為基于3種共享社會經濟路徑（SSP1.2.6、SSP2.4.5、SSP5.8.5）下，參加CMIP6的8個氣候模式的氣候模擬預測數據，見表1。

本研究所使用的作物數據來源于流域內各國家統計局、農業大數據中心和聯合國糧農組織（FAO，下同）；土壤數據提取自聯合國糧農組織提供的世界土壤類型數據庫；灌溉數據提取自FAO土地和水資源司聯合德國Rheinische Friedrich.Wilhems大學開發的全球灌溉渠地圖數據集Irrigation Areas v.5（Global Map of Irrigation Areas.Version 5）。

2.2 研究方法

2.2.1 AquaCrop模型

AquaCrop是FAO水土司開發的一種作物生長模型，廣泛被用于評估環境和管理對作物生產的影響。AquaCrop模擬草本作物對水的產量反應，特別適合于瀾湄流域下游水是作物生產中關鍵限制因素的情況。

模型將蒸發分為非生產性的土壤蒸發（E）和用于作物生產的作物蒸騰（Tr），而水稻種植時大部分時段地面是水田，模型只模擬生產性作物蒸騰，本研究將作物生長期內作物蒸騰總量定義為作物生長發育過程中的總用水量。模型計算作物蒸騰由作物系數KcTr乘以參考蒸散量 ，并考慮水分脅迫和冷應力脅迫得到，其計算公式如下：

Tr=（KsKsTr，xKcTr）ET0（1）

式中 Ks——水分脅迫系數；KsTr，x ——冷應力脅迫系數；KcTr——作物蒸騰系數；ET0——參考蒸散量。

AquaCrop設定了雨養和灌溉2種田間管理模式。在研究區域內，柬埔寨和泰國的水稻生產主要集中在雨季，超過80%農田僅在雨季種植水稻［26-28］；在老撾水稻種植主要為雨季（夏季）和冬春季［29］；越南年內種植雨季（夏季）、冬春季和晚季三季水稻［30］，根據各國水稻種植制度設定了不同的田間管理模式（表2）。

本研究中雨養模式不采取任何灌溉措施，而灌溉模式則設定為始終保持土壤有效水分為田間持水量的70%（AquaCrop模型推薦，這基本滿足所有作物的用水需求）。特別的是，越南的夏季稻分為雨養和灌溉兩種田間管理模式，而冬春稻、晚季稻和老撾冬春稻則僅種植在灌溉所覆蓋的地區，因此，采取全灌溉模式。

2.2.2 Penman.Monteith公式

依據近年來的研究結果，為了提高對作物生產用水量的計算精度，聯合國糧農組織（FAO）推薦使用參考蒸散發的概念，并提出了計算參考蒸散發的Penman.Monteith公式：

式中 ET0——參考蒸散量，mm；Rn——凈輻射量，MJ/（m2·d）；G——土壤熱通量； γ——濕度計常數，kPa/℃；T——空氣平均溫度，℃；u2——地面2 m高處的平均風速，m/s；es——空氣飽和水汽壓，kPa；ea——空氣實際水汽壓；Δ——飽和水汽壓隨溫度的變化率，kPa/℃。

2.2.3 模型評價指標

采用均方根誤差RMSE和相對誤差RE評價AquaCrop模型對瀾湄流域下游水稻產量的模擬結果，RMSE和RE計算公式如下：

式中 Si——模擬產量，t/hm2；Oi——觀測產量，t/hm2；i——時間序列；n——系列長度。

2.2.4 Mann.Kendall趨勢檢驗

Mann.Kendall趨勢檢驗法是一種非參數統計檢驗方法，廣泛應用于降水、蒸發、徑流等要素隨時間序列變化的情況［31］。本研究采用此方法檢驗水稻生產用水量的時間變化趨勢，所有結果在5%水平下檢驗趨勢的顯著性。

2.2.5 Pearson相關系數

Pearson相關系數被廣泛地應用于度量2個變量之間的相關程度，其值介于-1～1，由皮爾遜在19世紀80年代提出。本研究采用Pearson相關系數來評估水稻模擬產量和生長期內溫度和降水的相關關系，計算公式如下：

3 結果與分析

3.1 水稻生產用水量

采用收集到的各省實際水稻產量數據進行AquaCrop模型的率定和驗證，然后以省為單位對歷史時期和未來時期不同情景下的水稻生產用水量進行模擬，根據各國家內各省份的多年平均收獲面積加權得到各國的水稻生產用水量，多季種植的水稻則通過每季的多年平均收獲面積加權平均得到最終結果，代表平均每季水稻的生產用水量。

3.1.1 模型率定和驗證

各國平均統計產量、模擬產量、RMSE和RE見表3。在率定期間，老撾的RMSE最大值為0.342 t/hm2，泰國的RMSE最小值為0.057 t/hm2，分別占模擬平均產量的6.2%和2.4%。越南相對誤差最大，RE值為0.017，老撾的相對誤差最小，RE值為0.001。在驗證期內，泰國模擬水稻的均方根誤差和相對誤差均最小，RMSE值0.089 t/hm2，RE值為0.014，模型模擬效果最好；越南模擬水稻的均方根誤差最大，RMSE值為0.231 t/hm2，柬埔寨模擬水稻的相對誤差最大，RE值為0.044?？傮w而言，各國家在率定期和驗證期的RMSE和RE值分別小于0.342 t/hm2和0.044，模型模擬滿足精度要求。

3.1.2 歷史時期水稻用水量的模擬

在1982—2020年的歷史時期瀾湄流域下游水稻生產平均用水量為296.4 mm，其中柬埔寨、老撾、泰國和越南水稻生產平均用水量分別為302.9、347.2、308.0、258.8 mm。柬埔寨和泰國水稻生產用水量的年際波動相對更大，老撾水稻生產用水量年際變化相對平穩，而越南水稻生產用水量年際變化在初期較為平穩，后期波動逐漸增大（圖2）。水稻生產用水量在空間上總體表現為北高南低的分布特征（圖3）。

對每個省歷史時期的水稻生產模擬用水量進行Mann.Kendall趨勢檢驗（結果通過Z≥1.96顯著性水平檢驗，圖4中slope大于0表示呈上升趨勢，小于0表示呈下降趨勢），結果顯示：歷史時期內瀾湄流域下游水稻生產用水量總體表現為下降趨勢，水稻生產用水量表現出顯著下降趨勢的省份一共有55個，其中柬埔寨11個省，老撾11個省，泰國18個省，越南15個省。有1個省份（泰國Nakhon Ratchasima?。┧旧a用水在歷史時期表現出顯著的增加趨勢，其余20個省份水稻生產用水量在歷史時期均無明顯趨勢。

3.1.3 未來時期水稻生產用水量的模擬

基于未來預測數據，本研究以省為單位對瀾湄流域下游2個未來時期（未來近期：2021—2060年，未來遠期：2061—2100年）水稻生產的用水量進行了模擬，并對比了瀾湄流域下游歷史時期的水稻生產模擬用水量（1982—2020年），見圖5。

對比歷史時期的模擬水稻生產用水量，在未來不同情境下各國水稻生產用水量有所降低，其降低趨勢總體上表現為在高排放情景SSP5.8.5下更為明顯。在未來近期，SSP1.2.6情景、SSP2.4.5情景、SSP5.8.5情景下瀾湄流域下游水稻生產用水量將分別減少17.8%、21.4%、20.4%。其中，柬埔寨水稻生產用水量將分別減少約17.3%、20.5%、20%；老撾的水稻生產用水量分別減少約18.0%、19.5%、21.1%；泰國的水稻生產用水量分別減少約13.1%、18.6%、16.0%；越南的水稻生產用水量別減少27.3%、28.2%、29.4%。在未來遠期，SSP1.2.6情景、SSP2.4.5情景、SSP5.8.5情景下瀾湄流域下游水稻生產用水量將分別減少18.5%、25.2%、29.7%。其中，柬埔寨的水稻生產用水量將分別減少約17.7%、23.9%、28.9%；老撾的水稻生產用水量分別減少約18.5%、23.2%、30.9%；泰國的水稻生產用水量分別減少約13.8%、22.5%、25.2%；越南的水稻生產用水量別減少28.4%、32.3%、38.8%。同一排放情景下，未來近期和未來遠期2個時期水稻生產用水量在低排放情景SSP1.2.6和中等排放情景SSP2.4.5下的變化并不大，在高排放情景SSP5.8.5下，未來近期到未來遠期各國的水稻生產用水量明顯減少。

3.2 氣象要素與水稻用水量的關系

各氣象要素對水稻生產用水量有著重要的影響，通過計算各國家不同季節水稻生產用水量和氣象要素的相關系數，分析了不同時期內各氣象要素和水稻生產用水量的相關關系。

3.2.1 最低溫度與水稻生產用水量的關系

在AquaCrop模型中，當溫度低于作物生長所需的最低限制溫度時，作物氣孔導度將會降低，作物蒸騰將會減少，這一過程通過冷應力系數KsTr，x來模擬。本研究計算了水稻生產用水量和生長期內最低溫度的相關系，以探究最低溫度對水稻生產用水量的影響（圖6—12）。

在歷史時期，瀾湄流域下游各制度下水稻生長期內平均最低溫度為柬埔寨23.7 ℃、老撾夏季21.1 ℃、老撾冬春季15.9 ℃、泰國23.6 ℃、越南夏季23.6 ℃、越南冬春季21.7 ℃，越南晚季23 ℃。除柬埔寨和老撾水稻生產用水量與最低溫度之間無明顯相關外，其余各制度下水稻生產用水量和最低溫度之間都表現出一定的負相關，但相關程度并不高，其中越南冬春季水稻生產用水量和最低溫度的相關系數為-0.744（P<0.01）。在未來近期，幾乎所有制度下的水稻生產用水量和最低溫度之間都表現出顯著的負相關關系，其中越南晚季水稻的用水量和最低溫度相關性最高，相關系數為-0.947（P<0.01）。所有在夏季種植的水稻其用水量和最低溫度的相關性隨著排放情景的升高而增強。而老撾冬春季水稻在低排放情景下呈正相關，在中等排放情景下相關性減弱，在高排放情景下轉為負相關；越南冬春季水稻用水量和最低溫度的相關性隨著排放情景的升高而減弱。在未來遠期，水稻生產用水量和最低溫度之間幾乎都表現為負相關關系，在高排放情景下，柬埔寨水稻用水量和最低溫度相關系數達到了-0.951（P<0.01）。水稻生產用水量和最低溫度的相關性同樣隨排放情景的升高而增強，但二者的相關程度在低、中等排放情景下往往低于未來近期，在高排放情景下高于未來近期。此外，不同于其余國家，老撾冬春季水稻生長期內平均最低溫度遠低于其余種植制度下的平均最低溫度，在低排放情景SSP1.2.6下，水稻生產用水量與最低溫度在未來2個時期內均表現出顯著的正相關關系。

3.2.2 最高溫度與水稻生產用水量的關系

在作物的生長過程中，除低溫限制外，當溫度高于限制溫度時，作物葉面會因過多失水從而關閉氣孔，導致作物蒸騰減弱。通過計算水稻生產模擬用水量和水稻生長期內最高溫度的相關系數，分析二者之間的相關性，探究最高溫度對水稻生產用水量的影響，結果見圖13—19。

在歷史時期，瀾湄流域下游各制度下水稻生長期內平均最高溫度為柬埔寨30.2 ℃、老撾夏季27.5 ℃、老撾冬春季27 ℃、泰國30 ℃、越南夏季29.1 ℃、越南冬春季29.7 ℃，越南晚季28.8℃，水稻生產用水量和最高溫度之間無明顯的相關性，僅老撾冬春季水稻生產用水量和最高溫度的相關系數為-0.321（P<0.01）。在未來近期，總體來說，水稻生產用水量和最高溫度之間表現為明顯的負相關關系，相關性最高為老撾夏季水稻，其相關系數為-0.928（P<0.01）。特別的，老撾冬春季水稻生產用水量和最高溫有顯著的正相關，相關系數為0.824（P<0.01）。在未來遠期，老撾夏季和泰國水稻生產用水量和最高溫度在3種排放情景下為顯著的負相關，且相關程度在SSP5.8.5情景下高于SSP2.4.5高于SSP1.2.6，其中老撾夏季和泰國水稻生產用水量和最高溫度在SSP5.8.5情景下的相關系數分別為-0.925（P<0.01）和0.896（P<0.01）。老撾冬春、越南夏季、越南晚季和越南冬春季水稻生產用水量和最高溫度在SSP1.2.6情景下為顯著的正相關，相關系數分別為0.832、0.706、0.625、0.458，隨著排放情景升高，在中等排放情景下其正相關關系減弱，繼而在高等排放情景下轉為顯著的負相關關系，其相關系數分別為-0.79、-0.902、-0.89、-0.646。柬埔寨水稻生產用水量和最高溫度在SSP5.8.5情景下呈顯著負相關，相關系數為-0.926（P<0.01）。對比未來近期和未來遠期，各國水稻生產用水量和生長期內最高溫度的相關性在未來遠期普遍強于未來近期。

3.2.3 降水與水稻生產用水量的關系

在AquaCrop模型中，作物根區缺水和水分過量都可能減少作物的蒸騰作用，這一過程通過水分脅迫系數Ks來進行模擬。降水作為模型水分補給的主要途徑，會影響作物根區的含水量從而影響作物蒸騰。本研究通過計算水稻生產用水量和水稻生長期內降水的相關系數，分析其在不同時期不同情境下的相關關系（圖20—26）。

在歷史時期內，僅越南冬春季和晚季水稻生產用水量和降水之間有顯著的負相關（r=-0.744，r=-0.623）外，其余各國各季度水稻生產用水量和降水之間均無明顯相關關系。在未來近期，柬埔寨、越南夏季、越南冬春和越南晚季的水稻生產用水量和降水在SSP1.2.6、SSP2.4.5、SSP5.8.5情景下都表現出一定程度的負相關，其中越南晚季水稻產量和降水之間的相關性最高，在SSP1.2.6、SSP2.4.5、SSP5.8.5情景下的相關系數分別為-0.542、-0.438、-0.622。老撾冬春稻在SSP5-8.5情景夏表現出負相關，其相關系數為-0.537。在未來遠期，越南冬春稻和晚季稻在3種情境下均表現出一定的負相關，其中在SSP2.4.5情景下越南冬春季水稻生產用水量和降水之間的相關性最好，相關系數為-0.567，老撾冬春稻在SSP1.2.6和SSP2.4.5情景下為顯著負相關，相關系數分別為-0.536、-0.634。老撾夏季和泰國僅在SSP5.8.5情景下表現出一定的負相關，相關系數分別為-0.629、-0.409?？傮w來說，越南冬春稻和晚季稻的生產用水量和降水之間的相關性要強于其他國家。

3.2.4 CO2濃度與水稻生產用水量的關系

CO2是影響全球氣候變化主要溫室氣體，大氣中CO2濃度升高不僅會加劇全球變暖，還會減小作物、植被的氣孔開放度，導致氣孔阻力增加，從而減少作物、植被蒸騰。通過計算瀾湄流域下游模擬水稻生產用水量和年平均CO2濃度的相關系數，本研究分析了水稻生產用水量和CO2濃度之間的相關關系（圖27—33）。

在歷史時期，老撾夏季、老撾冬春季、越南夏季、越南晚季和越南冬春季水稻生產用水量和CO2濃度之間均表現出了顯著的負相關，其中老撾夏季相關性最強，相關系數為-0.657（P<0.01）。在未來近期，幾乎所有國家各季度的水稻生產用水量和CO2濃度之間都表現出顯著的負相關，其中越南晚季水稻生產用水量和CO2濃度的相關性達到了-0.949（P<0.01）。同時，通過對比同一季度不同情景下水稻生產用水量和CO2濃度的相關系數發現，其相關性隨著排放情景的升高而逐漸增強。在未來遠期，在SSP1.2.6情景下，柬埔寨、老撾夏季和泰國水稻生產用水量和CO2濃度存在一定的負相關，其相關系數分別為-0.409、-0.365、-0.365；在SSP2.4.5情景下，僅越南夏季水稻生產用水量和CO2濃度不相關外，其余各國水稻生產用水量和CO2濃度均呈現出顯著的負相關，其中越南晚季水稻生產用水量和CO2濃度的相關性最強，相關系數為0.714；在SSP5.8.5情景下，各國所有季節的水稻生產用水量和CO2濃度都表現出極為顯著的負相關關系，泰國的相關系數為-0.875（P<0.01），其余國家各季度下的相關系數都小于-0.9（P<0.01）。

3.3 討論

氣候變化對作物生產的影響在物候學上是復雜而模糊的，AquaCrop模型是一個相對簡化的模型，旨在為全球更廣泛區域內研究環境和管理對作物生產的影響提供更便利的方法。由于研究區內經濟技術以及數據訪問等限制，所能獲取的各方面數據資源有限，部分模型參數采用FAO在標準情況下給予的參考值。同時，本研究對水稻生產用水量的模擬是基于AquaCrop模型的Python開源版本AquaCrop.OSPy來實現，AquaCrop.OSPy主要考慮的是氣候變化對作物生產的影響，暫時沒有納入土壤肥力脅迫和鹽分脅迫，而氣候變暖致使海平面上升，瀾湄流域下游的鹽分入侵正嚴重威脅水稻生產環境［32-33］，盡管模擬結果表現較好，但模型在一定程度上是缺乏代表性的。

此外，在分析氣象要素對水稻生產用水量的影響時方法相對單一，忽略了各氣候要素之間的相互作用。目前，許多方法可用于作物生產用水量的模擬計算和分析，如Hargreaves模型、灰色神經網絡模型和利用遙感影像方法等。隨著理論、科技的不斷進步，對作物生產用水量的模擬和計算精度將進一步提升，灌溉制度和措施也會更加優化完善，水資源利用率將進一步提高。

4 結論

本研究以瀾滄江-湄公河流域下游為研究區，基于ERA5.Land數據集和最新的CMIP6氣候預測數據，選用SSP1.2.6、SSP2.4.5、SSP5.8.5三種排放情景，AquaCrop模型將非生產性土壤蒸發和生產性作物蒸騰分離，將水稻生長期內蒸騰總量作為水稻生產用水量，模擬了瀾滄江-湄公河流域下游歷史和未來時期水稻生產用水量，分析了水稻生產用水量的變化情況，及其與溫度、降水和CO2濃度之間的相關關系，主要結論如下。

a）在1982—2020年的歷史時期瀾湄流域下游水稻生產平均用水量為296.4 mm，其中柬埔寨、老撾、泰國和越南的水稻生產平均用水量分別為302.9、347.2、308.0、258.8 mm。瀾湄流域下游水稻生產用水量總體呈現逐年減少趨勢，水稻生產的用水量在空間上表現為北多南少。

b）未來時期，瀾湄流域下游水稻生產的用水量相對于歷史時期有所減少，在SSP5.8.5情景下減少趨勢更加明顯，到未來遠期瀾湄流域下游水稻生產用水量將減少29.7%，其中柬埔寨、老撾、泰國和越南水稻生產的用水量將分別減少28.9%、30.9%、25.5%、38.8%。

c）瀾湄流域下游水稻生產用水量在未來時期減少與4個氣象變量（最低溫度、最高溫度、降水和CO2濃度）有一定的相關關系。冬春季水稻生產用水量在低排放情景下和最低溫度、最高溫度表現為正相關，在其余種植制度下則總體表現為負相關關系。越南冬春季和晚季水稻生產用水量和降水之間的相關性強于其他種植制度。水稻生產用水量和CO2濃度呈現顯著的負相關關系，在SSP5.8.5情景下的未來遠期，泰國水稻生產用水量和CO2濃度相關系數為-0.875（P<0.01），其余國家各季度下水稻生產用水量和CO2濃度相關系數都小于-0.9（P<0.01）。

參考文獻：

［1］楊敏慎，劉曉雨，郭輝.氣候變暖和CO2濃度升高對農作物的影響［J］.江蘇農業學報，2021，37（1）：246-258.

［2］朱永彬，白冰，劉昌新，等.氣候變化對我國農作物產量影響及氣候損失估算［J］.統計觀察，2019，35（24）：85-89.

［3］LEE S，CHOI J，HUR S，et al.Food.centric interlinkages in agricultural food.energy.water nexus under climate change and irrigation management［J］.Resources Conservation and Recycling，2020，163.DOI：10.1016/j.resconrec.2020.105099.

［4］WANG J，WANG X，XIN L，et al.The Impact of Climate Change and Grain Planting Structure Change on Irrigation Water Requirement for Main Grain Crops in Mainland China［J］.Land，2022，11（12）：11-23.

［5］DERYNG D，CONWAY D，RAMANKUTTY N，et al.Global crop yield response to extreme heat stress under multiple climate change futures［J］.Environmental Research Letters，2014，9（3）：34-41.

［6］SáNCHEZ B，RASMUSSEN A，PORTER J R.Temperatures and the growth and development of maize and rice：a review［J］.Global Change Biology，2014，20（2）：408-417.

［7］熊玉江，喬偉，姚付啟，等.氣候變化和灌溉方式對稻田需水量的影響［J］.長江科學院院報，2018，35（4）：76-80.

［8］韓宇平，王朋，王富強.氣候變化下淮河區主要作物需水量變化特征［J］.灌溉排水學報，2013，32（5）：114-117.

［9］王衛光，彭世彰，孫風朝，等.氣候變化下長江中下游水稻灌溉需水量時空變化特征［J］.水科學進展，2012，23（5）：656-664.

［10］運曉博，湯秋鴻，徐錫蒙，等.氣候變化對瀾湄流域上下游水資源合作潛力的影響［J］.氣候變化研究進展，2020，16（5）：555-563.

［11］陳興茹，王興勇.湄公河流域洪旱災害損失分析［J］.水利經濟，2019，37（1）：54-58.

［12］MAINUDDIN M，KIRBY M，HOANH C T.Adaptation to climate change for food security in the lower Mekong Basin［J］.Food Security，2011，3（4）：433-450.

［13］ZHANG L，SONG W，SONG W.Assessment of Agricultural Drought Risk in the Lancang.Mekong Region，South East Asia［J］.International Journal of Environmental Research and Public Health，2020，17（17）：61-73.

［14］KANG H，SRIDHAR V，MAINUDDIN M，et al.Future rice farming threatened by drought in the Lower Mekong Basin［J］.Scientific Reports，2021，11（1）：93-103.

［15］DONG Z Q，LIU H，BAI YB，et al.Future projection of seasonal drought characteristics using CMIP6 in the Lancang.Mekong River Basin［J］.Journal of Hydrology，2022，610.DOI：10.1016/j.jhydrol.2022.127815.

［16］THILAKARATHNE M，SRIDHAR V.Characterization of future drought conditions in the Lower Mekong River Basin［J］.Weather and Climate Extremes，2017，17（10）：47-58.

［17］YAMAUCHI K.Climate change impacts on agriculture and irrigation in the Lower Mekong Basin［J］.Paddy and Water Environment，2013，12（S2）：227-240.

［18］PAIBOONVORACHAT C，KITYUTTACHAI K，MINAMI I，et al.Irrigation Database Improvement for the Lower Mekong Basin ［R］.2018.

［19］BALLESTER C，VOTE C，HORNBUCKLE J，et al.Effects of Frequency of Irrigation on Dry.Season Furrow.Irrigated Maize and Peanut Production in the Rice.Growing Lowlands of the Lower Mekong Basin［J］.Agriculture，2019，9（6）.DOI：10.3390/agriculture9060128.

［20］YOSHIDA Y，LEE H S，TRUNG B H，et al.Impacts of Mainstream Hydropower Dams on Fisheries and Agriculture in Lower Mekong Basin［J］.Sustainability，2020，12（6）.DOI：10.3390/su12062408.

［21］LY K，METTERNICHT G，MARSHALL L.Linking Changes in Land Cover and Land Use of the Lower Mekong Basin to Instream Nitrate and Total Suspended Solids Variations［J］.Sustainability，2020，12（7）：29-42.

［22］ZHANG B H，ZHANG L，GUO H D，et al.Drought impact on vegetation productivity in the Lower Mekong Basin［J］.Internitional Journal of remotesensing，2014，35（8）：2835-2856.

［23］何大明.瀾滄江：湄公河水文特征分析［J］.云南地理環境研究，1995，7（1）：58-74.

［24］Mekong River Commission，Ministry of Water Resourse of China.Joint Observation and Evaluation of the Emergency Water Supplement from China to the Mekong River［A］.2016.

［25］LIU H.Flood Prevention and Drought Relief in Mekong River Basin［M］.Springer，2021.

［26］CHUNG S，TAKEUCHI J，FUJIHARA M，et al.Flood damage assessment on rice crop in the Stung Sen River Basin of Cambodia［J］.Paddy and Water Environment，2019，17（2）：255-263.

［27］HUANG Z，LI Z，ZHOU K，et al.Rice production and hybrid rice development in thailand［J］.Hybrid Rice，2018，33（5）：68-71.

［28］SINNARONG N，CHEN C C，MCCARL B，et al.Estimating the potential efects of climate change on rice production in Thailand［J］.Paddy and Water Environment，2019，17（4）：761-769.

［29］侯儷.老撾大米出口現狀及競爭力研究［D］.昆明：云南大學，2019.

［30］章忠貴，高前寶，范凌，等.越南雜交水稻市場現狀及發展對策探討［J］.現代農業科技，2017（22）：42-43.

［31］TANG Y Z，HUO J J，ZHU D J，et al.Spatio.Temporal V ariation of Precipitation and Evaporation on the Tibetan Plateau and Their Influence on Regional Drought［J］.Atmosphere，2022，13.DOI：10.3390/atmos13081323.

［32］SMAJGL A，TOAN T Q，NHAN D K，et al.Responding to rising sea levels in the Mekong Delta［J］.Nature Climate Change，2015，5（2）：167-174.

［33］VU D T，YAMADA T，ISHIDAIRA H.Assessing the impact of sea level rise due to climate change on seawater intrusion in Mekong Delta，Vietnam［J］.Water Science and Technology，2018，77（5/6）：1632-1639.