臺風(fēng)過程對廣西紅樹林凈生態(tài)系統(tǒng)碳交換的影響研究

孫 明，謝 敏，莫偉華，陳燕麗，潘良浩

（1.廣西科學(xué)院廣西紅樹林研究中心，廣西紅樹林保護(hù)與利用重點(diǎn)實(shí)驗(yàn)室，廣西 北海 536000；2.廣西壯族自治區(qū)氣象科學(xué)研究所，南寧 530022；3.廣西壯族自治區(qū)氣候中心，南寧 530022；）

引言

氣候變化背景下，極端事件增多，特別是加劇了沿海地區(qū)臺風(fēng)、暴雨、風(fēng)暴潮等災(zāi)害［1］，對人類生命財(cái)產(chǎn)、自然生態(tài)系統(tǒng)均產(chǎn)生較為嚴(yán)重的影響。生長于海岸潮間帶的紅樹林生態(tài)系統(tǒng)具有減弱海岸帶侵蝕和臺風(fēng)影響等生態(tài)功能，是海洋向陸地的最后一道“生態(tài)屏障”［2］，既能抵抗風(fēng)浪、保護(hù)海岸，起到防災(zāi)減災(zāi)的作用，還是潛在的重要碳匯，在減緩氣候變化中有重要作用［3-7］。紅樹林生態(tài)系統(tǒng)碳匯功能的演變是全球氣候變暖特征下的研究熱點(diǎn)，在未來氣候變暖的前提下，臺風(fēng)的強(qiáng)度和活動頻率也會加強(qiáng)，更需要對臺風(fēng)干擾下紅樹林凈生態(tài)系統(tǒng)碳交換的變化特征進(jìn)行深入研究。

目前，渦度相關(guān)技術(shù)是全球公認(rèn)的碳循環(huán)標(biāo)準(zhǔn)觀測方法，其可對生態(tài)系統(tǒng)與大氣間物質(zhì)湍流通量進(jìn)行實(shí)時(shí)、非破壞性測量，實(shí)現(xiàn)較大空間尺度上大氣—植被—土壤間的交換過程中碳、水等通量的長期測定，具有直接、準(zhǔn)連續(xù)、無干擾等優(yōu)點(diǎn)。國內(nèi)有關(guān)通量的研究主要集中在森林［8-12］、草地［13-17］、農(nóng)田［18-20］等陸地生態(tài)系統(tǒng)，由于觀測環(huán)境惡劣及影響因素復(fù)雜等原因，國內(nèi)外有關(guān)紅樹林濕地碳通量研究較少，而關(guān)于臺風(fēng)過程對紅樹林碳通量影響的研究更少。本文以廣西紅樹林為研究對象，利用北海生態(tài)氣象觀測試驗(yàn)站的通量和氣象觀測數(shù)據(jù)，研究臺風(fēng)“韋帕”登陸前后紅樹林凈生態(tài)系統(tǒng)碳交換變化特征及其對關(guān)鍵氣象因子的響應(yīng)特征，據(jù)此評估極端天氣影響下紅樹林碳匯功能的變化特征。

1 材料與方法

1.1 研究區(qū)概況

北海紅樹林生態(tài)觀測試驗(yàn)站于2018 年11 月8日建成，位于廣西壯族自治區(qū)北海市國家海洋科技園區(qū)內(nèi)的竹林小海灣（21°27′5.57″N，109°18′4.39″E），是我國東南沿?！拔r塘—海堤—紅樹林”天然次生紅樹林小海灣的典型代表。海堤內(nèi)紅樹林總面積5.5hm2，土壤以礫質(zhì)細(xì)沙土至粗沙土為主，土壤有機(jī)質(zhì)含量較低。由于修建海堤、挖掘?yàn)┠?、毀林修塘等劇烈人為干擾活動，海灣內(nèi)紅樹林斑塊破碎、植株低矮，種類以白骨壤（Aricennia marina）群落，混生少量真紅樹植物秋茄（Kandelia obovata）、桐花樹（Aegiceras corniculata）。紅樹林群落平均高度2.0m，最高達(dá)3.0m，平均基徑8.8cm，最大基徑13.8cm，平均冠幅為1.6m×1.6m，平均蓋度為70%。

北海紅樹林生態(tài)氣象觀測試驗(yàn)站架設(shè)有通量觀測和氣象觀測設(shè)備，能夠?qū)t樹林區(qū)域的三維風(fēng)速、溫度、CO2和H2O 的脈動以及氣象因子進(jìn)行準(zhǔn)確、快速、連續(xù)地觀測。

