海南省南繁水稻制種水分供求情況分析

摘""要：基于海南省南繁水稻制種大市（縣）東方市和樂東縣水稻生育期資料和1971—2020年氣象數(shù)據(jù)，采用Penman-Monteith公式、作物系數(shù)法和美國農(nóng)業(yè)部土壤保持局推薦的方法，得出兩市（縣）南繁水稻制種大田期及其前、中和后期3個(gè)生育階段的需水量（ETc）、有效降水量（Pe）和缺水量（Dw），進(jìn)而分析兩市（縣）南繁水稻制種水分供求變化特征。采用多元回歸方法對影響南繁水稻制種ETc和Dw的主要?dú)庀笠蜃舆M(jìn)行定量研究。結(jié)果表明：東方市和樂東縣南繁水稻制種大田期ETc"50"a均值分別為400.2、322.3"mm，Pe均值分別為45.7、84.3"mm，Dw均值分別為354.5、238.0"mm。兩市（縣）南繁水稻制種大田期ETc和Dw均主要分布在中、前期，占大田期ETc和Dw的比例之和為72.4%～85.2%。兩市（縣）南繁水稻制種大田期Pe均主要分布在后、中期，占大田期Pe的比例之和分別為78.2%和85.6%。近50"a東方市南繁水稻制種大田期及其各生育階段ETc和Dw均呈增加趨勢（其中大田期增加顯著），Pe均呈減少趨勢。樂東縣南繁水稻制種ETc、Pe和Dw隨時(shí)間變化趨勢與東方市大致相反。近50"a東方市南繁水稻制種的缺水形勢顯著加重，而樂東縣南繁水稻制種的缺水形勢有所緩解但仍嚴(yán)峻，需積極采取有效措施減輕氣候變化給兩市（縣）南繁水稻制種產(chǎn)業(yè)帶來的影響。近50"a高溫不顯著升高且低溫極顯著升高或高溫不顯著升高且日照時(shí)數(shù)不顯著增加是引起東方市南繁水稻制種大多數(shù)生育階段ETc和Dw增加的主要原因；而日照時(shí)數(shù)不顯著減少且平均風(fēng)速極顯著降低或日照時(shí)數(shù)不顯著減少且降水量不顯著增加是引起樂東縣南繁水稻制種大多數(shù)生育階段ETc和Dw減少的主要原因。本研究結(jié)果為氣候變化背景下海南省南繁水稻制種合理用水、制定科學(xué)灌溉方案提供有力支撐。

關(guān)鍵詞：南繁水稻制種；需水量；有效降水量；缺水量；多元回歸；海南中圖分類號(hào)：P49；S511""""""文獻(xiàn)標(biāo)志碼：A

Analysis"of"Water"Supply"and"Demand"of"Rice"Seed"Reproduction"in"Hainan，"China

ZOU"Haiping1，2，"CHEN"Xiaomin1，2，"ZHANG"Mingjie1，2，"LYU"Run1，2，"LIN"Shaowu1，2，"LI"Weiguang1，2，"BAI"Rui1，2，"CHEN"Yanxi2，3

1."Hainan"Provincial"Climate"Center，"Haikou，"Hainan"570203，"China;"2."Hainan"Provincial"Key"Laboratory"of"South"China"Sea"Meteorological"Disaster"Prevention"and"Mitigation，"Haikou，"Hainan"570203，"China;"3."Hainan"Meteorological"Exploration"Center，"Haikou，"Hainan"570203，"China

Abstract："Based"on"the"data"of"rice"growth"period"in"Dongfang"and"Ledong，"main"producers"for"rice"seed"reproduction"in"Hainan，"and"the"meteorological"data"from"1971"to"2020，"the"Penman-Monteith"formula"and"crop"coefficient"method"andnbsp;the"method"recommended"by"the"Soil"Conservation"Service"of"the"United"States"Department"of"Agriculture"were"used"to"obtain"the"water"requirement"（ETc），"effective"precipitation"（Pe）"and"water"deficit"（Dw）"of"rice"seed"reproduction"in"field"period"and"early，"middle"and"late"stage"of"field"period."The"change"characteristics"of"water"supply"and"demand"of"rice"seed"reproduction"in"the"two"regions"were"then"analyzed."The"50-year"average"value"of"ETc"in"field"period"of"rice"seed"reproduction"in"Dongfang"and"Ledong"was"400.2"mm"and"322.3"mm，"respectively，"the"mean"value"of"Pe"was"45.7"mm"and"84.3"mm，"respectively，"and"the"mean"value"of"Dw"was"354.5"mm"and"238.0"mm，"respectively."The"ETc"and"Dw"in"field"period"of"rice"seed"reproduction"were"mainly"distributed"in"the"middle"and"early"stages，"accounting"for"72.4%"to"85.2%"of"ETc"and"Dw"during"the"field"period."The"Pe"in"field"period"of"rice"seed"reproduction"was"mainly"distributed"in"the"late"and"middle"stages，"accounting"for"78.2%"and"85.6%"of"Pe"during"the"field"period，"respectively."In"the"past"50"years，"the"ETc"and"Dw"of"rice"seed"reproduction"in"field"period"and"all"growth"stages"of"field"period"in"Dongfang"showed"an"increasing"trend"（the"increase"in"field"period"was"significant），"while"the"Pe"showed"a"decreasing"trend."The"variation"trend"of"ETc，"Pe"and"Dw"of"rice"seed"reproduction"with"time"in"Ledong"was"roughly"opposite"to"that"in"Dongfang."In"the"past"50"years，"the"water"shortage"situation"of"rice"seed"reproduction"in"Dongfang"was"significantly"worsening，"while"the"water"shortage"situation"in"Ledong"was"easing"to"some"extent"but"remained"severe."Effective"measures"should"be"taken"to"mitigate"the"impact"of"climate"change"on"the"rice"seed"reproduction."In"the"past"50"years，"the"main"reasons"for"the"increase"of"ETc"and"Dw"during"most"growth"stages"of"rice"seed"reproduction"in"Dongfang"were"the"insignificant"increase"in"high"temperature"and"extremely"significant"increase"in"low"temperature，"or"the"insignificant"increase"in"high"temperature"and"sunshine"hours."The"main"reasons"for"the"decrease"of"ETc"and"Dw"during"most"growth"stages"of"rice"seed"reproduction"in"Ledong"were"the"insignificant"decrease"in"sunshine"hours，"extremely"significant"decrease"in"average"wind"speed，"or"insignificant"decrease"in"sunshine"hours"and"insignificant"increase"in"precipitation."The"results"of"this"study"provide"strong"support"for"rational"water"use"and"scientific"irrigation"scheme"for"rice"seed"reproduction"in"Hainan"province"under"the"background"of"climate"change.

