蔡普默，趙夢婷，宋蘊(yùn)哲，孟李濤，季清娥，林 嘉，李焰焰，聶傳朋

（1.武夷學(xué)院茶與食品學(xué)院園藝系 福建武夷山 354300； 2.生物農(nóng)藥與化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室·閩臺(tái)作物有害生物生態(tài)防控國家重點(diǎn)實(shí)驗(yàn)室·福建農(nóng)林大學(xué)作物病蟲生物防治研究所 福州 350002）

以大氣CO等溫室氣體濃度升高而造成的環(huán)境溫度升高現(xiàn)象為主要特征的“全球氣候變化”已成為各國政府、科研工作者、新聞媒體及公眾最關(guān)注的環(huán)境問題。政府間氣候變化專門委員會(huì)（Intergovernmental Panel on Climate Change，IPCC）指出，從1880 年到2012 年，全球平均氣溫上升了0.72 ℃，平均上升速度為每10 年上升0.12 ℃。我國是全球氣候變暖最顯著的國家之一，1951 年至2021 年，我國平均氣溫每10 年的增速為0.26 ℃，高于同期全球平均升溫水平（0.15 ℃/10 年）。

昆蟲作為數(shù)量、種類最多的動(dòng)物，同樣也受到氣候變暖的影響。國內(nèi)外眾多學(xué)者的研究已明確，氣候變暖使害蟲的發(fā)育速度加快、世代增多、危害時(shí)間提前，導(dǎo)致害蟲暴發(fā)的強(qiáng)度和頻率增加，增加了害蟲防控的難度。此外，害蟲為適應(yīng)全球氣候變暖，向高緯度和高海拔地區(qū)遷移擴(kuò)散，導(dǎo)致其發(fā)生范圍擴(kuò)大。害蟲、寄主植物、天敵等不同營養(yǎng)級生物對全球氣候變暖的響應(yīng)不同，改變了植物-害蟲-天敵三者發(fā)生時(shí)間與空間格局，原有的三層營養(yǎng)級之間的內(nèi)在關(guān)系和平衡格局也發(fā)生變化，進(jìn)而引起寄主植物的抗性和天敵所發(fā)揮的作用發(fā)生改變，致使一些害蟲發(fā)生嚴(yán)重，而一些害蟲發(fā)生減緩甚至滅絕。由此可見，揭示全球氣候變暖對昆蟲影響的生態(tài)學(xué)機(jī)制，對氣候變暖大背景下的農(nóng)林重大害蟲可持續(xù)綜合防控新技術(shù)和新方法的研究具有重大意義。

瓜 實(shí) 蠅[（Coquillett）]屬雙翅目（Diptera）實(shí)蠅科（Tephritidae）鏃實(shí)蠅屬（），是我國進(jìn)境植物檢疫性害蟲和重點(diǎn)管控的外來入侵生物。瓜實(shí)蠅是一種重要的果蔬害蟲，嚴(yán)重危害苦瓜、黃瓜、西瓜、南瓜等100 余種寄主作物。目前，瓜實(shí)蠅廣泛分布于亞熱帶30 多個(gè)國家和地區(qū)。據(jù)報(bào)道，瓜實(shí)蠅最早于1985 年在我國深圳口岸由香港輸入內(nèi)地的白瓜中截獲，后在深圳、上海、?？凇⒗ッ骱徒T等口岸被多次截獲。目前，瓜實(shí)蠅在我國主要分布于福建、廣東、廣西、江蘇、海南、云南、四川等地，對當(dāng)?shù)氐墓呱a(chǎn)造成較大威脅。研究表明，未來的氣候變化會(huì)促使該害蟲在國內(nèi)的適生范圍北移，且適生區(qū)面積增加，危害程度加重。孫宏禹等基于@RISK 軟件和隨機(jī)模擬方法對瓜實(shí)蠅給我國苦瓜產(chǎn)業(yè)造成的經(jīng)濟(jì)損失進(jìn)行評估，結(jié)果表明，在防治的情況下，瓜實(shí)蠅每年給我國苦瓜產(chǎn)業(yè)可能造成的直接經(jīng)濟(jì)損失達(dá)141.6 億元；而在不防治的情況下，損失則高達(dá)234.8 億元。