1.2 影響廣西臺風(fēng)概況

廣西位于華南沿海，地形復(fù)雜造成了降水時(shí)空分布不均，災(zāi)害性臺風(fēng)天氣頻發(fā)，是我國臺風(fēng)多發(fā)地區(qū)之一。據(jù)統(tǒng)計(jì)，2013—2019 年，影響廣西的臺風(fēng)（含熱帶低壓）共計(jì)32 個(gè)，平均每年4.5 個(gè)，2013 年達(dá)到9 個(gè)，為近年來臺風(fēng)影響最頻繁的年份，而2019 年僅有3 個(gè)臺風(fēng)過程，其中“韋帕”是對廣西沿海造成損失最大的一次臺風(fēng)過程。臺風(fēng)“韋帕”概況：2019 年7 月30 日，南海中部海面有熱帶擾動云團(tuán)發(fā)展，17 時(shí)加強(qiáng)為熱帶低壓，隨后于31 日08 時(shí)加強(qiáng)為2019 第7 號臺風(fēng)“韋帕”（熱帶風(fēng)暴級），隨后向西北移動，在發(fā)展為熱帶風(fēng)暴17h 后，于8 月1 日凌晨1 點(diǎn)50 分前后在海南省文昌市沿海登陸，登陸后在文昌市停留了將近7h，于1 日10 時(shí)前后重新移到文昌市東北部近海海面上，之后熱帶風(fēng)暴北上7h后西折，于17 時(shí)40 分前后在廣東省湛江市坡頭區(qū)沿海再次登陸，繼續(xù)西行，于22 時(shí)前后移入北部灣北部海面，之后沿廣西海岸線緩慢移動，于2 日21時(shí)20 分在廣西防城港沿海再次登陸，并于3 日下午14 時(shí)在越南河西省境內(nèi)減弱為熱帶低壓，20 時(shí)中央氣象臺對其停止編號。

1.3 數(shù)據(jù)獲取與處理

1.3.1 氣象數(shù)據(jù)預(yù)處理

氣象觀測要素主要包括空氣溫度（以下簡稱Ta）、土壤表面5cm 溫度（以下簡稱Ts）、降雨量、飽和水汽壓差（Vapor pressure deficit，以下簡稱VPD）以及光合有效輻射（Photosynthetically Active Radiation，以下簡稱PAR），其中，VPD 由開路式紅外氣體分析儀（LI-7500DS）根據(jù)采集的空氣相對濕度和氣溫計(jì)算而來，采集頻率為30min；其他氣象要素由氣象觀測系統(tǒng)采集，采集頻率為10min；為了便于與通量數(shù)據(jù)進(jìn)行相關(guān)性分析以及紅樹林氣象特征分析，我們將氣象數(shù)據(jù)處理為30min 時(shí)間尺度。氣象數(shù)據(jù)分析和圖表繪制主要在Excel 中進(jìn)行，通量與氣象因子之間的響應(yīng)主要利用SPASS 22.0 軟件進(jìn)行擬合和相關(guān)性分析。

1.3.2 通量數(shù)據(jù)預(yù)處理及質(zhì)量控制

選取2019 年7 月29 日—8 月6 日的通量原始觀測數(shù)據(jù)（ghg 格式），采用Li-cor 公司研發(fā)的EddyPro 軟件（版本7.0.6）和TOVI 軟件（版本2.8.1）對通量數(shù)據(jù)進(jìn)行預(yù)處理及質(zhì)量控制。

EddyPro 軟件已成為歐洲綜合碳觀測系統(tǒng)ICOS、美洲通量網(wǎng)AmeriFlux 以及中國生態(tài)系統(tǒng)研究網(wǎng)絡(luò)CERN 處理通量原始數(shù)據(jù)的標(biāo)準(zhǔn)軟件，主要用來將渦度觀測系統(tǒng)采集到的10Hz 原始觀測數(shù)據(jù)處理為30min 間隔的通量數(shù)據(jù)，處理過程主要包括：野點(diǎn)去除、消除延時(shí)、坐標(biāo)旋轉(zhuǎn)、計(jì)算通量原始通量值、超聲溫度修正、頻譜修正、WPL 修正、計(jì)算最終通量值以及Footprint 計(jì)算；數(shù)據(jù)處理采用默認(rèn)的Express Mode 模式，輸出結(jié)果包括各類通量信息和數(shù)據(jù)質(zhì)量控制標(biāo)識等［21］。

