曹玉平 袁熱林 焦樹林 張倩 鄧飛艷

摘要：為探究筑壩對河流溶解無機碳、溶解無機碳氮空間分布的影響，2016年9月（豐水期）在巖溶區(qū)深水水庫（龍灘水庫）進行分層采樣，獲取水體主要的理化參數(shù)，運用經(jīng)典統(tǒng)計學(xué)和地統(tǒng)計學(xué)方法揭示了水中溶解無機碳（DIC）、氮（D/Ⅳ）及溶解無機碳氮比（C/N）的空間變異特征。結(jié)果表明：龍灘水庫夏季出現(xiàn)明顯熱分層現(xiàn)象，導(dǎo)致總?cè)芙夤腆w、酸堿度、堿度、DIC、DIN出現(xiàn)空間分層;除水溫、pH外，水中主要理化參數(shù)在空間上變異較大，DIC、DIN、C/N的空間變異總體表現(xiàn)為DIC>C/N>DIN;庫區(qū)DIC、DIN、C/N的水平變異呈現(xiàn)從壩前向庫尾逐漸遞減的趨勢;由試驗隨機誤差產(chǎn)生的空間異質(zhì)性較小，自然結(jié)構(gòu)因素對DIC、DIN、C/N引起的空間異質(zhì)性起主導(dǎo)作用。

關(guān)鍵詞：溶解無機碳;溶解無機氮;半方差變異函數(shù);空間異質(zhì)性;熱分層;巖溶區(qū)

中圖分類號：X524，P95

文獻標志碼：A

doi：10.3969/j .issn.1000- 13 79.2019.02.019

近年來全球變暖受到普遍重視[1-2]，全球碳循環(huán)問題成為研究熱點，河流碳循環(huán)作為全球碳循環(huán)的重要組成部分受到廣泛關(guān)注[3-5]。水電作為清潔能源在我國不斷推廣和發(fā)展，截至2015年年底我國共有水庫97 988座，其中西南地區(qū)19 774座，占全國的20.18%[6]。水庫修建給河流水化學(xué)環(huán)境、水動力條件帶來了嚴重影響，由筑壩產(chǎn)生的“水庫效應(yīng)”成為當前研究的熱點之一。筑壩阻斷了河流天然連續(xù)性，改變了水體自然性質(zhì)，使之成為介于湖泊和河流之間的半自然態(tài)水體[7-8]，改變了自然狀態(tài)下河流生態(tài)系統(tǒng)的碳、氮循環(huán)，形成新的水域生態(tài)系統(tǒng)和碳、氮循環(huán)模式[9-10]。目前關(guān)于溶解無機碳、溶解無機氮空間分布特征的研究主要集中在一些河流河口或大中型水庫[11-13]。焦樹林等[12]認為西江河口段溶解無機碳穩(wěn)定同位素組成平均值夏秋季小于冬春季，存在明顯的季節(jié)性差異：吳起鑫等[13]認為筑壩沒有改變河流溶解無機碳的地球化學(xué)行為過程：趙敏等[14]認為，土地利用方式改變影響水中溶解無機碳、氮的組成和變化，在巖溶地區(qū)水中溶解無機碳受石漠化影響比受原始森林影響的比重大，土壤巖溶作用對水中溶解無機碳起控制作用。但關(guān)于水庫中溶解無機碳、氮空間分異格局的相關(guān)研究較少。

龍灘水庫是一座兼具防洪、發(fā)電、航運、養(yǎng)殖、灌溉等多功能的大型水利樞紐工程，在國內(nèi)其裝機容量和庫容僅次于長江三峽水庫。龍灘水庫水溫呈明顯季節(jié)性分層，目前對龍灘水庫的研究主要集中在水體富營養(yǎng)化和浮游植物群落空間分布兩方面[15-16]。本文運用半方差變異函數(shù)理論模型探討了巖溶區(qū)龍灘水庫夏季溶解無機碳、氮的空間分異格局。

1 材料與方法

1.1 試驗采樣與分析