Keywords："rice"seed"reproduction;"water"requirement;"effective"precipitation;"water"deficit;"multiple"regression;"Hainan

DOI："10.3969/j.issn.1000-2561.2025.04.024

海南省是我國重要的雜交水稻南繁制種基地，年制種面積達(dá)到13"333.33"hm2左右，生產(chǎn)的水稻種子可供200萬～267萬"hm2大田種植。水資源作為限制農(nóng)業(yè)發(fā)展的重要因素之一，是農(nóng)業(yè)生產(chǎn)中最重要的資源支撐和保證，與農(nóng)業(yè)生產(chǎn)和糧食安全息息相關(guān)[1]。海南省南繁水稻制種主要是春季制種，生育期集中在冬、春季[2]。而冬、春季為海南省非汛期，其降水量僅占年降水量的23.4%，加之南繁水稻制種主要分布在海南省年降水量最少和次少的西部、南部[3-4]，導(dǎo)致南繁水稻制種水分供需矛盾突出，尤其是移栽期和授粉期，幾乎每年均有部分制種區(qū)域前期因缺水導(dǎo)致不能移栽，授粉期缺水導(dǎo)致嚴(yán)重減產(chǎn)甚至絕收。如2016年因干旱時(shí)間過長、水庫缺水、地下水位降低，樂東縣佛羅鎮(zhèn)青山村有33.3"hm2制種田絕收[5]。因此，研究如何提高海南省南繁水稻制種灌溉用水利用率，對保障南繁水稻制種科學(xué)高效生產(chǎn)意義重大?？茖W(xué)確定作物需水量、有效降水量和缺水量，分析作物水分供求情況是進(jìn)行水資源優(yōu)化管理、合理制定作物灌溉方案的重要依據(jù)。全球氣候變暖以及隨之改變的降水時(shí)空格局將會(huì)直接影響作物生長發(fā)育和耗水過程，進(jìn)而影響作物需水量和缺水量[6]。因此很有必要對南繁水稻制種水分供求變化特征及氣候歸因進(jìn)行分析。

國內(nèi)外很多學(xué)者分析了不同作物需水量、缺水量和與其密切相關(guān)的參考作物蒸散量的變化特征及歸因[7-19]。在研究方法上，主要利用聯(lián)合國糧食及農(nóng)業(yè)組織（Food"and"Agriculture"Organization"of"the"United"Nations，"FAO）推薦的Penman-"Monteith公式和作物系數(shù)法計(jì)算需水量，該方法考慮了影響蒸散的大氣物理特性和植物生理機(jī)制，具有很好的物理基礎(chǔ)[20]；缺水量為需水量與同時(shí)期有效降水量之差[21]，有效降水量的計(jì)算主要采用美國農(nóng)業(yè)部土壤保持局（SCS）推薦的方法，該方法是目前眾多的有效降水量計(jì)算方法中得到公認(rèn)和普遍推廣的方法之一，其有效性已在許多學(xué)者的研究中予以證明[22]；作物水分供求變化歸因分析主要采用多元回歸方法，該方法能夠定量不同氣象要素對作物水分供求變化的影響[17]。目前尚無專門針對海南省南繁水稻制種水分供求變化特征及氣候歸因的研究報(bào)道。海南省南繁水稻制種面積最大的市（縣）為東方市，約6000"hm2，其次為樂東縣，這兩個(gè)市（縣）的南繁水稻制種面積占全省的比例高達(dá)八成左右[4]。因此，本研究基于東方市和樂東縣南繁水稻制種生育期資料和氣象數(shù)據(jù)，采用Penman-Monteith公式和作物系數(shù)法計(jì)算南繁水稻制種需水量，同時(shí)結(jié)合有效降水量獲得南繁水稻制種缺水量，分析兩市（縣）南繁水稻制種水分供求變化特征，并對其主要影響因素進(jìn)行探討，為氣候變化背景下海南省南繁水稻制種合理用水、制定科學(xué)灌溉方案提供有力支撐。

1""材料與方法

1.1""材料

氣象資料來自海南省氣象信息中心，包括東方市和樂東縣1971—2020年逐日的均溫、高溫、低溫、日照時(shí)數(shù)、平均風(fēng)速、水汽壓、氣壓及經(jīng)緯度和海拔高度。生育期資料來自海南水稻制種產(chǎn)量保險(xiǎn)調(diào)查數(shù)據(jù)，包括東方市和樂東縣2021年春季水稻制種移栽期、孕穗期、開花期和成熟期。

1.2""方法