由于受到天氣、人為和機(jī)器故障等因素影響［22］，Eddypro 數(shù)據(jù)處理完成后，需要對Eddypro 輸出數(shù)據(jù)進(jìn)行通量數(shù)據(jù)篩選及補(bǔ)插，TOVI 軟件能夠顯著簡化復(fù)雜的通量算法，實(shí)現(xiàn)通量貢獻(xiàn)區(qū)FootPrint 成圖，以及氣象和通量缺失數(shù)據(jù)補(bǔ)插等，本文采用的主要處理過程包括：（1）數(shù)據(jù)完整度檢查：本研究所用數(shù)據(jù)全年數(shù)據(jù)完整度接近100%，其中CO2通量完整度為98.92%；（2）數(shù)據(jù)質(zhì)量控制（Data QA／QC）［23］：剔除因?yàn)闄C(jī)器故障和天氣因素（降雨等）導(dǎo)致的異常數(shù)據(jù)，包含剔除軟件設(shè)定合理閾值范圍之外的數(shù)據(jù)、剔除質(zhì)量控制標(biāo)識為“2”（通量原始數(shù)據(jù)經(jīng)過Eddypro Express 處理后，輸出結(jié)果有一欄為數(shù)據(jù)質(zhì)量評價(jià)，把數(shù)據(jù)質(zhì)量分為0、1、2 三個(gè)等級，0 代表數(shù)據(jù)質(zhì)量優(yōu)異，1 代表數(shù)據(jù)質(zhì)量良好，2 代表數(shù)據(jù)質(zhì)量較差，需剔除）的數(shù)據(jù)［21］、依據(jù)MPT 方法計(jì)算不同季節(jié)摩擦風(fēng)速臨界值并對數(shù)據(jù)剔除，經(jīng)過數(shù)據(jù)質(zhì)量控制后，2019 年全年CO2通量數(shù)據(jù)質(zhì)量較好占比65%，一般情況下數(shù)據(jù)缺失和拒絕率小于50%則具有代表性［24］；（3）通量貢獻(xiàn)源區(qū)分析（Footprint Analysis）：通量觀測的空間代表性指空間某一點(diǎn)的測量值能在多大程度上反映下墊面的平均或累積情況，而對觀測點(diǎn)的通量產(chǎn)生影響的區(qū)域稱為通量貢獻(xiàn)的源區(qū)（source area），可以用足跡（footprint）函數(shù)來定量分析。通量源區(qū)分析是檢查一個(gè)通量站數(shù)據(jù)質(zhì)量的重要方面，以此確定CO2通量足跡的范圍；經(jīng)過質(zhì)量控制后全年的CO2通量足跡范圍均分布于紅樹林群落內(nèi)，無需再進(jìn)行非紅樹林區(qū)域數(shù)據(jù)剔除；（4）數(shù)據(jù)插補(bǔ)（Gap filling）：TOVI 軟件采用MDS 法對缺失部分?jǐn)?shù)據(jù)進(jìn)行插補(bǔ)，得到30min 為間隔、連續(xù)的、質(zhì)量可靠的CO2通量數(shù)據(jù)集。

由渦度觀測系統(tǒng)直接觀測計(jì)算出的CO2通量即為凈生態(tài)系統(tǒng)碳交換量（Net Ecosystem Exchange，以下簡稱NEE），主要是指生態(tài)系統(tǒng)中植物光合作用、冠層空氣中的碳儲存和生物及非生物呼吸消耗的碳排放引起的生態(tài)系統(tǒng)碳蓄積的變化［20］，當(dāng)NEE＜0 時(shí)，生態(tài)系統(tǒng)表現(xiàn)為碳匯；當(dāng)NEE＞0 時(shí)，生態(tài)系統(tǒng)表現(xiàn)為碳源。

2 結(jié)果與分析

2.1 氣象因子變化特征

臺風(fēng)“韋帕”登陸北部灣后對觀測站區(qū)域紅樹林產(chǎn)生影響主要集中在8 月1 日—3 日，為了分析臺風(fēng)登陸前后氣象因子的變化特征，本文選取7 月29日—8 月6 日的氣象觀測數(shù)據(jù)，分為臺風(fēng)登陸前（7月29 日—31 日）、臺風(fēng)登陸（8 月1 日—3 日）及臺風(fēng)登陸后（8 月4 日—6 日）三段進(jìn)行定量的對比分析，選取的氣象要素主要包括：最大風(fēng)速、降雨量、相對濕度、光合有效輻射、氣溫、地表5cm 溫度以及飽和水汽壓差（圖1 所示為各要素變化特征曲線圖）。