龍灘水庫位于紅水河上游廣西天峨縣境內(nèi)，壩址距離天峨縣城15 km.壩址以上流域面積98 500 km，占紅水河流域面積的71%，其庫容273億m，回水長度120 km，正常蓄水位378 m，最大水深173.5 m，平均水深92.7 m，裝機容量為540萬kW。龍灘水庫回水區(qū)域狹長，為保證采樣點均勻分布，從壩前逆流向上游設(shè)置7個采樣點，并用GPS確定采樣點位置（見圖1）。2016年9月（豐水期）在巖溶區(qū)深水水庫（龍灘水庫）進行分層采樣。采用自制雙通采水器采集垂直斷面水樣，從表層至底層間隔10 m分層采樣，直到采取庫底沉積物為止?，F(xiàn)場利用美國麥隆公司生產(chǎn)的Ultru -Ⅱ（6P）水質(zhì)參數(shù)儀測定水樣的溫度（T）、氧化還原電位（ORP）、酸堿度（pH值）、電導(dǎo)率（Cond）、總?cè)芙夤腆w（TDS）含量，并用稀鹽酸滴定水樣堿度（ALK）。每個采樣點用水樣清洗3次后裝瓶，密封帶回實驗室用HI83200色譜儀測定樣品營養(yǎng)鹽含量。

1.2 半方差變異函數(shù)模型

在水文學(xué)、地質(zhì)學(xué)、生態(tài)學(xué)中許多變量都具有空間分布的特點，當一個變量呈現(xiàn)為空間分布時，就稱之為區(qū)域化變量[17]。水中溶解無機碳、氮可視為區(qū)域化變量，具有結(jié)構(gòu)性（空間自相關(guān)性）和隨機性特征，同時庫區(qū)中溶解無機碳、氮變量滿足二階平穩(wěn)假設(shè)條件，樣本空間足夠大時其半方差變異函數(shù)γ（h）的計算式為[18]Kolmogorov-Smironov（ K-S）在SPSS19.0中完成，用地統(tǒng)計學(xué)軟件包GS+ for Windows 9.1進行半方差函數(shù)γ（h）的計算，運用Surferll.0軟件進行克里金插值，用Origin9.0軟件進行基本統(tǒng)計分析。

2 結(jié)果與討論

2.1 統(tǒng)計分析與正態(tài)分布檢驗

水庫水體主要理化性質(zhì)描述性統(tǒng)計結(jié)果見表1（其中“*”表示通過顯著性水平為0，05的顯著性檢驗;變異系數(shù)、偏度、峰度、檢驗概率P為無量綱）。巖溶區(qū)水體受巖溶作用的影響，pH值為7.63～ 8.30，平均值為8.02，水體總體上表現(xiàn)為弱堿性。夏季水體容易出現(xiàn)熱分層現(xiàn)象[15]，龍灘水庫出現(xiàn)明顯的水溫熱分層現(xiàn)象。0—10 m水體水溫形成溫躍層，10 m以下水體水溫相對較穩(wěn)定，見圖2（a）。表層水溫變化梯度大于庫底層的，兩者溫度相差5.8℃。壩前水體的水溫相對變化小，分析其原因，壩前水體的涵熱作用使水溫的變幅較小[21]。由于夏季存在熱分層效應(yīng)，因此水庫水體溫度、水生生物的分布及新陳代謝的強度和方向在不同水層存在顯著差異，水溫形成的密度分層限制上下水體交換，致使水庫水體總?cè)芙夤腆w、pH、堿度等物理化學(xué)性質(zhì)差異較大且出現(xiàn)明顯的分層。從空間分布看，總?cè)芙夤腆w濃度表層小于底層，隨著水深的增加濃度隨之升高，在距離大壩40 km的水體中間形成一個高濃度區(qū)，網(wǎng)箱養(yǎng)魚眾多是該區(qū)域總?cè)芙夤腆w濃度大的主要原因，見圖2（b）。pH值在整個水體空間分布相對均勻，在受支流匯人擾動的條件下，pH值變化幅度較大，見圖2（c）。堿度也出現(xiàn)分層結(jié)構(gòu)，壩前水體的堿度變化大于庫區(qū)上游地區(qū)的。在距大壩20～30 km處，堿度形成高值區(qū)域，分析認為主要受濛江和漕渡河支流匯人的影響，夏季上游巖溶溶蝕作用帶來豐富的碳酸鹽類溶解物，造成該區(qū)域堿度濃度升高，見圖2（d）。從各理化參數(shù)變異系數(shù)分析可知，除水溫、pH外，其余理化參數(shù)的空間變異較大，總體表現(xiàn)為NH4-N>ORP>TDS=ALK>DIC>C/N>DIN.