1.2.1""生育期確定""本文研究的是南繁水稻制種大田期（移栽至成熟）的水分供求情況，根據(jù)水稻需水規(guī)律將其劃分為3個(gè)生育階段，即前期（移栽至孕穗）、中期（孕穗至開花）和后期（開花至成熟），各生育期出現(xiàn)時(shí)間為多個(gè)品種實(shí)測平均值。其中移栽期涉及父本和母本移栽期，二者相差2～3"d，取二者平均值作為移栽期。確定的水稻生育期見表1。

1.2.2""需水量計(jì)算""作物需水量是指作物在適宜的土壤水分和肥力水平下，經(jīng)過正常生長發(fā)育，獲得高產(chǎn)時(shí)的植株蒸騰、株間蒸發(fā)以及構(gòu)成植株體的水量之和。由于構(gòu)成植株體的水量很少，可忽略不計(jì)，因此在實(shí)際計(jì)算中作物需水量為植株蒸騰量與株間蒸發(fā)量之和。南繁水稻制種需水量采用參考作物蒸散量和作物系數(shù)法進(jìn)行逐日計(jì)算，該方法充分考慮了植株蒸騰量與株間蒸發(fā)量，其中參考作物蒸散量是指生長高度為0.12"m、固定葉面阻力為70"s/m、反射率為0.23的參考作物的蒸發(fā)蒸騰量，類似于面積廣闊、供水充分、高度一致、生長旺盛、完全遮蔽地面的草地的蒸發(fā)蒸騰量[20]。南繁水稻制種各生育階段的需水量由生育階段內(nèi)逐日需水量累加得出：

式中，ETc為南繁水稻制種某生育階段需水量（mm）；ETci為日需水量（mm/d）；i為生育階段天數(shù)（d）；ET0i為日參考作物蒸散量（mm/d），采用Penman-Monteith公式計(jì)算[20]；Kc為作物系數(shù)，不同發(fā)育階段Kc不同，根據(jù)FAO推薦的Kc值確定南繁水稻制種大田期前期、中期和后期的作物系數(shù)分別為1.05、1.20、0.75。由于海南實(shí)際情況與FAO推薦的Kc值需滿足的標(biāo)準(zhǔn)條件相差較遠(yuǎn)，因此根據(jù)當(dāng)?shù)赝寥篮蜌夂驐l件參考文獻(xiàn)[23]的方法對Kc值進(jìn)行修正。

1.2.3""有效降水量計(jì)算""有效降水量是指總降水量中能夠被作物有效利用的降水量[16]。南繁水稻制種有效降水量采用美國農(nóng)業(yè)部土壤保持局（SCS）推薦的方法進(jìn)行逐日計(jì)算，利用該方法計(jì)算的有效降水量充分考慮了無效降水量（含徑流量、蒸發(fā)量和深層滲漏量），其有效性已在許多學(xué)者的研究中予以證明[22]，南繁水稻制種各生育階段的有效降水量由生育階段內(nèi)逐日有效降水量累加得出：

式中，Pe為南繁水稻制種某生育階段有效降水量（mm）；Pei為日有效降水量（mm/d）；Pi為日降水量（mm/d）。

1.2.4""缺水量計(jì)算""缺水量為作物生長過程中需依靠灌溉補(bǔ)充的水量，南繁水稻制種各生育階段缺水量等于需水量減去有效降水量：，式中，Dw為南繁水稻制種某生育階段缺水量（mm）。Dwgt;0，表示該生育階段作物缺水；Dw=0，表示水分供需平衡；Dwlt;0，表示水分盈余。

1.2.5""氣候傾向率""氣候傾向率是指某一要素每10"a的變化速率，用于表征該要素隨時(shí)間的變化趨勢，正值表示該要素呈增加趨勢，負(fù)值表示呈減少趨勢。本研究用Xk表示樣本量為k的南繁水稻制種ETc或Pe或Dw或各氣象因子，用tk表示Xk所對應(yīng)的年份（即1971，1972，…2020），利用Excel軟件建立Xk與tk之間的一元線性回歸方程，并采用F檢驗(yàn)法進(jìn)行顯著性檢驗(yàn)[24]，以方程

回歸系數(shù)（a）的10倍作為南繁水稻制種ETc或Pe或Dw或各氣象因子的氣候傾向率：Xk=atk+b"（k=1，"2，"3，…50）。式中，b為回歸常數(shù)。

1.2.6""多元回歸分析""由于影響南繁水稻制種ETc和Dw的氣象因子較多，因此利用SPSS"16.0軟件多元回歸方法計(jì)算各氣象因子對ETc和Dw的貢獻(xiàn)率，確定各氣象因子對南繁水稻制種ETc和Dw變化的影響。因各氣象因子的單位、量綱和范圍不盡相同，對各個(gè)氣象因子及南繁水稻制種ETc和Dw進(jìn)行標(biāo)準(zhǔn)化處理，標(biāo)準(zhǔn)化后的數(shù)據(jù)序列按照以下公式計(jì)算貢獻(xiàn)率[19]：

式中，為南繁水稻制種ETc或Dw的標(biāo)準(zhǔn)化值，X1，X2，X3…Xn為各氣象因子的標(biāo)準(zhǔn)化值，a1，a2，a3…an為標(biāo)準(zhǔn)化后回歸系數(shù)，Wj為某一氣象因子對南繁水稻制種ETc或Dw的貢獻(xiàn)率。

