尚廣春，李亮，馬春梅,2，弋雙文，趙琳，崔安寧

1. 南京大學(xué)地理與海洋科學(xué)學(xué)院，南京 210023

2. 江蘇省氣候變化協(xié)同創(chuàng)新中心，南京 210023

沉積物可以有效地反映地質(zhì)與歷史時(shí)期的氣候變化和沉積環(huán)境[1-4]。為了對(duì)未來(lái)的環(huán)境變化進(jìn)行預(yù)測(cè)，重建地質(zhì)時(shí)期尤其是全新世以來(lái)的環(huán)境演變變得十分重要[5-6]。過(guò)去全球變化（PAGES）的基本任務(wù)旨在研究過(guò)去氣候環(huán)境、生態(tài)過(guò)程和人類(lèi)活動(dòng)的長(zhǎng)期相互作用。湖泊水量是控制粒度分布的重要因素，粒度的粗細(xì)代表水動(dòng)力的強(qiáng)弱和湖泊水量的相對(duì)大小[7-9]；在黃土沉積物中，粗顆粒的含量和中值粒徑的大小可以指示歷史時(shí)期沙塵天氣發(fā)生的頻率和強(qiáng)度[10-11]；而在古洪水的判別中，粒度頻率分布和累積曲線(xiàn)也是很重要的判別方法[12-13]。隨著數(shù)學(xué)在沉積學(xué)中的應(yīng)用愈加廣泛和深入，利用因子分析法[14-15]、Weibull 分布函數(shù)擬合[16]、粒級(jí)-標(biāo)準(zhǔn)偏差[17-18]、端元模型法[19-23]等數(shù)學(xué)方法可以從多峰態(tài)的粒度頻率分布曲線(xiàn)中分離出單一粒度組分特征，提取敏感組分并對(duì)沉積物的來(lái)源進(jìn)行分析。Liu[23]等在蘇凱湖通過(guò)端元法得出4 個(gè)端元，取其中較有意義的3 個(gè)做對(duì)數(shù)運(yùn)算，得出的HI 指數(shù)與季風(fēng)有很好的對(duì)應(yīng)關(guān)系；Yu[22]等在端元分析法基礎(chǔ)上改進(jìn)的貝葉斯端元分析法，應(yīng)用于西部的沙漠和湖泊。但運(yùn)用端元分析法對(duì)泥炭地的探究未有實(shí)例，因此本文希望對(duì)此進(jìn)行研究。

過(guò)去2 000 a 的氣候變化由來(lái)已久，最早對(duì)過(guò)去2 000 a 的氣候研究起源于歐洲學(xué)者，Lamb[24]定義了羅馬暖期（Roman W.）、黑暗時(shí)代冷期（Dark Age Cold）、中世紀(jì)暖期（MWP）、小冰期（LIA）的大致起訖年代。但歐洲的氣候情況與我國(guó)不完全相同，因此，我國(guó)也有很多學(xué)者致力于研究近2 000 a 的氣候研究，最早的有竺可楨先生通過(guò)物候法對(duì)我國(guó)近五千年來(lái)氣候變化得出溫度曲線(xiàn)，劃分出歷朝的大致冷暖時(shí)期；葛全勝[25]等在前人研究基礎(chǔ)上得出，中國(guó)的溫度變化存在準(zhǔn)700 a 和200～300 a 的周期；降水變化則存在20～30 a 和準(zhǔn)70 a 的周期。鄭景云[26]等的研究表明中國(guó)東部降水同樣存在數(shù)百年的趨勢(shì)變化，西部中世紀(jì)暖期與小冰期不明顯；中世紀(jì)暖期時(shí)，中國(guó)東部的華北地區(qū)相對(duì)干旱，江南則相對(duì)濕潤(rùn)；而在小冰期，華北地區(qū)則相對(duì)濕潤(rùn)，且整個(gè)東部地區(qū)的降水變率增大。

本文基于江西玉華山泥炭鉆孔沉積物高分辨率AMS14C 測(cè)年的年代框架，結(jié)合粒度特征、沉積環(huán)境分析，運(yùn)用粒級(jí)-標(biāo)準(zhǔn)偏差提取的敏感組分和EMMA 端元分析法[20]的端元進(jìn)行比較，將歷史資料與石筍記錄等數(shù)據(jù)進(jìn)行對(duì)比分析，明確了粒度各端元組分的環(huán)境指示意義，揭示了研究區(qū)內(nèi)近2 000 a的環(huán)境變化及其控制因素。