當(dāng)前實(shí)蠅的防治手段主要有農(nóng)藝防治、化學(xué)防治、誘殺成蟲和生物防治等，實(shí)蠅物候?qū)W是科學(xué)合理制定及實(shí)施不同防治方案的重要參考。在全球氣候變化的背景下，實(shí)蠅物候可能會(huì)根據(jù)氣候變化做出適應(yīng)性的改變，現(xiàn)有防控措施的防治效果也會(huì)受到相應(yīng)的影響。因此，探究實(shí)蠅對氣候變暖的物候響應(yīng)對制定科學(xué)、有效的防控措施具有指導(dǎo)意義。由于氣候變暖的長期性和緩慢性，使得這項(xiàng)研究工作通過田間實(shí)地實(shí)時(shí)調(diào)查的難度和工作量增大，可行性不高。在地中海地區(qū)，Gordo 等基于橄欖實(shí)蠅50 年的監(jiān)測數(shù)據(jù)和同期的氣候數(shù)據(jù)，通過多元線性回歸分析發(fā)現(xiàn)在1975 年后橄欖實(shí)蠅的始見期顯著提前，而1975 年前始見期則延遲，但橄欖實(shí)蠅的始見期整體呈提前趨勢。借鑒前人的研究結(jié)果和方法，筆者收集了瓜實(shí)蠅歷年來在國內(nèi)的發(fā)生情況與發(fā)生地區(qū)的溫度等資料，通過分析田間成蟲始見期和種群高峰期年際變化情況、發(fā)生地區(qū)溫度變化情況及二者的線性關(guān)系，實(shí)證研究了瓜實(shí)蠅對氣候變暖的田間物候響應(yīng)。該研究旨在深入了解氣候變暖對害蟲物候?qū)W的影響，為優(yōu)化瓜實(shí)蠅的防控方案提供科學(xué)依據(jù)。

1 材料與方法

1.1 瓜實(shí)蠅發(fā)生數(shù)據(jù)收集與處理

中國知網(wǎng)數(shù)據(jù)庫是目前收錄最廣泛、最全面的中文期刊數(shù)據(jù)庫。因此，筆者研究中的瓜實(shí)蠅發(fā)生數(shù)據(jù)來源于中國知網(wǎng)數(shù)據(jù)庫的文獻(xiàn)記錄。在檢索知網(wǎng)數(shù)據(jù)庫時(shí)，以物種拉丁學(xué)名或中文名為主檢索詞，如瓜實(shí)蠅、針蜂、瓜蛆、、；再分別結(jié)合其他檢索詞如種群動(dòng)態(tài)、高峰、越冬和發(fā)生。為保證檢索結(jié)果的全面性，多次使用以上詞條的不同組合進(jìn)行檢索，將所需的瓜實(shí)蠅發(fā)生數(shù)據(jù)從文獻(xiàn)中提取出，匯總整理后，建立瓜實(shí)蠅發(fā)生數(shù)據(jù)庫。瓜實(shí)蠅為多化性昆蟲，世代重疊嚴(yán)重，以老熟幼蟲化蛹越冬。筆者將田間首次監(jiān)測到成蟲的日期作為始見期，以田間誘集蟲量達(dá)最大值的日期為種群高峰期。由于國內(nèi)文獻(xiàn)記錄的瓜實(shí)蠅始見期和種群高峰期皆始于2001 年，因此本研究瓜實(shí)蠅物候數(shù)據(jù)研究周期設(shè)定為2000 至2020 年。