半方差函數(shù)主要通過區(qū)域化變量分割等距離采樣點間的差異來研究變量的空間結(jié)構(gòu)和空間相關(guān)性。但空間相關(guān)分析的變量必需滿足正態(tài)分布，并且由隨機抽樣的方式獲得[22]。為了檢驗試驗數(shù)據(jù)是否符合半方差函數(shù)的分析條件，采用SPSS中Kolmogorov - Smi-ronov（K-S）正態(tài)分布檢驗概率（P）進行檢驗。檢驗時顯著性水平設(shè)定為α= 0.05，若P >0.05.則認為數(shù)據(jù)服從正態(tài)分布。由表1可知，水中DIC、DIN、C/N均大于0.05，說明數(shù)據(jù)服從正態(tài)分布，可進行空間自相關(guān)分析，并且DIC、C/N的P明顯優(yōu)于DIN的。

2.2 庫區(qū)DIC、DIN、C/N各向同性下空間結(jié)構(gòu)性

如果區(qū)域化變量在采樣尺度上具有空間依賴性（或者空間相關(guān)性），那么半方差函數(shù)值會隨著滯后距的增加而增大，并且在超過變程后逐漸趨近于基臺值或圍繞基臺值波動[22]。從圖3和表2可知，水體中DIC、DIN、C/N具有明顯的基臺值，說明水中溶解性DIC、DIN、C/N的分布具有一定的空間依賴性和結(jié)構(gòu)性;基臺值DIN>C/N>DIC。本研究中水中溶解性C、N的分析采用不同的擬合模型，DIC適合球狀模型、DIN適合高斯模型、C/N適合系數(shù)模型。確定系數(shù)（R2）和殘差（RSS）可以評價模型擬合對研究對象的解釋效率。從表2可知，DIC、C/N的半方差函數(shù)γ（h）與空間距離h的理論模型解釋率比DIⅣ的高，DIC、C/N、DIN的解釋率分別為48%、56%、28%。

通常認為區(qū)域化變量的空間異質(zhì)性（SHZ）由空間自相關(guān)（SHA）和隨機誤差（SHR）兩部分組成。B表示空間自相關(guān)引起的空間異質(zhì)性，塊金值Bo表示隨機誤差的空間異質(zhì)性，所以基臺值Bo+B越大，區(qū)域化變量的空間異質(zhì)性越高。B/（B0+B）反映了結(jié)構(gòu)因素SHA對SHZ的貢獻程度，Bo/（B0+B）反映了隨機誤差SHR引起的空間異質(zhì)性占SHZ的比例。由表2可知，水中DIC、DIN、C/N的B/（B0+B）值均較高，說明結(jié)構(gòu)因素對其空間異質(zhì)性的貢獻率分別高達99.97%、99.91%、76.56%。按照區(qū)域化變量空間相關(guān)性的分級標準[23-24]，當B0/（B0+B）≥0.75時，變量的空間相關(guān)性很弱;當0.25≤Bo/（Bo+B）<0.75時，變量具有中等程度的空間相關(guān)性;當B0/（ B0+B）<0.25時，變量具有強烈的空間相關(guān)性。結(jié)合該分級標準，庫區(qū)水中溶解性DIC、DIN、C/N的B0/（B0 +B）值都小于<0.25%，說明變量具有強烈的空間相關(guān)性。隨機試驗產(chǎn)生的誤差（B0）對水中溶解性C、N的影響較小，受水體本身的性質(zhì)和流域物質(zhì)輸出的影響較大，研究表明水體中溶解性離子濃度是區(qū)域氣候、地層、巖石、土壤、農(nóng)業(yè)活動和人為排放等環(huán)境因子的函數(shù)[10]。同時，水生生物的光合和呼吸作用是引起水體C、N空間分布差異的主要原因之一。

變程反映的是區(qū)域化變量的空間影響范圍。從表2和圖3可知，水中C/N值的變程為0.15 km，是最小的，說明水中C、N比值在較短的距離內(nèi)存在空間異質(zhì)性。相比而言，水中溶解無機碳和溶解無機氮的變程分別為2.00、6.47 km，說明水中溶解無機碳在空間距離超過2 km后空間相關(guān)性消失，而水中溶解無機氮在較長的距離內(nèi)存在空間結(jié)構(gòu)異質(zhì)性，當超過變程6.47km后空間相關(guān)性消失。