1 區(qū)域地理背景

玉華山位于江西省中部（圖1），最高峰1 169 m，其水系整體呈北北東向分布。該區(qū)屬贛中南亞熱帶潮濕多雨季風(fēng)氣候區(qū)，四季交替明顯，全年冷暖氣團(tuán)出現(xiàn)頻繁，天氣復(fù)雜且多變。春季多雨，來(lái)自南方的暖濕氣流與北方的干冷氣流相遇，出現(xiàn)強(qiáng)降雨、大風(fēng)等強(qiáng)對(duì)流天氣，4 月降水集中；夏季高溫炎熱，6 月降水集中；秋季由于副熱帶高壓的影響，降水較少；冬季較寒冷，溫度下降明顯，并伴有雨雪等天氣。年平均氣溫18 ℃，夏季最高氣溫達(dá)到39.5 ℃，冬季最低氣溫為-9.2 ℃。年降水量達(dá)到1 500～2 000 mm，3—6 月份為雨季，年蒸發(fā)量為1 100～1 600 mm[27-28]。

圖1 研究區(qū)地形圖Fig.1 Map showing the location of study area

2 材料與方法

2.1 樣品采集

通過(guò)對(duì)玉華山的野外考察，研究組于2014 年用半刀式俄羅斯泥炭鉆鉆取了鉆孔，命名為YHS2（27°50′29.44″N、115°38′54.98″E），YHS2深度為3.5 m。

2.2 實(shí)驗(yàn)方法

本文選取了17 個(gè)YHS2 鉆孔的樣品挑選陸生植物殘?bào)w或提取孢粉濃縮物在美國(guó)Beta 實(shí)驗(yàn)室進(jìn)行AMS14C 年代測(cè)定。

YSH2 鉆孔共分得350 個(gè)樣品，粒度實(shí)驗(yàn)分析步驟以前人研究得出的黃土粒度前處理方法為基礎(chǔ)[29-30]，根據(jù)泥炭的特點(diǎn)做出適當(dāng)調(diào)整，具體步驟如下：取大約0.5 g 樣品置于200 mL 燒杯中，靜置一夜后，加10%雙氧水（去除樣品中有機(jī)質(zhì)）浸泡并加熱促使反應(yīng)，直到無(wú)氣泡產(chǎn)生為止（因泥炭地的樣品含有機(jī)質(zhì)較多，反應(yīng)時(shí)間較長(zhǎng)）。然后，加入10%的HCl 適量（去除樣品中碳酸鹽類(lèi)物質(zhì)），于電熱板上加熱煮沸，直到基本無(wú)氣泡產(chǎn)生為止。加入超純水浸泡樣品靜置一夜，重復(fù)兩遍，樣品水洗至中性，加入0.006 mol/L 的六偏磷酸鈉[（NaPO3）6] 試劑，超聲波振蕩器中震蕩10 min，使細(xì)顆粒物質(zhì)充分散開(kāi)，用Mastersizer2000 激光粒度儀進(jìn)行粒度檢測(cè)。

3 結(jié)果

3.1 剖面巖性特征

YHS2 鉆孔根據(jù)地層巖性差異自上而下劃分為14 個(gè)巖性單元（圖2），分別為：0～30 cm 為棕色沼澤土；30～70 cm 為灰黃色沼澤土；70～108 cm 為灰黑色泥炭土；108～130 cm 為灰黃色沼澤土；130～145 cm為灰黑色泥炭土；145～165 cm 為灰黃色沼澤；165～172 cm 為灰色沼澤土；172～217 cm 為灰色黏土；217～225 cm 為灰色粗砂；225～265 cm 為灰色黏土；265～304 cm 為褐色黏土；304～312 cm 為灰黑色泥炭土；312～323 cm 為灰色泥炭土；323～350 cm為褐色粉砂質(zhì)黏土。0～172 cm 為沼澤-泥炭沉積，172～304 cm 為黏土或砂，為湖相沉積，304～323 cm為泥炭沉積，323～350 cm 為褐色粉砂到黏土。

3.2 年代的校正