為了使物候指標(biāo)數(shù)據(jù)標(biāo)準(zhǔn)化，在文獻(xiàn)中常使用一些較為模糊的時(shí)間概念，如7 月下旬、8 月中旬等，以取中間值的方式將模糊的時(shí)間量化。例如，上旬為每月前10 d，取5 日；中旬為每月11—20日，取15 日；下旬為每月最后10 d 或9 d，取25日。若時(shí)間跨度較廣取中間值，如4 月20 日至5月6 日，取中間值4 月28 日。在部分文獻(xiàn)中，記錄了同一地區(qū)不同觀測點(diǎn)對瓜實(shí)蠅的觀測值，則取這些值的均值作為該地區(qū)的觀測值。筆者的研究將瓜實(shí)蠅物候數(shù)據(jù)量化為數(shù)字，以發(fā)生時(shí)間與當(dāng)年1月1 日之間相差的天數(shù)表示。例如1990 年在廣州地區(qū)記錄的瓜實(shí)蠅種群高峰期為10 月1 日，該日期與當(dāng)年1 月1 日之間相差的天數(shù)為273 d，則1990 年廣州地區(qū)瓜實(shí)蠅的種群高峰期以273 表示。

1.2 氣象數(shù)據(jù)獲取與分析

筆者所使用的氣象數(shù)據(jù)均來自于中國氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn/），數(shù)據(jù)集為中國國際交換站地面氣候標(biāo)準(zhǔn)值日值數(shù)據(jù)集、月值數(shù)據(jù)集及年值數(shù)據(jù)集。根據(jù)匯總的瓜實(shí)蠅發(fā)生年限和地區(qū)，選擇瓜實(shí)蠅發(fā)生較為集中的區(qū)域，將發(fā)生區(qū)域劃分為華南（廣東省、廣西壯族自治區(qū)和海南?。?、華東（江蘇省、浙江省、福建省、江西?。⑷A中（河南省、湖南省和湖北?。?、西南地區(qū)（貴州省、四川省和云南?。┘叭珖秶ｋS后，計(jì)算各地區(qū)的年平均溫度，并通過線性回歸分析各地區(qū)氣溫隨年份的變化趨勢，是氣候變暖還是變冷，這種變化的趨勢是否顯著。

1.3 數(shù)據(jù)處理和分析

分別對華南、華東、華中、西南地區(qū)及全國范圍內(nèi)瓜實(shí)蠅成蟲的始見期和種群高峰期進(jìn)行線性回歸分析，得出各個(gè)地區(qū)幾個(gè)物候指標(biāo)的年際變化情況。其中，回歸系數(shù)為負(fù)值表明物候期在年際間有提前趨勢，若為正值則說明有延遲趨勢，且回歸關(guān)系顯著則說明提前或延遲趨勢顯著，反之則表示趨勢不顯著。再用同樣的方法對不同地區(qū)瓜實(shí)蠅物候期與涉及地區(qū)年均溫進(jìn)行線性回歸分析，得出物候期變化與氣候變暖之間的關(guān)系。數(shù)據(jù)分析使用SPSS 軟件（Windows 20.0 版）。

2 結(jié)果與分析

2.1 基于中國知網(wǎng)瓜實(shí)蠅在我國的物候記錄情況

經(jīng)檢索并剔除不相關(guān)文獻(xiàn)后，共得記錄有瓜實(shí)蠅物候數(shù)據(jù)的相關(guān)中文文獻(xiàn)54 篇，文獻(xiàn)主要來源于《植物保護(hù)學(xué)報(bào)》《植物保護(hù)》《環(huán)境昆蟲學(xué)報(bào)》《中國植保導(dǎo)刊》《植物檢疫》《廣東農(nóng)業(yè)科學(xué)》《湖南農(nóng)業(yè)科學(xué)》《安徽農(nóng)業(yè)科學(xué)》《西南農(nóng)業(yè)學(xué)報(bào)》等期刊及福建農(nóng)林大學(xué)、廣西大學(xué)等高校碩士和博士畢業(yè)論文，經(jīng)整理后獲瓜實(shí)蠅物候記錄情況見表1。