強(qiáng)風(fēng)和降雨是臺風(fēng)登陸期間最顯著的特征：如圖1a 所示，臺風(fēng)登陸期間平均最大風(fēng)速為6.28m·s-1，其中，8 月2 日的日均最大風(fēng)速達(dá)到8.07m·s-1，最大風(fēng)速極值出現(xiàn)在8 月2 日12∶00，風(fēng)速達(dá)到11.7m·s-1；與登陸前相比，最大風(fēng)速增大了86.35%，臺風(fēng)登陸后，最大風(fēng)速均值迅速降低為2.52m·s-1；根據(jù)圖4b分析，8 月1 日晚間20∶00 開始雨量逐漸加強(qiáng)，至8月2 日00∶00，小時(shí)降雨量達(dá)到極值50.6mm，為臺風(fēng)期間最大降雨量；之后降水逐漸減弱，其中8 月2 日11∶00、20∶00 以及8 月3 日5∶00 為降水量次峰值點(diǎn)，小時(shí)雨量分別達(dá)到12.6mm、15.2mm 以及21.7mm；根據(jù)日累計(jì)雨量統(tǒng)計(jì)，8 月1 日—3 日的日降水量分別為48.6mm、190.4mm 及76.3mm，其中2 日的總降水量占整個(gè)臺風(fēng)期間降水量的58.67%，而臺風(fēng)登陸前后的總降水量分別為6.9mm 和6.4mm，差異顯著；同時(shí)，對比圖1a 和圖1b，發(fā)現(xiàn)降雨量和最大風(fēng)速出現(xiàn)的時(shí)間及變化趨勢一致。

圖1 臺風(fēng)登陸前后各氣象因子變化特征曲線圖

光合有效輻射是太陽輻射能中能被綠色植物用來進(jìn)行光合作用的那部分能量，與太陽光照直接相關(guān)，一般在中午12 時(shí)達(dá)到最大值；臺風(fēng)登陸前，光合有效輻射的變化范圍在0～2 088μmol·（m2·s）-1浮動，極值呈逐漸上升的趨勢，臺風(fēng)登陸前達(dá)到最大值；臺風(fēng)登陸過程帶來的強(qiáng)風(fēng)和暴雨顯著降低了光合有效輻射，光合有效輻射值在0～644μmol·（m2·s）-1浮動，其中2 日光合有效輻射極值最低，為127μmol·（m2·s）-1；從4 日起隨著臺風(fēng)的消散，光合有效輻射值開始逐漸遞增，但午間極值仍低于登陸前；臺風(fēng)登錄后，平均溫度和最高溫度（包括氣溫和地表5cm 溫度）均會有顯著的下降，如圖1d 所示，7 月29 日—31 日，Ta 和Ts 均值分別為29.7℃和30.6℃，Ts 略高于Ta，但Ts 的峰值溫度變化較Ta 穩(wěn)定，維持在33℃左右，Ta 則在31 日升高至34.1℃，臺風(fēng)登陸前，氣溫顯著升高，而土壤溫度則保持平穩(wěn)；8 月1 日—3 日，Ta 和Ts 均值均降至27.2℃，臺風(fēng)登陸后期，氣溫迅速升高，而Ts 的均值和極值均顯著大于Ta，這可能由于臺風(fēng)的登陸破壞了紅樹林群落結(jié)構(gòu)，冠層的枝干和樹葉大量掉落，太陽輻射較臺風(fēng)登陸前能夠更多地到達(dá)地面，從而使得Ts 快速升溫。

VPD 的變化趨勢與光合有效輻射相似，如圖1（e）所示，臺風(fēng)登陸前的31 日，VPD 峰值達(dá)到最大，達(dá)到2341.88Pa，臺風(fēng)登陸則顯著降低了VPD 的均值和極值（353.8Pa 和711.82Pa），登陸后期VPD 迅速回升，但峰值仍遠(yuǎn)低于登陸前。

2.2 NEE 的變化特征