1.2.7""偏相關(guān)分析""在分析氣象因子與南繁水稻制種ETc和Dw的關(guān)系時(shí)，氣象因子之間相互影響，存在一定的相關(guān)性，可能掩蓋氣象因子與南繁水稻制種ETc和Dw的真實(shí)關(guān)系。偏相關(guān)分析能在有效控制其他變量影響條件下揭示兩變量之間的線性相關(guān)程度。因此利用SPSS"16.0軟件偏相關(guān)方法計(jì)算偏相關(guān)系數(shù)分別確定南繁水稻制種ETc和Dw與各氣象因子的相關(guān)性，并采用t檢驗(yàn)法進(jìn)行顯著性檢驗(yàn)。

2""結(jié)果與分析"

2.1""南繁水稻制種水分供求的時(shí)間分布

1971—2020年東方市和樂東縣南繁水稻制種大田期各生育階段需水量（ETc）、有效降水量（Pe）和缺水量（Dw）均值及占大田期的比例如圖1所示。由圖1A可看出，東方市南繁水稻制種ETc大田期中期（164.7"mm）gt;前期（125.4"mm）gt;后期（110.2"mm），大田期ETc為400.2"mm；Pe表現(xiàn)為后期（22.3"mm）gt;中期（13.7"mm）gt;前期（9.7"mm）的規(guī)律，大田期Pe為45.7"mm；大田期各生育階段Dw分布情況與ETc一致，大田期Dw為354.5"mm。樂東縣大田期各生育階段ETc和Pe分布情況與東方市一致，大田期ETc和Pe分別為322.3、84.3"mm；大田期各生育階段Dw分布情況與東方市略有差異，表現(xiàn)為前期（101.6"mm）gt;中期（99.6"mm）gt;后期（36.8"mm），大田期Dw為238.0"mm。為更直觀地體現(xiàn)南繁水稻制種大田期ETc、Pe和Dw在各生育階段分布情況，計(jì)算了各生育階段各指標(biāo)值占大田期的比例（圖1B）。由圖1B可看出，東方市南繁水稻制種大田期中期、前期ETc占大田期ETc的比例之和為72.4%，Dw相應(yīng)的占比和為75.3%，說明東方市南繁水稻制種大田期ETc和Dw均主要分布在中期、前期；同理，大田期Pe主要分布在后期、中期，二者占比和為78.2%。樂東縣大田期ETc和Dw亦主要分布在中期、前期，二者占比和分別為74.9%和85.2%；大田期Pe亦主要分布在后期、中期，二者占比和為85.6%。

1971—2020年南繁水稻制種大田期及其各生育階段Dw的正負(fù)情況如圖2所示。由圖2可知，1971—2020年東方市南繁水稻制種大田期及其所有生育階段的Dw均為正值。而樂東縣南繁水稻制種除大田期后期1986年、1999年、2000年和2019年的Dw為負(fù)值外（分別為–2.6、–10.5、–5.0、–10.4"mm），后期其余年份和大田期及其前期、中期每年的Dw均為正值，說明1971—2020年東方市南繁水稻制種大田期所有生育階段每年均缺水，樂東縣稍有差異，僅南繁水稻制種大田期后期極少數(shù)年份水分有盈余。

2.2""南繁水稻制種水分供求變化趨勢

1971—2020年南繁水稻制種水分供求變化趨勢見表2。由表2可知，東方市南繁水稻制種大田期前、中和后期的ETc均呈不顯著增加趨勢，其氣候傾向率分別為2.4、1.9、1.4"mm/10"a，而三者累加致大田期ETc呈顯著增加趨勢，氣候傾向率為5.7"mm/10"a；Pe的變化趨勢與ETc相反，大田期及其前、中和后期均不顯著減少；Dw為ETc與Pe之差，因ETc呈增加趨勢而Pe呈減少趨勢，導(dǎo)致大田期及其前、中和后期的Dw均呈增加趨勢，其中大田期的Dw顯著增加，對應(yīng)的氣候傾向率分別為8.6、2.5、2.7、3.4"mm/10"a。

樂東縣南繁水稻制種ETc、Pe和Dw隨時(shí)間變化趨勢與東方市相反。樂東縣大田期及其前、中和后期ETc均呈不顯著減少趨勢，其氣候傾向率分別為–4.3、–1.8、–1.5、–1.0"mm/10"a；除大田期后期Pe呈減少趨勢外（氣候傾向率為–0.3"mm/10"a），大田期前、中期和大田期Pe均呈增加趨勢；大田期及其前、中和后期Dw均呈不顯著減少趨勢，其氣候傾向率分別為–4.4、–2.1、–1.6、–0.7"mm/10"a。

2.3""南繁水稻制種水分供求變化的歸因分析

利用多元回歸方法計(jì)算各氣象因子對南繁水稻制種ETc和Dw的貢獻(xiàn)率（表3）。東方市和樂東縣南繁水稻制種大田期及其各生育階段的ETc、Dw與相應(yīng)階段的氣象因子（降水量、均溫、高溫、低溫、日照時(shí)數(shù)、水汽壓和平均風(fēng)速）的多元回歸方程的決定系數(shù)（R2）為0.918～0.995，均大于0.900，說明構(gòu)建的方程可靠性高，采用上述7個(gè)氣象因子解釋ETc和Dw的變化原因是科學(xué)合理的。但多元回歸方法計(jì)算的貢獻(xiàn)率均為正值，無法區(qū)分各氣象因子對南繁水稻制種ETc和Dw的正負(fù)貢獻(xiàn)，進(jìn)而無法確定影響ETc和Dw變化的主要?dú)庀笠蜃?。為解決該問題，需結(jié)合各氣象因子與南繁水稻制種ETc和Dw的相關(guān)性（即偏相關(guān)系數(shù)，表3）及各氣象因子的變化趨勢（即氣候傾向率，表4）進(jìn)行分析。