表1 瓜實(shí)蠅在我國的物候記錄情況

2.2 我國瓜實(shí)蠅發(fā)生危害區(qū)氣候變化情況

自1990 年以來，華東地區(qū)（除山東省和安徽?。⑷A南地區(qū)、華中地區(qū)、西南地區(qū)（除西藏地區(qū)）及全國范圍的年均溫均呈上升趨勢。其中，華東地區(qū)升溫趨勢最為明顯，1990 年該地區(qū)的年平均氣溫為17.748 ℃，2020 年為19.244 ℃，30 年間上升了1.496 ℃，年平均升溫速率達(dá)（0.035±0.006）℃（圖1-d）。華南地區(qū)的年平均升溫速率為（0.018±0.007）℃（圖1-a），華中地區(qū)的年平均升溫速率為（0.023±0.007）℃（圖1-b），而西南地區(qū)的年平均升溫速率為（0.026±0.006）℃（圖1-c）。在這30 年間，我國年均氣溫以（0.026±0.006）℃的速率顯著上升（圖1-e）。

圖1 1990—2020 年間我國瓜實(shí)蠅物候記錄區(qū)的年均溫變化情況

2.3 瓜實(shí)蠅成蟲田間物候變化情況

2.3.1 我國不同地區(qū)瓜實(shí)蠅成蟲田間始見期變化情況 從圖2-a～c 可以看出，華南、華中、西南地區(qū)瓜實(shí)蠅歷年平均始見期分別為4 月28 日、5 月24日和4 月8 日。通過一元線性回歸分析，發(fā)現(xiàn)這3個(gè)地區(qū)瓜實(shí)蠅始見期和年份皆呈正回歸關(guān)系，說明年際間瓜實(shí)蠅始見期在這些地區(qū)延遲了，僅華中地區(qū)延遲趨勢達(dá)顯著水平。華南、華中、西南地區(qū)瓜實(shí)蠅成蟲始見期平均每年分別延遲3.489、5.270、1.987 d。由圖2-d 可以看出，華東地區(qū)瓜實(shí)蠅歷年平均始見期為5 月7 日，該地區(qū)始見期和年份呈負(fù)回歸關(guān)系，說明年際間華東地區(qū)的瓜實(shí)蠅成蟲始見期提前了，但未達(dá)顯著水平，平均每年約提前3.251 d。從圖2-e 看出，全國范圍瓜實(shí)蠅歷年平均始見期為5 月4 日，且始見期隨年份的推移顯著延遲，平均每年推遲約3.054 d。

圖2 我國瓜實(shí)蠅成蟲田間始見期年際變化情況

2.3.2 我國不同地區(qū)瓜實(shí)蠅成蟲田間種群高峰期變化情況從圖3-a～d 可以看出，華南、華中、西南、華東地區(qū)瓜實(shí)蠅歷年平均種群高峰期分別為7月15日、7月18日、7月18日和8月8日。通過一元線性回歸分析，發(fā)現(xiàn)這4個(gè)地區(qū)瓜實(shí)蠅成蟲種群高峰期和年份皆呈正回歸關(guān)系，說明年際間瓜實(shí)蠅成蟲種群高峰期在這些地區(qū)延遲了，且皆未達(dá)顯著水平。華南、華中、西南和華東地區(qū)瓜實(shí)蠅種群高峰期平均每年分別延遲4.371、3.544、5.184、0.036d。從圖4-e可以看出，全國范圍瓜實(shí)蠅歷年平均種群高峰期為7月26日，且年際間瓜實(shí)蠅種群高峰期在我國顯著延遲了，平均每年約延遲3.560d。

圖3 我國瓜實(shí)蠅成蟲田間種群高峰期年際變化情況