圖2 為7 月29 日00∶00-8 月6 日23∶00 的紅樹林NEE 的動態(tài)變化曲線，從圖中可以看出NEE 日動態(tài)變化總體呈“U”形分布，日間NEE 小于0，紅樹林表現(xiàn)為碳匯，碳匯峰值一般出現(xiàn)在日間11∶00—14∶00；夜間NEE 大于0，紅樹林表現(xiàn)為碳源；夜間NEE 的碳源峰值在研究期間波動幅度較小，而碳匯峰值則出現(xiàn)較大的波動，臺風(fēng)登陸前，碳匯峰值均值維持在-9.78μmol·（m2·s）-1左右，登陸期間，NEE 碳匯峰值迅速降低，8 月2 日達(dá)到最低值-2.57μmol·（m2·s）-1，隨后8 月3—6 日隨著臺風(fēng)影響的逐漸消散又開始逐漸增大，臺風(fēng)登陸后期NEE 碳匯均值維持在-7.57μmol·（m2·s）-1，和登陸前期相比下降22.60%，紅樹林的碳匯能力受臺風(fēng)過程影響而出現(xiàn)較大幅度的下降。圖3 為紅樹林NEE 的日累積變化特征圖，從圖中可以看出，臺風(fēng)登陸前，NEE 的日累積值逐漸增大，到7.31 日達(dá)到-1.576gC·（m2·d）-1；臺風(fēng)登陸期間，8 月1 日—2 日，NEE 的日累積由負(fù)轉(zhuǎn)正，表明紅樹林生態(tài)系統(tǒng)在此期間由碳匯轉(zhuǎn)變?yōu)樘荚?，向大氣中釋放CO2，其中8 月2 日達(dá)到最大1.076gC·（m2·d）-1，隨后8 月3 日，NEE 日累積值轉(zhuǎn)為負(fù)值，在臺風(fēng)登陸后期，NEE 的日累積又往負(fù)值方向迅速增大，表明紅樹林的碳匯能力逐漸恢復(fù)。

圖2 臺風(fēng)登陸前后NEE 變化特征曲線圖

2.3 臺風(fēng)登陸過程中NEE 對氣象因子的響應(yīng)

為了分析臺風(fēng)期間各氣象因子與NEE 的響應(yīng)關(guān)系，本文利用NEE 與PAR、降雨量、VPD、最大風(fēng)速以及溫度（Ta 和Ts）進(jìn)行相關(guān)分析，由于紅樹林生態(tài)系統(tǒng)在白天主要進(jìn)行光合作用，夜間進(jìn)行呼吸作用，這是兩個(gè)不同的生態(tài)過程，受控的環(huán)境因子也不盡相同，因此本文利用NEE 與PAR、VPD、Ta 及Ts四個(gè)因子按照晝夜分別進(jìn)行小時(shí)尺度的擬合分析；由于NEE 與降雨量及最大風(fēng)速在半小時(shí)尺度上不存在明顯的相關(guān)性，本文對降雨量和最大風(fēng)速進(jìn)行日尺度的擬合分析。

圖3 NEE 日累積值變化特征圖

PAR 是紅樹林生態(tài)系統(tǒng)日間NEE 的重要影響因子，因此僅利用日間NEE 與PAR 進(jìn)行擬合，由圖4可以看出NEE 與PAR 呈現(xiàn)出較好的對數(shù)函數(shù)關(guān)系，臺風(fēng)登陸期間，日間PAR 大幅降低，因此導(dǎo)致NEE 的碳匯峰值也隨之減小。

小時(shí)尺度日間：VPD 與NEE 的相關(guān)系數(shù)最大，NEE 隨著VPD 的增大而往負(fù)值方向變化，即紅樹林碳匯能力逐漸增強(qiáng)，表明VPD 對NEE 具有一定的促進(jìn)作用；VPD 越高，NEE 碳匯峰值也往負(fù)值方向移動，碳匯能力增強(qiáng)；臺風(fēng)登陸期間，VPD 由31 日的2 341.88Pa 迅速降低到1 日的165.4Pa，陡然降低的VPD 導(dǎo)致NEE 迅速往正值方向變化，及VPD的降低顯著降低了紅樹林碳匯能力；日間Ta 對NEE 的直接作用稍強(qiáng)于Ts，臺風(fēng)登陸期間，溫度增加會增強(qiáng)會使NEE 往負(fù)值移動從而碳匯能力得到增強(qiáng)，臺風(fēng)期間的降溫也在一定程度上對NEE 起抑制作用。