AMS14C 測(cè)年結(jié)果見(jiàn)表1[31]所示，為了更好地處理年代-深度關(guān)系，我們基于貝葉斯公式年代模型[32]，結(jié)合R 軟件為不同的年代控制點(diǎn)選擇合適的函數(shù)[33]，從而為YHS2 鉆孔建立準(zhǔn)確完整的年代框架。

3.3 粒度組成與參數(shù)特征

沉積物粒度組成的粒級(jí)分成黏土（＜4 μm）、細(xì)粉砂（4～16 μm）和中粗粉砂（16～64 μm）、砂（＞64 μm）4 個(gè)級(jí)別；粒度參數(shù)是綜合反映沉積物粒度特征及沉積環(huán)境的量化指標(biāo)，本文分析采用Fork、Ward 提出的平均粒徑（Mz）、分選系數(shù)（r1）、偏度（SK）、峰態(tài)（Ku）4 項(xiàng)粒度參數(shù)，結(jié)合平均粒徑輔以比較。

圖2 年代-深度模型Fig.2 Age-depth model

YSH2 鉆孔的粒度組成與參數(shù)見(jiàn)圖3 所示，平均粒徑6.57～3.18Φ，平均值5.77Φ；黏土含量6.46%～20.07%，平均值14.48%；細(xì)粉砂含量13.42%～52.46%，平均值34.36%；粗粉砂含量18.50%～48.85%，平均值34.67%；砂含量1.57%～61.60%，平均值16.46%，分選系數(shù)1.06～2.97，平均值1.99；偏度-0.30～0.38，平均值0.013；峰態(tài)0.70～1.58，平均值1.11。

3.4 敏感組分的提取

粒級(jí)-標(biāo)準(zhǔn)偏差法用于提取敏感組分、進(jìn)行古氣候變化的研究由來(lái)已久，對(duì)于追溯物源，探討季風(fēng)海流等提供重要信息[34-35]。本方法是依據(jù)每一粒級(jí)對(duì)應(yīng)含量的標(biāo)準(zhǔn)偏差變化來(lái)獲取環(huán)境敏感因子，它的基本思路是計(jì)算每一粒級(jí)在柱狀沉積物中的標(biāo)準(zhǔn)偏差值，將標(biāo)準(zhǔn)偏差值最大的粒級(jí)作為環(huán)境敏感因子，其敏感粒級(jí)往往是單一的。由此可以看出YHS2 鉆孔存在3 個(gè)敏感的組分，分別為10.473～16.535、29.264 ～72.951 、404.412 ～453.326 μm，占比最高的10.473～16.535 μm 組分是氣候的最敏感組分。

3.5 端元組分的提取

EMMA 端元模型法[20]能夠區(qū)分粒度的不同組分，在MATLAB 環(huán)境下運(yùn)行Paterson[21]等提供的Analysize 軟件包，然后通過(guò)綜合分析各自線(xiàn)性相關(guān)性、角度偏差和端元相關(guān)度確定了最佳端元數(shù)量。對(duì)粒度進(jìn)行計(jì)算，此地屬于泥炭地，多為粉砂與黏土，沉積物來(lái)源穩(wěn)定單一，基于上述討論，本文采用3 個(gè)端元，兩個(gè)單峰模式和有一個(gè)高低峰的模式，分別標(biāo)記為EM1、EM2、EM3。EM1 的峰值大約為13 μm、EM2 為32 μm，EM3 低峰為8 μm，而高峰為453 μm。

表 1 年代-深度模型[32]Table 1 List of Age-depth model[32]

圖3 粒度組分與參數(shù)圖Fig.3 Particle size composition and parameters

4 討論

4.1 年代及粒度組成參數(shù)的討論

YHS2 鉆孔雖然在最底部的測(cè)年達(dá)到7 ka，但是從沉積速率來(lái)看，在312 cm 以上沉積速率較快，為0.23 cm/a，312 cm 以下速率緩慢，其中312～322 cm沉 積 速 率 為0.004 cm/a，322～337 cm 沉 積 速 率為0.0087 cm/a，337～350 cm 沉積速率為0.0203 cm/a。我們認(rèn)為在312 cm 以下分辨率較低，與前部差距較大，可另做其他方面研究?；诖送茢?，本文主要對(duì)該鉆孔312 cm 以上部分的粒度進(jìn)行環(huán)境替代指標(biāo)變化分析，討論過(guò)去2 ka 期間的氣候環(huán)境變化。