以東方市南繁水稻制種大田期前期ETc為例，由表3可知，各氣象因子中水汽壓對大田期前期ETc變化的貢獻(xiàn)率最高為25.0%，水汽壓與ETc的偏相關(guān)系數(shù)為–0.846且通過極顯著檢驗(yàn)，說明水汽壓與ETc呈極顯著負(fù)相關(guān)。而從表4可知，1971—2020年東方市南繁水稻制種大田期前期水汽壓呈不顯著增加趨勢（氣候傾向率為0.1"hPa/10"a）。水汽壓增加導(dǎo)致ETc減少，而實(shí)際ETc在增加（表2），說明東方市大田期前期水汽壓對ETc變化的貢獻(xiàn)為負(fù)貢獻(xiàn)。各氣象因子中高溫對ETc變化的貢獻(xiàn)率為23.4%，高溫與ETc呈極顯著正相關(guān)（偏相關(guān)系數(shù)為0.422），1971—2020年東方市南繁水稻制種大田期前期高溫呈不顯著升高趨勢（氣候傾向率為0.2"℃/10"a）。高溫升高導(dǎo)致ETc增加，實(shí)際ETc也在增加，說明東方市大田期前期高溫對ETc變化的貢獻(xiàn)為正貢獻(xiàn)。以此類推，可確定東方市和樂東縣南繁水稻制種大田期及其各生育階段各氣象因子對其ETc和Dw變化貢獻(xiàn)的正負(fù)情況，同一類別貢獻(xiàn)按照貢獻(xiàn)率大小由高至低排序，結(jié)果見表5。

通過計(jì)算發(fā)現(xiàn)，正貢獻(xiàn)率排在前2位的氣象因子貢獻(xiàn)率之和占所有正貢獻(xiàn)率之和的47.7%～"91.6%，普遍高于50.0%（表5），說明這2個(gè)氣象因子可以解釋南繁水稻制種ETc和Dw變化的大部分原因，因此定義為主要因子。由此可知，引起東方市南繁水稻制種大田期及其前、中期ETc和Dw增加的主要因子為高溫和低溫或高溫和日照時(shí)數(shù)。從表4可知，上述生育階段的高溫均不顯著升高，而低溫均極顯著升高，日照時(shí)數(shù)均不顯著增加，說明高溫不顯著升高且低溫極顯著升高或高溫不顯著升高且日照時(shí)數(shù)不顯著增加是引起上述生育階段ETc和Dw增加的主要原因；大田期后期的主要因子為低溫和水汽壓或降水量和水汽壓，該生育階段低溫極顯著升高，降水量不顯著減少，水汽壓不顯著降低，說明低溫極顯著升高且水汽壓不顯著降低或降水量不顯著減少且水汽壓不顯著降低是導(dǎo)致大田期后期ETc和Dw增加的主要原因。而樂東縣的情形有所差異，引起樂東縣南繁水稻制種大田期及其各生育階段ETc減少的主要因子均為日照時(shí)數(shù)和平均風(fēng)速，而上述階段的日照時(shí)數(shù)均不顯著減少，平均風(fēng)速均極顯著降低，說明日照時(shí)數(shù)不顯著減少且平均風(fēng)速極顯著降低是導(dǎo)致大田期及其各生育階段ETc減少的主要原因；大田期及其前、中期Dw減少的主要因子為日照時(shí)數(shù)和降水量或日照時(shí)數(shù)和低溫，上述生育階段的日照時(shí)數(shù)均不顯著減少，而降水量不顯著增加，低溫均不顯著升高，說明日照時(shí)數(shù)不顯著減少且降水量不顯著增加或日照時(shí)數(shù)不顯著減少且低溫不顯著升高是引起上述生育階段Dw減少的主要原因。大田期后期Dw減少的主要因子為降水量和低溫，該生育階段降水量不顯著增加，低溫不顯著升高，說明降水量不顯著增加和低溫不顯著升高是引起大田期后期Dw減少的主要原因。

3""結(jié)論與討論

研究表明，1971—2020年東方市和樂東縣南繁水稻制種大田期ETc均值分別為400.2、322.3"mm，大田期各生育階段ETc中期gt;前期gt;后期，符合水稻生育過程中ETc先由小到大，再由大到小的變化規(guī)律[25]。兩市（縣）大田期ETc均主要分布在中、前期，占大田期ETc的比例之和分別為72.4%和74.9%；東方市和樂東縣大田期Pe均值分別為45.7、84.3"mm，大田期各生育階段Pe后期gt;中期gt;前期。這是因?yàn)閮墒校h）南繁水稻制種大田期前、中、后期對應(yīng)的月份均大致分別為3、4、5月，3—5月降水量逐月明顯增加[3]，進(jìn)而導(dǎo)致Pe亦逐月增加。東方市和樂東縣大田期Pe均主要分布在后、中期，占大田期的比例之和分別為78.2%和85.6%；大田期ETc主要分布在中、前期，而中、前期Pe相對較少，占大田期Pe的比例之和分別為52.6%和48.6%，導(dǎo)致南繁水稻制種中、前期缺水嚴(yán)重，兩市（縣）南繁水稻制種大田期中、前期Dw占大田期Dw的比例之和分別為75.3%和85.2%。實(shí)際生產(chǎn)中，樂東縣南繁水稻制種缺水階段主要為移栽期-授粉期[5]，本研究結(jié)果與實(shí)際情況相符。1971—2020年東方市南繁水稻制種大田期所有生育階段每年均缺水，而樂東縣稍有差異，南繁水稻制種大田期后期有極少數(shù)年份水分有盈余。因此，需加強(qiáng)這兩市（縣）南繁水稻制種大田期各生育階段尤其是中、前期的節(jié)水灌溉管理。