2.3 庫區(qū)DIC、DIN、C/N各向異性下空間結(jié)構(gòu)性

為了研究水體中DIC、DIN和C/N在不同方向上的特點，即各向異性，對不同方向的半方差函數(shù)進行計算。計算時將全方位平均分為4個方向，即0°、45°、90°、135°，分別代表東一西、東北一西南、南一北、西北一東南方向。從表3可知，水體中DIC、DIN和C/N的基臺值遠大于塊金值，表明在不同方向上三者同樣具有明顯的空間結(jié)構(gòu)。其中DIC的半方差函數(shù)γ（h）與空間距離h的關(guān)系符合指數(shù)模型.DIN與C/N的符合高斯模型。從確定系數(shù)R2和殘差RSS對模型的解釋率可知，與各向同性半方差函數(shù)模型相比，DIC、DIN和C/N各向異性半方差函數(shù)模型的解釋率都較高，解釋率分別為71%、62%、54%。另外，從表3可以看出，水體中DIC、DIⅣ和C/N的B/（B0+B）值相差不大，分別為0.778、0.776、0.736，說明DIC、DIⅣ和C/N在各方向上的空間變異程度差異明顯。各向異性下，水中DIC、DIN和C/N的B0/（B0 +B）值分別為0.222、0.224、0.264，相較于各向同性的空間相關(guān)性來說，水中C、N各向異性的空間相關(guān)性降低，C/N的空間相關(guān)性由強烈的空間相關(guān)性變?yōu)橹械葟姸鹊目臻g相關(guān)性。綜上所述，結(jié)構(gòu)因素是庫區(qū)DIC、DIN、C/N的空間異質(zhì)性的主要影響因素，與各向同性下的分析結(jié)果一致。

DIC、DIN和C/N各向異性的半方差變異函數(shù)與空間距離關(guān)系見圖4，結(jié)合表3，水體中DIC、DIN和C/Ⅳ三者在空間上各方向的變程各不相同。DIC、DIN在空間各方向上的變異以0°方向變異為主，也就是說從壩前到庫區(qū)上游方向DIC、DIN空間變化最大;C/N在0°、45°兩個方向上變化最大。

2.4 庫區(qū)DIC、DIN、C/N及NH4 -N的空間分布格局

水中溶解性C、N的組成和變化是一個十分復(fù)雜的綜合過程，受流域侵蝕、地形、氣候、植被、人為干擾活動，以及水中生物光合、呼吸作用及硝化細菌硝化作用等影響[11]。水庫蓄水后改變河流原始狀態(tài)下C、N的地球化學(xué)過程，水庫對水體營養(yǎng)物質(zhì)具有明顯截留效應(yīng)[25]。水中溶解性C、N具有空間變異性?；诎敕讲詈瘮?shù)的建立，通過正態(tài)分布檢驗，水中DIC、DIN、C/N、NH4 -N均服從正態(tài)分布，采用Surferll.0軟件中克里金插值法繪制水中DIC、DIN、C/N、NH4 -N空間等值線分布圖，以直觀反映水中C、N在空間上的分布特征。龍灘水庫水溫在夏季出現(xiàn)明顯熱分層效應(yīng)，水庫的熱分層阻礙了水中其他物質(zhì)在上下層之間的交換。從圖5可以看出，DIC、DIN、C/N、NH4 -N呈明顯分層。

在垂直方向上.DIC表層水體濃度變化梯度大于底層水體的，隨著水深的增加DIC濃度升高，夏季水生生物光合作用強，受光透深度、水體熱分層等因素影響，上層水體是水庫初級生產(chǎn)力的主要產(chǎn)生區(qū)域。上層水體光合作用一定程度上降低了水體的DIC濃度[26]，庫底DIC長期滯留，補充表層DIC，是表層DIC的重要來源。受水庫蓄熱作用影響[21]，在壩前深度為10、30、40 m處出現(xiàn)了半閉合DIC等值線，見圖5（a）。