2.4瓜實(shí)蠅成蟲田間物候?qū)夂蜃兣捻憫?yīng)

2.4.1瓜實(shí)蠅成蟲田間始見期對氣候變暖的響應(yīng)通過對4個(gè)地區(qū)瓜實(shí)蠅成蟲始見期和對應(yīng)年份的年平均溫度進(jìn)行回歸分析，發(fā)現(xiàn)華南和華中地區(qū)瓜實(shí)蠅始見期和年均溫呈正回歸關(guān)系，說明瓜實(shí)蠅始見期在這2 個(gè)地區(qū)隨年均溫提升呈現(xiàn)延遲的趨勢，但未達(dá)顯著水平（圖4-a～b）。若華南和華中地區(qū)年均溫提高0.1℃，當(dāng)?shù)毓蠈?shí)蠅始見期分別推遲2.242 d 和6.388d。此外，西南和華東地區(qū)瓜實(shí)蠅始見期與年均溫呈負(fù)回歸關(guān)系，表示隨年均溫增加，這2 個(gè)地區(qū)的瓜實(shí)蠅始見期呈提前趨勢，但也未達(dá)顯著水平（圖4-c～d）。若西南和華東地區(qū)年均溫提高0.1℃，當(dāng)?shù)毓蠈?shí)蠅始見期分別提前1.021d 和3.082 d。從圖5 e 可以看出，瓜實(shí)蠅始見期在全國范圍內(nèi)與年均溫呈顯著正回歸關(guān)系，表明瓜實(shí)蠅始見期在我國隨年均溫增加呈顯著延遲趨勢。若全國年均溫提升0.1℃，瓜實(shí)蠅成蟲始見期推遲4.913d。

2.4.2瓜實(shí)蠅成蟲田間種群高峰期對氣候變暖的響應(yīng)由圖5-a～b可得，華南和華中地區(qū)瓜實(shí)蠅種群高峰期與年均溫呈正回歸關(guān)系，即瓜實(shí)蠅種群高峰期在這2個(gè)地區(qū)隨年均溫提高呈延遲的趨勢，但華南地區(qū)的變化未達(dá)顯著水平，而華中地區(qū)則達(dá)顯著水平。若華中和華南地區(qū)年均溫提高0.1℃，當(dāng)?shù)毓蠈?shí)蠅種群高峰期推遲8.723d 和0.630d。隨年均溫提高，西南和華東地區(qū)的瓜實(shí)蠅種群高峰期和年均溫呈負(fù)回歸關(guān)系，即這2個(gè)地區(qū)瓜實(shí)蠅種群高峰表現(xiàn)出不顯著的提前趨勢（圖5-c～d）。若華東和西南地區(qū)年均溫提高0.1℃，當(dāng)?shù)毓蠈?shí)蠅成蟲種群高峰期分別提前2.492d 和4.188 d。從圖5-e 可以看出，瓜實(shí)蠅種群高峰期在全國范圍內(nèi)與年均溫呈負(fù)回歸關(guān)系，說明瓜實(shí)蠅種群高峰期在我國隨年均溫提升呈提前趨勢，但趨勢不顯著。若全國年均溫提高0.1℃，則瓜實(shí)蠅種群高峰期在我國提前2.230d。

圖5 我國瓜實(shí)蠅成蟲種群高峰期對氣候變暖的響應(yīng)

圖5 我國瓜實(shí)蠅成蟲種群高峰期對氣候變暖的響應(yīng)

3 討論與結(jié)論