1971—2020年東方市南繁水稻制種大田期各生育階段ETc均不顯著增加，大田期ETc顯著增加，氣候傾向率為5.7"mm/10"a；大田期及其各生育階段Pe均不顯著減少；大田期各生育階段Dw均不顯著增加，大田期Dw顯著增加，氣候傾向率為8.6"mm/10"a。樂東縣南繁水稻制種大田期及其各生育階段ETc均不顯著減少；Pe除大田期后期不顯著減少外，大田期及其前、中期均不顯著增加；大田期及其各生育階段Dw均不顯著減少。東方市南繁水稻制種大田期各生育階段Dw在增加，且大田期增加顯著，表明近50"a東方市南繁水稻制種的缺水形勢日益嚴(yán)峻。樂東縣南繁水稻制種大田期及其各生育階段Dw均不顯著減少，說明樂東縣南繁水稻制種缺水形勢有所緩解但仍較為嚴(yán)峻。因此，需采取有效措施，如加強(qiáng)農(nóng)業(yè)灌溉設(shè)施建設(shè)、選用抗旱品種等，以應(yīng)對氣候變化給兩市（縣）南繁水稻制種產(chǎn)業(yè)帶來的影響。

歸因分析表明，高溫不顯著升高且低溫極顯著升高或高溫不顯著升高且日照時(shí)數(shù)不顯著增加是引起東方市南繁水稻制種大田期及其前、中期ETc和Dw增加的主要原因；低溫極顯著升高且水汽壓不顯著降低或降水量不顯著減少且水汽壓不顯著降低是導(dǎo)致東方市大田期后期ETc和Dw增加的主要原因。日照時(shí)數(shù)不顯著減少且平均風(fēng)速極顯著降低是導(dǎo)致樂東縣南繁水稻制種大田期及其各生育階段ETc減少的主要原因；日照時(shí)數(shù)不顯著減少且降水量不顯著增加或日照時(shí)數(shù)不顯著減少且低溫不顯著升高是引起樂東縣南繁水稻制種大田期及其前、中期Dw減少的主要原因，降水量不顯著增加且低溫不顯著升高是引起樂東縣大田期后期Dw減少的主要原因?？傮w而言，高溫不顯著升高且低溫極顯著升高或高溫不顯著升高且日照時(shí)數(shù)不顯著增加是引起東方市南繁水稻制種大多數(shù)生育階段ETc和Dw增加的主要原因；而日照時(shí)數(shù)不顯著減少且平均風(fēng)速極顯著降低或日照時(shí)數(shù)不顯著減少且降水量不顯著增加是引起樂東縣南繁水稻制種大多數(shù)生育階段ETc和Dw減少的主要原因。南繁水稻制種ETc和Dw均是在參考作物蒸散量（ET0）的基礎(chǔ)上通過計(jì)算而得到，且與ET0呈正相關(guān)，陳小敏等[26]采用敏感系數(shù)與氣象因子相對變化率相結(jié)合的方法，研究得出東方市春季（3—5月）ET0增加的主要原因是最低溫和最高溫升高，而樂東縣春季ET0減少的主要原因是日照時(shí)數(shù)減少和平均風(fēng)速降低。本研究結(jié)果與其大致相符，說明采用多元回歸方法確定南繁水稻制種ETc和Dw變化歸因切實(shí)可行。本研究結(jié)果為氣候變化背景下海南省南繁水稻制種合理用水、制定科學(xué)灌溉方案提供有力支撐。

參考文獻(xiàn)

• 段愛旺，"張寄陽."中國灌溉農(nóng)田糧食作物水分利用效率的研究[J]."農(nóng)業(yè)工程學(xué)報(bào)，"2000，"16（4）："41-44.DUAN"A"W，"ZHANG"J"Y."Water"use"efficiency"of"grain"crops"in"irrigated"farmland"in"China[J]."Transactions"of"the"Chinese"Society"of"Agricultural"Engineering，"2000，"16（4）："41-44."（in"Chinese）
• 呂青，"蘇保忠，"王靖，"梅子思."海南（南繁）水稻制種天氣指數(shù)保險(xiǎn)可行性研究[J]."現(xiàn)代農(nóng)業(yè)研究，"2018（1）："5-9，"16.LYU"Q，"SU"B"Z，"WANG"J，"MEI"Z"S."A"feasibility"study"on"the"weather"index"insurance"of"rice"in"Hainan[J]."Modern"Agriculture"Research，"2018（1）："5-9，"16."（in"Chinese）
• 王春乙."海南氣候[M]."北京："氣象出版社，"2014："18-21.WANG"C"Y."Hainan"climate[M]."Beijing："China"Meteorological"Press，"2014："18-21."（in"Chinese）
• 呂青，"柯用春，"何志軍，"王仕明，"唐萍，"陳冠銘."南繁制種水稻基地現(xiàn)狀以及問題分析[J]."農(nóng)村經(jīng)濟(jì)與科技，"2017，"28（20）："24-25.LYU"Q，"KE"Y"C，"HE"Z"J，"WANG"S"M，"TANG"P，"CHEN"G"M."Analysis"of"current"situation"and"problems"of"seed"southern-production"rice"bases[J]."Rural"Economy"and"Science，"2017，"28（20）："24-25."（in"Chinese）

[5]"張維亮，"陳同明，"肖層林."兩系雜交水稻新組合隆兩優(yōu)華占在海南高產(chǎn)高質(zhì)制種技術(shù)[J]."雜交水稻，"2019，"34（1）："33-34，"48."ZHANG"W"L，"CHEN"T"M，"XIAO"C"L."High-yielding"and"high-quality"seed"production"techniques"for"new"two-line"hybrid"rice"combination"Longliangyou"Huazhan"in"Hainan[J]."Hybrid"Rice，"2019，"34（1）："33-34，"48."（in"Chinese）