玉華山泥炭位于山間盆地，我國(guó)西部的泥炭主要來(lái)源物質(zhì)以風(fēng)成為主，而東部以水成為主，且受到東亞季風(fēng)影響強(qiáng)烈[11,36]。物質(zhì)輸入粒度的粗細(xì)受到碎屑物質(zhì)來(lái)源和水位波動(dòng)相對(duì)變化的影響，暴雨和大暴雨會(huì)夾帶少量大顆粒，周邊的落石也會(huì)滾落而至，分選性較差。分選性和平均粒徑有一定的相關(guān)性，即較粗顆粒時(shí)，分選較差，較細(xì)顆粒時(shí)分選較好。偏態(tài)基本都在對(duì)稱(chēng)態(tài)左右波動(dòng)，沒(méi)有明顯的趨向。峰度的變化范圍基本都為0.9～1.4，峰態(tài)較寬?？偟膩?lái)說(shuō)是變細(xì)的過(guò)程，砂的含量在減少、粉砂的含量在增加。

4.2 端元組分的探討

懸浮微粒沉積過(guò)程與其輸導(dǎo)機(jī)制和沉積環(huán)境密切相關(guān)，湖泊沉積物粒度組成的控制因素主要有3 個(gè)：（1）碎屑物質(zhì)來(lái)源，（2）湖泊水位變化，（3）湖泊水動(dòng)力[37]。研究區(qū)位于沼澤泥炭，東亞季風(fēng)帶來(lái)大部分降水，山間盆地蓄水形成泥炭沼澤或是小型湖泊；同時(shí)，物質(zhì)輸入粒度的粗細(xì)受到碎屑物質(zhì)來(lái)源和水位波動(dòng)的相對(duì)變化的影響。

EM1 組分的眾數(shù)為13 μm（圖4），占比最大，對(duì)應(yīng)的敏感組分為10.473～16.535 μm，是對(duì)氣候干濕變化最為敏感的部分，屬于細(xì)粉砂，粒度特征分布為正態(tài)分布，與平均粒徑的變化也較為一致，此地是山間集水盆地，粒徑含量的變化主要反映氣候干濕對(duì)粒徑的影響，當(dāng)東亞季風(fēng)增強(qiáng)、降水增多、氣候濕潤(rùn)、侵蝕能力加強(qiáng)，沖入更多的水流易帶走的細(xì)粉砂粒徑的物質(zhì)；反過(guò)來(lái)沉積物中EM1 的含量增多，指示水動(dòng)力變強(qiáng)、降水增多、氣候濕潤(rùn)、東亞夏季風(fēng)增強(qiáng)。同理，EM1 減少，水動(dòng)力變?nèi)?、降水減少、氣候干燥。EM1 指代的是降水的變化，即氣候的干濕狀況的反映，與葛全勝重建的江南地區(qū)的干濕指數(shù)和董哥洞石筍的同位素記錄有很好的對(duì)應(yīng)關(guān)系。

EM2 組分的眾數(shù)為32 μm（圖4），對(duì)應(yīng)的偏差法的敏感組分為29.264～72.951 μm，屬于粗粉砂，單峰結(jié)構(gòu)反映較好的分選性，同此鉆孔的腐殖化度曲線(xiàn)有較好的對(duì)應(yīng)，腐殖化度[38-39]是在微生物作用下，通過(guò)生物與化學(xué)作用形成腐殖質(zhì)的過(guò)程，腐殖化度高代表成壤作用強(qiáng)，反之代表弱，因此EM2 可能是土壤的成壤改造作用，受環(huán)境的氧化還原影響，此地泥炭發(fā)育，沉積環(huán)境較穩(wěn)定，微生物作用持續(xù)而強(qiáng)烈，對(duì)沉積環(huán)境同樣起到更加適應(yīng)自己發(fā)展的改造作用，因此EM2 代表的是一定條件下土壤的成壤改造作用。

圖4 標(biāo)準(zhǔn)偏差法(上圖)和端元分析法(下圖)Fig.4 Standard deviation method （a） and end element analysis method （b）