DIN空間分布大致上是中間水層濃度大，上下兩層濃度小，呈斑塊狀分布，空間上相關(guān)性不強，具有高度空間異質(zhì)性，見圖5（b）。這與張壘[15]等研究的龍灘水庫總氮分布特征一致，中表層水體氮含量最大。龍灘水庫屬于深水水庫，庫內(nèi)水體更新交換的時間較長，從圖5（b）可以得出這樣的結(jié)論，水庫底層長期滯留N元素，形成DIN富集區(qū)。

水體中C/N空間分布受N03 -N、NH4 -N、DIC、DIⅣ等的綜合影響。研究區(qū)C/N的空間分布與DIC的具有相似性，在縱向上出現(xiàn)高值斑塊區(qū)，表層水體的C/N出現(xiàn)明顯分層，見圖5（c），主要是受庫區(qū)一支流交互作用的影響。夏季降水充沛，流域侵蝕作用強烈，上游來水攜帶了大量的營養(yǎng)鹽，加上受水體頂托作用的影響，在干、支流交匯處容易出現(xiàn)營養(yǎng)鹽富集。在深水水庫中，水體達到一定深度時，水中光合作用減弱，水生生物以呼吸作用為主，造成水庫底層缺氧的環(huán)境。在缺氧的條件下，反硝化細菌反硝化作用明顯，還原硝酸鹽，釋放出分子態(tài)氮（N2）或一氧化二氮（N20），造成水中溶解性N濃度降低。NH4 -N在水中分布比較紊亂，高、低值斑塊交錯出現(xiàn)，上游地區(qū)NH4 -N分層現(xiàn)象比壩前水體明顯（圖5（d））。在整個水體中NH4 -N/DIN為1%，水中DIN的形態(tài)主要是N02 -N，比例達98%。

3 結(jié)語

筑壩改變原自然狀態(tài)河流水環(huán)境、水動力條件，導(dǎo)致水中溶解無機碳、氮的生物地球化學(xué)行為過程發(fā)生改變。筆者運用經(jīng)典統(tǒng)計學(xué)，結(jié)合地統(tǒng)計學(xué)中半方差變異函數(shù)理論模型，探討了巖溶區(qū)深水水庫溶解無機碳、氮的空間分異特征。龍灘水庫夏季水溫出現(xiàn)明顯熱分層效應(yīng)，導(dǎo)致水中TDS、pH值、ALK、DIC、DIN、C/Ⅳ及NH4 -N出現(xiàn)空間分層現(xiàn)象，除水溫、pH值外，水中主要理化參數(shù)在空間上變異性相差較大，總體表現(xiàn)為NH4 -N>ORP> TDS= ALK> DIC> C/N> DIⅣ;庫區(qū)DIC、DIN、C/N的水平變異具有明顯的水平空間結(jié)構(gòu)，水平變異呈現(xiàn)從壩前向庫尾逐漸遞減的趨勢;DIC、DIN、C/N在空間上具有強烈的空間自相關(guān)性，由試驗隨機誤差（Bo）產(chǎn)生的空間異質(zhì)性較小，自然結(jié)構(gòu)因素（B）對DIC、DIN、C/N引起的空間異質(zhì)性起主導(dǎo)作用;DIN在空間分布上出現(xiàn)明顯的高、低值斑塊區(qū)，表現(xiàn)為高異質(zhì)性，總體上呈現(xiàn)中間高、兩端低的分布格局。

參考文獻：

[1]常跟應(yīng)，黃夫朋，李曼，中國公眾對全球氣候變化認知與支持減緩氣候變化政策研究：基于全球調(diào)查數(shù)據(jù)和與美國比較視角[J].地理科學(xué)，2012，32（ 12）：1481-1487.

[2]于洋，陳才，中國二氧化碳凈排放的時空特性分析[J].地理科學(xué)，2013，33（ 10）：1173-1179.

[3] 高全洲，沈承德，河流碳通量與陸地侵蝕研究[J].地球科學(xué)進展，1998，13（4）：369-375.

[4] 魏秀國，李寧利，沈承德，等，西江水體有機碳含量變化及懸浮物碳同位素的意義[J].地理科學(xué)，2011，31（2）：166-171.

[5] 張永領(lǐng)，楊小林，張東，小浪底水庫影響下的黃河花園口站和小浪底站pC02特征及擴散通量[J].環(huán)境科學(xué)，2015，36（1）：40-48.