[6]"張華，"王浩."1967—2017年甘肅省小麥需水量和缺水量時(shí)空特征[J]."干旱區(qū)地理，"2019，"42（5）："1094-1104.ZHANG"H，"WANG"H."Spatial"and"temporal"characteristics"of"water"requirement"and"water"deficit"of"wheat"in"Gansu"province"from"1967"to"2017[J]."Arid"Land"Geography，"2019，"42（5）："1094-1104."（in"Chinese）

[7]"劉鈺，"汪林，"倪廣恒，"叢振濤."中國主要作物灌溉需水量空間分布特征[J]."農(nóng)業(yè)工程學(xué)報(bào)，"2009，"25（12）："6-12.LIU"Y，"WANG"L，"NI"G"H，"CONG"Z"T."Spatial"distribution"characteristics"of"irrigation"water"requirement"for"main"crops"in"China[J]."Transactions"of"the"Chinese"Society"of"Agricultural"Engineering，"2009，"25（12）："6-12."（in"Chinese）

[8]"任修琳，"李宏亮，"張玉虎，"蒲曉，"張立林."2000—2015年三江平原主要作物需水量特征及影響因素分析[J]."干旱區(qū)地理，"2019，"42（4）："854-866.REN"X"L，"LI"H"L，"ZHANG"Y"H，"PU"X，"ZHANG"L"L."Water"requirement"characteristics"and"influencing"factors"of"main"crops"in"the"Sanjiang"plain"from"2000"to"2015[J]."Arid"Land"Geography，"2019，"42（4）："854-866."（in"Chinese）

• 張鳳怡，"遲道才，"陳濤濤."遼寧主要糧食作物生長季需水與降水耦合度分析[J]."中國農(nóng)業(yè)氣象，"2021，"42（9）："746-"760.ZHANG"F"Y，"CHI"D"C，"CHEN"T"T."Assessment"of"coupling"degree"between"water"requirement"of"main"cereal"crops"and"precipitation"in"growing"season"in"Liaoning"province[J]."Chinese"Journal"of"Agrometeorology，"2021，"42（9）："746-760."（in"Chinese）

[10]"KUMAR"S."Reference"evapotranspiration"（ET0）"and"irrigation"water"requirement"of"different"crops"in"Bihar[J]."Journal"of"Agrometeorology，"2017，"19（3）："238-241.

• 曹言，"王杰，"王樹鵬，"黃英，"張雷，"李尤亮."氣候變化下滇中地區(qū)水稻需水量與灌溉需水指數(shù)時(shí)空變化研究[J]."干旱地區(qū)農(nóng)業(yè)研究，"2020，"38（5）："226-235.CAO"Y，"WANG"J，"WANG"S"P，"HUANG"Y，"ZHANG"L，"LI"Y"L."Temporal"and"spatial"variation"analysis"of"water"and"irrigation"requirement"index"of"rice"in"central"Yunnan"province"under"climate"change[J]."Agricultural"Research"in"the"Arid"Areas，"2020，"38（5）："226-235."（in"Chinese）
• 鄒海平，"張京紅，"李偉光，"陳小敏，"白蕤，"呂潤."海南島水稻需水量與缺水量的時(shí)空變化特征[J]."熱帶生物學(xué)報(bào)，"2023，"14（5）："569-576.ZOU"H"P，"ZHANG"J"H，"LI"W"G，"CHEN"X"M，"BAI"R，"LYU"R."Spatial"and"temporal"variations"of"water"requirement"and"water"deficit"of"rice"in"Hainan"island[J]."Journal"of"Tropical"Biology，"2023，"14（5）："569-576."（in"Chinese）

[13]"朱士江，"李凱凱，"徐文，"丁學(xué)軍，"張國紅，"朱瑾，"貴樹彪，"張禮杰."氣候變化對三湖連江灌區(qū)水稻生育期和灌溉需水量變化的影響[J]."安徽農(nóng)業(yè)大學(xué)學(xué)報(bào)，"2023，"50（2）："191-198.ZHU"S"J，"LI"K"K，"XU"W，"DING"X"J，"ZHANG"G"H，"ZHU"J，"GUI"S"B，"ZHANG"L"J."Effects"of"climate"change"on"rice"growth"period"and"variation"of"irrigation"water"requirement"in"Sanhu"Lianjiang"irrigation"district[J]."Journal"of"Anhui"Agricultural"University，"2023，"50（2）："191-198."（in"Chinese）

[14]"武夢晗，"王弋."新疆莎車灌區(qū)農(nóng)作物需水量和產(chǎn)量特征及其影響因素[J]."生態(tài)學(xué)雜志，"2023，"42（8）："1858-1868.WU"M"H，"WANG"Y."Characteristics"and"influencing"factors"of"crop"water"demand"and"yield"in"Shache"irrigation"area"of"Xinjiang[J]."Chinese"Journal"of"Ecology，"2023，"42（8）："1858-"1868."（in"Chinese）

[15]"王潔，"姚懷柱，"梁森，"侯新蕊，"張脈然，"胡家珍，"楊士紅."江蘇省冬小麥需水量時(shí)空變化特征及氣候影響因素分析[J]."節(jié)水灌溉，"2024（5）："107-114.WANG"J，"YAO"H"Z，"LIANG"S，"HOU"X"R，"ZHANG"M"R，"HU"J"Z，"YANG"S"H."Characterization"of"spatial"and"temporal"variations"in"water"demand"of"winter"wheat"in"Jiangsu"province"and"analysis"of"climate"influencing"factors[J]."Water"Saving"Irrigation，"2024（5）："107-114."（in"Chinese）