年均溫顯著升高是氣候變暖最明顯的特征，自國內(nèi)有瓜實(shí)蠅物候數(shù)據(jù)記錄以來，我國及該害蟲危害的幾個(gè)重要區(qū)域的年均溫以不同速率上升。近30 年來，我國增溫趨勢明顯，氣溫年均增速為0.026 ℃，高于同期全球氣溫年均增速。由此可見，我國是全球氣候變暖特征最顯著的國家之一。昆蟲物候變化是其響應(yīng)氣候變化最明顯的現(xiàn)象之一，國內(nèi)外眾多研究證實(shí)絕大多數(shù)害蟲由于氣候變暖生長發(fā)育速率加快，羽化和首次飛行提前，導(dǎo)致危害期提前，并且每年代數(shù)增多，使得有害生物發(fā)生強(qiáng)度和暴發(fā)頻率增加，害蟲防治的難度亦增加。然而不同種類昆蟲或同種昆蟲不同性別對氣候變暖的響應(yīng)存在差異。在1976—1998 年間，紅襟粉蝶（）和紅蛺蝶（）田間始見期分別提前17.5、36.3 d。

基于文獻(xiàn)記載的瓜實(shí)蠅田間物候數(shù)據(jù)及對應(yīng)年份的溫度數(shù)據(jù)，筆者研究發(fā)現(xiàn)，隨著年份的推移，國內(nèi)瓜實(shí)蠅的田間種群始見期和高峰期顯著延遲，且種群高峰期比始見期平均每年多推遲0.506 d。由于瓜實(shí)蠅種群終見期的文獻(xiàn)記錄少，因此無法對瓜實(shí)蠅種群終末期進(jìn)行分析，進(jìn)而無法確定瓜實(shí)蠅在田間的整個(gè)發(fā)生期是延長或縮短。不過可以明確的是，瓜實(shí)蠅在田間始見成蟲活動(dòng)到種群增長至峰值的這段時(shí)間（可稱為種群增長階段）延長，且總體上推遲了。通常，每年春季3—4 月份氣溫達(dá)10 ℃以上，越冬瓜實(shí)蠅蛹便開始羽化成成蟲，7～14 d后成蟲便開始交配，在寄主作物上產(chǎn)卵繁衍后代，即瓜實(shí)蠅在羽化7 d 后內(nèi)便可開始為害寄主植物，瓜實(shí)蠅始見期至種群高峰期的這一階段是危害瓜果的重要時(shí)期。因此，筆者認(rèn)為近20 年來國內(nèi)瓜實(shí)蠅在田間的危害期有延長且加重的態(tài)勢。

此外，瓜實(shí)蠅在我國不同地區(qū)田間物候的年際變化情況存在差異，比如，華南、華中、西南地區(qū)瓜實(shí)蠅始見期平均每年分別延遲3.489、5.270、1.987 d，而華東地區(qū)的瓜實(shí)蠅始見期平均每年約提前3.251 d；華南、華中、西南和華東地區(qū)瓜實(shí)蠅種群高峰期平均每年分別延遲4.371、3.544、5.184、0.036 d。由此可見，華南、華中和西南地區(qū)瓜實(shí)蠅種群增長階段整體呈現(xiàn)不同程度的推遲，而華中地區(qū)瓜實(shí)蠅種群增長時(shí)期縮短。原因可能是國內(nèi)不同地區(qū)影響瓜實(shí)蠅生長發(fā)育及種群動(dòng)態(tài)的生物條件（食物資源、天敵資源等）和非生物條件（溫度、光照、濕度等）不同。這也是在我國年均溫逐年升高的背景下，瓜實(shí)蠅種群高峰期隨年際變化情況與隨溫度升高變化情況不一致的主要因素。例如，不同地區(qū)種植的主要瓜果種類和品種不同或瓜實(shí)蠅在不同地區(qū)偏好危害的作物不同，不同寄主瓜果的物候期存在差異，進(jìn)而也影響瓜實(shí)蠅的田間發(fā)生情況。又比如，在海南省儋州市，李磊等發(fā)現(xiàn)瓜實(shí)蠅在苦瓜園發(fā)生高峰期集中于8 月，而瓜實(shí)蠅在混種有杧果、黃皮、人心果及零星的楊桃和蓮霧的果園中于7 月達(dá)到種群高峰期。