EM3 是高低峰的形態(tài)，分為主峰和次峰，多峰結(jié)構(gòu)反映弱分選主峰的眾數(shù)在453 μm，次峰在8 μm，峰值在400 μm 以上，且位于東部地區(qū)，只能是水成的砂，與河流和湖濱相砂的累積概率曲線(xiàn)峰值相差不大，但符合這種雙峰形式粒度頻率曲線(xiàn)的是近源坡面沖積和洪積泥石流等突發(fā)性的事件，分選較差，且EM3 占比較小，不是主要的物質(zhì)來(lái)源，不能代表時(shí)間段的穩(wěn)定的氣候變化，可能對(duì)東部地區(qū)的洪澇有指示意義，一定程度上代表著較大的降水事件引發(fā)的山洪、滑坡或者泥石流的過(guò)程，指示著突發(fā)性氣候變化對(duì)于沉積的影響。

4.3 過(guò)去2 000 a 氣候重建

圖5 端元各組分與古氣候之間的比較Fig.5 Comparison of the components of terminal element and paleoclimate

通過(guò)粒度與前人所用指標(biāo)的對(duì)比來(lái)探討端元組分所反映的干濕狀況變化，端元法EM1 與偏差法得出的最敏感組分10.473～16.535 μm、董哥洞石筍δ18O，與葛全勝[40]等利用歷史文獻(xiàn)重建的江南地區(qū)干濕狀況和我國(guó)東中部溫度距平做了對(duì)比分析（圖5），在0—300 AD 階段，EM1 增多，EM2 減少，δ18O 含量相對(duì)較高，表明此時(shí)降水較多，沉積顆粒較細(xì)，為較為濕潤(rùn)的時(shí)期，東亞季風(fēng)強(qiáng)盛帶來(lái)的雨水較為充沛，溫度總體處于相對(duì)較暖的階段，但在此階段后期，也就是漢末三國(guó)時(shí)期，氣候開(kāi)始波動(dòng)較大；在300—500 AD 階段，EM1、EM2 都有緩慢下降的過(guò)程，腐殖質(zhì)處于低值，土壤的生產(chǎn)力較低，還原性較強(qiáng)，溫度逐漸轉(zhuǎn)冷，沉積粒徑相對(duì)變粗，整個(gè)江南地區(qū)都較為干旱，此時(shí)屬于由濕潤(rùn)轉(zhuǎn)入較干旱的時(shí)期，降水相對(duì)變少，東亞季風(fēng)減弱，東亞的暖濕氣流減少，江南地區(qū)整體偏干；500—650 AD 時(shí)期，EM1、EM2 再次升高，腐殖化度依然較低，δ18O 相對(duì)升高，細(xì)顆粒再次上升，代表不穩(wěn)定的EM3 減少，江南地區(qū)的大澇次數(shù)也較少，較為穩(wěn)定，而且此時(shí)也有了較為完備的記錄，總體比較可信，溫度也有了上升的趨勢(shì)，這是穩(wěn)定、持續(xù)的濕潤(rùn)時(shí)期，這個(gè)時(shí)間段處于Lamb 定義的黑暗時(shí)期中，但此地緯度較低，氣溫較高，因此與歐洲的結(jié)果有些不同，與百年尺度的冷暖階段存在區(qū)域的差異及北半球溫度變化重建結(jié)果尚存在一定不確定性有關(guān)；650—900 AD 階段，各曲線(xiàn)有良好的一致性，EM1降低，EM2、EM3 升高，細(xì)粒變少，粗粒含量增多，與δ18O 有很好的對(duì)應(yīng)關(guān)系，氣候開(kāi)始由濕潤(rùn)轉(zhuǎn)為干旱，其中出現(xiàn)轉(zhuǎn)折，出現(xiàn)大澇或是大旱次數(shù)也增多，說(shuō)明氣候較不穩(wěn)定，與郝志新[42]等統(tǒng)計(jì)的唐朝后期氣候較為干旱相符。氣候變化也可能是唐朝盛衰變化的轉(zhuǎn)折，大約從900 AD 起，EM2、EM3 變低，腐殖化度升高，沉積顆粒細(xì)粒部分升高，持續(xù)到1100 AD 前后，都是比較濕潤(rùn)的時(shí)期，這階段降水變化比較頻繁，EM1 波動(dòng)較大，水動(dòng)力的變化也較頻繁，但總體是長(zhǎng)期的高值，氣候濕潤(rùn)，季風(fēng)帶來(lái)較多的降水，氣候總體較為溫暖；1100—1200 AD 各個(gè)曲線(xiàn)依舊波動(dòng)比較頻繁，總體的趨勢(shì)是由濕潤(rùn)轉(zhuǎn)向干旱，溫度則是緩慢上升，1200—1400 AD 時(shí)期，EM1處于高值，細(xì)顆粒占比上升，δ18O 含量也上升，反映了江南地區(qū)氣候已經(jīng)由干旱轉(zhuǎn)入濕潤(rùn)，此時(shí)溫度也慢慢轉(zhuǎn)涼，EM2 與腐殖化度的占比都較高，土壤的發(fā)育較好，此時(shí)是比較濕潤(rùn)的涼爽時(shí)期；與鄭景云[26]等得出中世紀(jì)暖期時(shí)，中國(guó)東部的華北地區(qū)相對(duì)干旱，江南則相對(duì)濕潤(rùn)的結(jié)論相似。1400—1800 AD，總體也分為3 個(gè)小階段：EM1 與敏感組分都由高值緩慢轉(zhuǎn)入低值，江南地區(qū)整體分為濕-干-濕的3 個(gè)階段，δ18O 也有相應(yīng)的波動(dòng)，全都反映了小冰期時(shí)期前期較為濕潤(rùn)，中期緩慢變干，后期又慢慢變?yōu)闈駶?rùn)，但溫度總體還是偏低。葛全勝[40]等認(rèn)為1430—1530 年和1630—1740 年是兩個(gè)最為干旱的百年，這在偏差法和EM1 的曲線(xiàn)中都有體現(xiàn)；1800 AD 后受到人類(lèi)活動(dòng)的影響越來(lái)越大。