[16]"姬興杰，"成林，"朱業(yè)玉，"宋妮."河南省冬小麥需水量和缺水量的時(shí)空格局[J]."生態(tài)學(xué)雜志，"2014，"33（12）："3268-3277.JI"X"J，"CHENG"L，"ZHU"Y"Y，"SONG"N."Spatial"and"temporal"characteristics"of"water"requirement"and"water"deficit"of"winter"wheat"in"Henan"province[J]."Chinese"Journal"of"Ecology，"2014，"33（12）："3268-3277."（in"Chinese）

[17]"宋進(jìn)喜，"高雋清，"李曉鑫，"毛瑞晨，"史阿瑩."近20年來黃土高原蒸散發(fā)變化規(guī)律及其驅(qū)動(dòng)因素[J]."西北大學(xué)學(xué)報(bào)（自然科學(xué)版），"2023，"53（6）："974-990.SONG"J"X，"GAO"J"Q，"LI"X"X，"MAO"R"C，"SHI"A"Y."Changes"of"evapotranspiration"and"its"driving"factors"in"the"Loess"Plateau"in"recent"20"years[J]."Journal"of"Northwest"University"（Natural"Science"Edition），"2023，"53（6）："974-990."（in"Chinese）

[18]"譚嬌，"丁建麗，"董煜，"楊愛霞，"張喆."新疆艾比湖綠洲潛在蒸散量年代際變化特征[J]."農(nóng)業(yè)工程學(xué)報(bào)，"2017，"33（5）："143-148.TAN"J，"DING"J"L，"DONG"Y，"YANG"A"X，"ZHANG"Z."Decadal"variation"of"potential"evapotranspiration"in"Ebinur"lake"oasis"of"Xinjiang[J]."Transactions"of"the"Chinese"Society"of"Agricultural"Engineering，"2017，"33（5）："143-148."（in"Chinese）

[19]"劉泉杉，"吳宗俊，"崔寧博，"朱世丹，"鄭順生."云貴高原參考作物蒸散量時(shí)空特征及成因分析[J]."排灌機(jī)械工程學(xué)報(bào)，"2022，"40（3）："302-310.LIU"Q"S，"WU"Z"J，"CUI"N"B，"ZHU"S"D，"ZHENG"S"S."Spatial-temporal"distribution"characteristics"and"attribution"analysis"of"reference"crop"evapotranspiration"in"Yunnan-Kweichow"Plateau[J]."Journal"of"Drainage"and"Irrigation"Machinery"Engineering，"2022，"40（3）："302-310."（in"Chinese）

[20]"ALLEN"R"G，"PEREIRA"L"S，"RAES"D，"SMITH"M."Crop"evapotranspiration："guidelines"for"computing"crop"water"requirements"-"FAO"irrigation"and"drainage"paper"56[M]."Rome："Food"and"Agriculture"Organization"of"the"United"Nations，"1998："37-58．

[21]"BROUWER"C，"HEIBLOEM"M."Irrigation"water"management："irrigation"water"needs，"irrigation"water"management"training"manual"No."3[M]."Rome："Food"and"Agriculture"Organization"of"the"United"Nations，"1986："63-68.

[22]"SMITH"M."CROPWAT："a"computer"program"for"irrigation"planning"and"management"-"FAO"irrigation"and"drainage"paper"No."46）[M]."Rome："Food"and"Agriculture"Organization"of"the"United"Nations，"1992："20-21.

[23]"王明田，"曲輝輝，"楊曉光，"張曉煜，"李茂松."基于降水保證指數(shù)的四川省種植制度優(yōu)化研究[J]."西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(xué)版），"2012，"40（10）："82-92.WANG"M"T，"QU"H"H，"YANG"X"G，"ZHANG"X"Y，"LI"M"S."Optimized"cropping"systems"based"on"ensure"index"of"precipitation"in"Sichuan"province[J]."Journal"of"Northwest"A"amp;"F"University"（Natural"Science"Edition），"2012，"40（10）："82-92."（in"Chinese）

[24]"鄒海平，"張京紅，"陳小敏，"劉少軍，"李偉光."海南島農(nóng)業(yè)氣候資源的時(shí)空變化特征[J]."中國農(nóng)業(yè)氣象，"2015，"36（4）："417-427.ZOU"H"P，"ZHANG"J"H，"CHEN"X"M，"LIU"S"J，"LI"W"G."Spatiotemporal"change"characteristics"of"agricultural"climate"resources"in"Hainan"island[J]."Chinese"Journal"of"Agrome teorology，"2015，"36（4）："417-427."（in"Chinese）

[25]"陳玉民，"郭國雙，"王廣興，"康紹忠，"羅懷彬，"張大中."中國主要作物需水量與灌溉[M]."北京："水利電力出版社，"1995："115.CHEN"Y"M，"GUO"G"S，"WANG"G"X，"KANG"S"Z，"LUO"H"B，"ZHANG"D"Z."Main"crop"water"requirement"and"irrigation"of"China[M]."Beijing："Water"Resources"and"Electric"Power"Press，"1995："115."（in"Chinese）

[26]"陳小敏，"鄒海平，"張京紅，"劉少軍，"蔡大鑫."海南農(nóng)業(yè)氣候資源與主要作物區(qū)劃[M]."北京："氣象出版社，"2020："13-24.CHEN"X"M，"ZOU"H"P，"ZHANG"J"H，"LIU"S"J，"CAI"D"X."Agroclimatic"resources"and"main"crop"regionalization"in"Hainan[M]."Beijing："China"Meteorological"Press，"2020："13-"24."（in"Chinese）