氣候變暖對昆蟲發(fā)生期和種群動(dòng)態(tài)造成直接影響的生態(tài)學(xué)基礎(chǔ)是溫度升高改變了昆蟲發(fā)育歷期、存活率及繁殖能力等重要生活史參數(shù)。已有研究表明氣候變暖加速昆蟲各蟲態(tài)的發(fā)育，使得物候期提前，如始見期、遷飛期和種群高峰期。但也有例外，如Gordo 等發(fā)現(xiàn)與瓜實(shí)蠅同屬實(shí)蠅科的橄欖實(shí)蠅始見期于1975 年前呈延遲的態(tài)勢，隨后則提前，但50 年間該害蟲初始發(fā)生期總體上提前，西方蜜蜂（）和菜粉蝶（）的始見期也表現(xiàn)出相同的氣候變化響應(yīng)模式。筆者研究發(fā)現(xiàn)隨著年均溫升高，瓜實(shí)蠅在華南和華中地區(qū)的始見期和種群高峰期也出現(xiàn)延遲趨勢。Dingemanse 等認(rèn)為昆蟲物候變化是由溫度和影響其幼蟲發(fā)育或成蟲存活的其他環(huán)境因素的改變而引起的綜合效應(yīng)，且除溫度之外的環(huán)境因子可能會(huì)削弱由溫度提升造成的物候提前，以致于昆蟲物候延遲。此外，在華東和西南地區(qū)瓜實(shí)蠅的始見期和種群高峰期皆提前，且瓜實(shí)蠅在全國范圍內(nèi)的始見期顯著延遲了，但種群高峰期提前。由此可見，瓜實(shí)蠅在不同地區(qū)對氣候變暖的物候響應(yīng)模式也不相同。

不同地區(qū)瓜實(shí)蠅對氣候變暖的物候響應(yīng)差異還體現(xiàn)在相同年增溫率下不同地區(qū)的瓜實(shí)蠅物候期改變速率不一致。例如，各地區(qū)的年均溫提高0.1 ℃，瓜實(shí)蠅在華東地區(qū)始見期提前3.082 d 且種群高峰期提前2.492 d，說明在該地區(qū)瓜實(shí)蠅發(fā)生提前且種群增長期延長；在華南地區(qū)，瓜實(shí)蠅始見期延遲2.242 d 且種群高峰期延遲0.630 d，表明其在華南地區(qū)發(fā)生延遲且種群增長期縮短；在華中地區(qū)，瓜實(shí)蠅始見期延遲6.388 d 且種群高峰期延遲8.723 d，說明該害蟲在該地區(qū)發(fā)生延遲且種群增長期延長；瓜實(shí)蠅在西南地區(qū)始見期提前1.021 d 且種群高峰期提前4.188 d，說明瓜實(shí)蠅在該地區(qū)發(fā)生提前且種群增長期縮短；在我國瓜實(shí)蠅發(fā)生區(qū)，其始見期延遲4.913 d 且種群高峰期提前2.230 d，說明瓜實(shí)蠅在我國發(fā)生延遲且種群增長期縮短。綜上，氣候變暖對我國瓜實(shí)蠅物候產(chǎn)生的影響具有空間異質(zhì)性，與前人的研究結(jié)果一致。Wu 等發(fā)現(xiàn)我國新疆南部近60 年的增溫速率高于新疆北部，使得桃蚜（）、棉蚜（）和麥長管蚜（）在我國新疆南部的初見期和遷移期早于新疆北部，他們認(rèn)為氣候變暖對蚜蟲物候的影響具有空間異質(zhì)性。