整體看來(lái)，近2 000 a 來(lái)我國(guó)南方的溫度是總體偏暖，即使是在小冰期，溫度距平在0 以下的也不多見(jiàn)，較暖的空氣增加蒸發(fā)的強(qiáng)度，空氣濕度大的情況下遇冷凝結(jié)，這也是南方地區(qū)降水變率大的一個(gè)重要的因素，而區(qū)域內(nèi)的降水總體偏濕潤(rùn)，整體有略微上升的趨勢(shì)，南方地區(qū)在歷史時(shí)期的水熱變化不均一，有模擬[44]表明這可能與太陽(yáng)輻射和降水帶的移動(dòng)等因素有關(guān)。

5 結(jié)論

（1）玉華山泥炭地沉積物粒度組成以細(xì)粉砂（4～16 μm）為主、中粗粉砂次之（16～64 μm）、黏土（＜4 μm）再次之，砂（＞64 μm）的含量最少，結(jié)合巖性分析表明在過(guò)去2 000 a 沉積環(huán)境表現(xiàn)為湖泊向泥炭演化。

（2）粒度端元分析結(jié)果顯示3 個(gè)端元能夠最大程度地解釋玉華山鉆孔沉積，EM1 代表流水搬運(yùn)而來(lái)的較細(xì)組分，EM2 是反映土壤的成壤改造作用的組分，EM3 代表突發(fā)性的氣候事件帶來(lái)的組分。其中EM1 的變化可以較為準(zhǔn)確地反映該區(qū)域2 000 a以來(lái)的干濕狀況，與各條對(duì)比曲線(xiàn)都有很好的一致性。

（3）玉華山地區(qū)近2 000 a 來(lái)干濕變化的規(guī)律為：0—300 AD 氣候由干轉(zhuǎn)濕，300—500 AD，氣候由濕潤(rùn)轉(zhuǎn)入干旱，整體偏干，500—650 AD 是持續(xù)穩(wěn)定的濕潤(rùn)期，650—900 AD，氣候波動(dòng)較大，變化頻繁，總體偏干；900—1400 AD 是總體較為濕潤(rùn)的時(shí)期，其中有100 a 的波動(dòng)期，1400—1800 AD 為小冰期，氣候變化為濕-干-濕的3 個(gè)階段。

致謝：野外采樣中得到江西師范大學(xué)賈玉連副教授、朱笑虹博士、凌超豪碩士生、凌光久碩士生、南京大學(xué)地理與海洋科學(xué)學(xué)院的本科生李圓圓等人的幫助，在此一并感謝！