[6] 中華人民共和國國家統(tǒng)計局.2016中國統(tǒng)計年鑒[M].北京：中國統(tǒng)計出版社，2016：1- 56.

[7]HOOPER R P，AULENBACH B T，KELLY V J.The Na-tional Streamquality Accounting Network：a Flux-Based Ap-proach to Monitoring the Water Quality of Large Rivers[J].Hydrological Processes， 2001， 15（7）： 1089-1106.

[8]NAIMAN R J，MELILLO J M， LOCK M A，et al.Longitu-dinal Pattems of Ecosystem Processes and Community Struc-ture in a Subarctic River Continuum[J].Ecology， 1987， 68（5）：1139-1156.

[9]李干蓉，劉叢強，陳椽，等，貓?zhí)恿饔蛱菁壦畮煜囊磺锛竟?jié)溶解無機碳（ DIC）含量及其同位素組成的分布特征[J].環(huán)境科學(xué)，2009，30（10）：2891-2897.

[10] 劉叢強，汪福順，王雨春，等，河流筑壩攔截的水環(huán)境響應(yīng)：來自地球化學(xué)的視角[J].長江流域資源與環(huán)境，2009，18（4）：384-396.

[11] 范志偉，郝慶菊，黃哲，等，三峽庫區(qū)水體中可溶性C、N變化及影響因素[J].環(huán)境科學(xué)，2017，38（1）：129-137.

[12] 焦樹林，陶貞，高全洲，等，西江河口段溶解無機碳穩(wěn)定同位素組成的時空變化[J].地理學(xué)報，2008，63（5）：526-533.

[13] 吳起鑫，韓貴琳，唐楊，三峽水庫壩前水體水化學(xué)及溶解無機碳時空分布特征[J].環(huán)境科學(xué)學(xué)報，2012，32（3）：654-661.

[14] 趙敏，曾成，劉再華，土地利用變化對巖溶地下水溶解無機碳及其穩(wěn)定同位素組成的影響[J].地球化學(xué)，2009，38（6）：565-572.

[15] 張壘，李秋華，黃國佳，等，亞熱帶深水水庫——龍灘水庫季節(jié)性分層與富營養(yǎng)化特征分析[J].環(huán)境科學(xué)，2015，36（2）：438-447.

[16] 王崇，憨雪瑩，黃道明，等，紅水河龍灘水庫浮游植物的群落結(jié)構(gòu)特征[J].水生態(tài)學(xué)雜志，2013，34（6）：34-40.

[17] MATHERON G.Principles of Geostatistics[J]. EconomicGeology， 1963， 58（8）：1246-1266.

[18]徐建華，計量地理學(xué)[M].2版，北京：高等教育出版社，2014：171-177.

[19] ALLAN J D， CASⅡLLO M M.Stream Ecology： Structureand Function of Running Waters [M]. 2nd Edition.Berlin： Springer Press， 2007： 256-264.

[20]李哲，郭勁松，方芳，等，三峽水庫小江回水區(qū)不同TN/TP水平下氮素形態(tài)分布和循環(huán)特點[J].湖泊科學(xué)，2009，21（4）：509-517.

[21] 李冰凍，李嘉，李克鋒，等，二灘水庫壩前及下泄水體水溫分布現(xiàn)場觀測與分析[J].水利水電科技進展，2009，29（4）：21-23.

[22]王政權(quán)，地統(tǒng)計學(xué)及其在生態(tài)學(xué)中的應(yīng)用[M].北京：科學(xué)出版社，1999：35-43.

[23] 劉璐，曾馥平，宋同清，等，喀斯特木論自然保護區(qū)土壤養(yǎng)分的空間變異特征[J].應(yīng)用生態(tài)學(xué)報，2010，21（7）：1667-1673.

[24] 牟曉杰，孫志高，劉興土，黃河口濱岸潮灘濕地土壤碳、氮的空間分異特征[J].地理科學(xué)，2012，32（12）：1521-1529.

[25] 張恩仁，張經(jīng)，三峽水庫對長江N、P營養(yǎng)鹽截留效應(yīng)的模型分析[J].湖泊科學(xué)，2003，15（1）：41-48.

[26] 喻元秀，劉叢強，汪福順，等，烏江流域梯級水庫中溶解無機碳及其同位素分異特征[J].科學(xué)通報，2008，53 （16）：1935 - 1941.