全球氣候變化對昆蟲的作用是長期的生態(tài)效應(yīng)，除了能改變瓜實(shí)蠅物候期之外，還可能擴(kuò)大害蟲地理分布區(qū)，增加向兩極或高海拔地區(qū)擴(kuò)散的機(jī)會(huì)，協(xié)助其向新領(lǐng)地入侵定殖，擴(kuò)大害蟲的寄主范圍。前人的研究表明，未來氣候變暖促使瓜實(shí)蠅在我國的適生區(qū)向北擴(kuò)散，且適生面積擴(kuò)大，潛在的危害更為嚴(yán)重。此外，氣候變暖可能通過食物鏈擾亂瓜實(shí)蠅-寄主植物-天敵的種間物候同步性，有可能會(huì)產(chǎn)生以下結(jié)果：氣候變暖降低瓜實(shí)蠅和其對應(yīng)天敵的物候匹配性，瓜實(shí)蠅失去天敵的制約而暴發(fā)成災(zāi)；氣候變暖可能加強(qiáng)相應(yīng)天敵的控害效能而進(jìn)一步限制瓜實(shí)蠅種群增長；氣候變暖削弱瓜實(shí)蠅與寄主植物的同步性，若瓜實(shí)蠅對與寄主同步性削弱有很強(qiáng)的適應(yīng)能力，通過調(diào)節(jié)其田間種群物候發(fā)生期來恢復(fù)與寄主植物的同步性，這可能是瓜實(shí)蠅物候期隨氣候變暖改變的原因之一；若瓜實(shí)蠅對與寄主同步性削弱的適應(yīng)能力不強(qiáng)的話，有可能致使瓜實(shí)蠅轉(zhuǎn)移寄主危害，擴(kuò)大瓜實(shí)蠅的寄主范圍；氣候變暖通過影響昆蟲的取食總量和取食速率，進(jìn)而改變瓜實(shí)蠅的危害程度；氣候變暖可能會(huì)影響害蟲的發(fā)育速率、生殖能力和存活率，使得其田間種群密度發(fā)生改變。上述都是在全球氣候變暖背景下防控瓜實(shí)蠅所需考慮的，也是未來昆蟲應(yīng)對氣候變暖響應(yīng)研究的重點(diǎn)。

筆者基于文獻(xiàn)中的瓜實(shí)蠅發(fā)生數(shù)據(jù)和同期氣候數(shù)據(jù)建立回歸模型進(jìn)行的分析，有助于探究氣候變暖對昆蟲的長期效應(yīng)。然而筆者的研究具一定的局限性，例如，不同文獻(xiàn)中瓜實(shí)蠅田間調(diào)查方法不同而造成系統(tǒng)誤差大，數(shù)據(jù)質(zhì)量可能不高，如謝思偉等發(fā)現(xiàn)利用誘蠅酮和誘蠅醚實(shí)蠅引誘劑在苦瓜園中6 月上旬即可誘捕瓜實(shí)蠅成蟲，而黃板則需在7 月上旬。此外，某些歷史文獻(xiàn)提供的時(shí)間信息較模糊，時(shí)間的量化存在不確定性。因此，為減少文獻(xiàn)數(shù)據(jù)不確定性帶來的影響，必須仔細(xì)收集、標(biāo)準(zhǔn)化處理和分析數(shù)據(jù)。全球氣候變暖的特征不僅是均溫提升，還包括最高溫升高、日溫度積累、晝夜溫差變化，以及伴隨著如CO濃度升高、O濃度提升等其他因素。筆者的研究僅考慮瓜實(shí)蠅對年均溫提升的物候響應(yīng)，氣候變暖的其他特征引起昆蟲的物候響應(yīng)仍缺乏有效探索。因此，今后需綜合分析多個(gè)環(huán)境氣候因素對昆蟲的影響，才能準(zhǔn)確掌握未來全球氣候變化背景下昆蟲災(zāi)變規(guī)律及對害蟲的發(fā)生提出準(zhǔn)確預(yù)警。盡管如此，筆者的研究將為今后更多其他生物類群和地理區(qū)域探索相似的科學(xué)問題提供參考。