小時(shí)尺度夜間（圖5）：NEE 與各因子的相關(guān)性均大幅下降，表明紅樹林夜間的呼吸作用受氣象因子的影響較小。

日尺度上（圖6）：日尺度上最大風(fēng)速和降雨量對NEE 的影響顯著強(qiáng)于半小時(shí)尺度，R2均達(dá)到0.6以上，最大風(fēng)速與NEE 日累積量的相關(guān)性略強(qiáng)于降雨量，R2達(dá)到0.64，并且二者對NEE 日累積的影響均為負(fù)向，即最大風(fēng)速越強(qiáng)，降雨量越大，NEE 的日累積量會往正值方向變化，即紅樹林的日碳匯由負(fù)值轉(zhuǎn)變?yōu)檎?，紅樹林生態(tài)系統(tǒng)由原來的吸收CO2轉(zhuǎn)變?yōu)獒尫臗O2，由此我們可以看出，臺風(fēng)登陸過程中的強(qiáng)風(fēng)和強(qiáng)降水會對紅樹林NEE 的日累積起抑制作用，為NEE 日累積的限制因子。

圖4 日間NEE 與各氣象因子響應(yīng)圖

圖5 夜間NEE 與各氣象因子響應(yīng)圖

圖6 NEE 與降雨量及最大風(fēng)速響應(yīng)圖

3 結(jié)論與討論

由于廣西沿海臺風(fēng)頻發(fā)，紅樹林經(jīng)常遭受臺風(fēng)影響，本文主要以“韋帕”臺風(fēng)為例，分析了臺風(fēng)登陸過程中廣西紅樹林NEE 的變化特征及其對主要?dú)庀笠蜃拥捻憫?yīng)關(guān)系，研究結(jié)論主要有以下幾點(diǎn)：

（1）臺風(fēng)登陸前后氣象因子的變化特征：登陸前，Ta、Ts、PAR 及VPD 均呈逐日遞增的趨勢，在31日達(dá)到極大值；登陸時(shí)，Ta、Ts、PAR 及VPD 均大幅度下降，同時(shí)，降雨量和最大風(fēng)速迅速增大；登陸后期，隨著降雨和強(qiáng)風(fēng)的消散，Ta、Ts、PAR 及VPD 四個(gè)因子開始逐漸升高，除Ts 外，其他三個(gè)氣象因子臺風(fēng)后的極值均低于臺風(fēng)前，這可能由于臺風(fēng)登陸時(shí)強(qiáng)風(fēng)對紅樹林冠層的強(qiáng)烈擾動，造成大量枝干損傷，植被覆蓋度下降，從而導(dǎo)致臺風(fēng)后期太陽輻射能夠更多地到達(dá)地表，造成土壤溫度比于臺風(fēng)登陸前更高的現(xiàn)象。

（2）紅樹林NEE 的變化特征：紅樹林NEE 的日動態(tài)變化特征變化曲線與PAR 等氣象因子變化特征趨勢一致，29—31 日碳匯峰值呈逐漸遞增并在31日達(dá)到碳匯峰值，隨著臺風(fēng)的登陸碳匯峰值迅速降低，臺風(fēng)后期開始逐漸遞增，但碳匯峰值與臺風(fēng)登陸前相比顯著降低；在紅樹林NEE 的日累積量上，29—31 日，NEE 的日累積量逐日遞增，而在8 月1—2 日，在強(qiáng)風(fēng)和降雨的作用下，紅樹林NEE 日累積由負(fù)值轉(zhuǎn)為正值，碳匯能力減弱，呼吸作用增強(qiáng)；這與Barr 等［25-26］研究結(jié)論一致，臺風(fēng)的干擾會使生態(tài)系統(tǒng)的生產(chǎn)力降低及呼吸作用增加。8 月3 日往后，隨著臺風(fēng)的消散，NEE 日累積量開始開始由正轉(zhuǎn)負(fù)并逐漸遞增，紅樹林碳匯能力逐漸恢復(fù)，表明紅樹林抗臺風(fēng)擾動的能力較強(qiáng)。

（3）紅樹林NEE 與各氣象因子的響應(yīng)關(guān)系：小時(shí)尺度上，VPD 是紅樹林生態(tài)系統(tǒng)日間NEE 的主要影響因子，其次為PAR；在日尺度上，最大風(fēng)速和降雨量對NEE 的影響顯著，二者均為NEE 日累積的限制因子，能夠顯著降低NEE 的日累積量。