小型水體（< 1 km²）占全球內(nèi)陸水體總數(shù)量的99%以上，其蒸發(fā)估算是預(yù)測(cè)未來(lái)水資源儲(chǔ)量的關(guān)鍵環(huán)節(jié)（Verpoorter et al., 2014; Messager et al., 2016)。在小型水體上，傳統(tǒng)蒸發(fā)方法的應(yīng)用準(zhǔn)確度大大降低，例如Priestly-Taylor模型的模擬結(jié)果因平流效應(yīng)的增強(qiáng)而存在較大偏差（Bello and Smith, 1990），通量梯度法和渦度相關(guān)法因小型水體的風(fēng)浪區(qū)有限和儀器安裝高度的限制導(dǎo)致數(shù)據(jù)有效性較低（Zhao et al., 2019），然而穩(wěn)定同位素質(zhì)量守恒法在小型水體上的應(yīng)用凸顯了它獨(dú)特的優(yōu)勢(shì)，為蒸發(fā)估算提供了新的思路和方法，但是該方法的技術(shù)難題在于其耦合的Craig-Gordon模型中同位素動(dòng)力分餾系數(shù)ε_k的準(zhǔn)確量化。

同位素動(dòng)力分餾效應(yīng)是指在水循環(huán)的相變過(guò)程中（例如蒸發(fā)、降水等），因互為同位素的不同原子（例如¹⁶O和¹⁸O）的分子擴(kuò)散速率不同造成水體中各同位素富集程度發(fā)生變化的一種過(guò)程。同位素動(dòng)力分餾系數(shù)ε_k表征該效應(yīng)的強(qiáng)弱，代表大氣湍流削弱動(dòng)力分餾效應(yīng)的強(qiáng)弱，目前如何量化這種削弱作用仍不明確，致使動(dòng)力分餾效應(yīng)的估算不準(zhǔn)確（Craig and Gordon, 1965；Horita et al., 2008；Gonfiantini et al., 2018）。傳統(tǒng)ε_k算法方法主要有兩種：一種為基于實(shí)驗(yàn)室蒸發(fā)皿觀測(cè)總結(jié)的經(jīng)驗(yàn)值，ε_k分別為12.5‰（²H）和14.2‰（¹⁸O），簡(jiǎn)稱(chēng)LK算法；另一種基于海洋觀測(cè)數(shù)據(jù)和物理模型推導(dǎo)，通過(guò)風(fēng)速進(jìn)行求算 (Merlivat and Jouzel, 1979) ，簡(jiǎn)稱(chēng)OS算法。有研究對(duì)比發(fā)現(xiàn)基于OS算法（ε_k¹⁸O = 6.2‰）計(jì)算的太湖年蒸發(fā)量比基于LK算法估算的蒸發(fā)量高72%（Xiao et al., 2017）。因此，ε_k的兩種算法在實(shí)際應(yīng)用中產(chǎn)生的分歧極大，自然水體的ε_k真實(shí)值需要野外實(shí)驗(yàn)進(jìn)一步觀測(cè)研究。迄今為止研究ε_k的野外實(shí)驗(yàn)僅有Gonfiantini等人(2018)和Xiao等人（2017）。前者主要通過(guò)Unified Craig-Gordon模型基于前人觀測(cè)的水體同位素組分?jǐn)?shù)據(jù)推導(dǎo)ε_k，缺點(diǎn)是模型中大氣水汽同位素組分δ_V輸入的是一個(gè)固定值；后者是基于通量梯度法在太湖（2400 km²）觀測(cè)δ_E計(jì)算ε_k，但是太湖屬于超大型水體（>1000 km²），風(fēng)浪區(qū)大，風(fēng)速或者湍流強(qiáng)度也大，估算的ε_k較小。所以ε_k的真實(shí)值需要在分布最為廣泛的小型水體上進(jìn)一步得到檢驗(yàn)。為探討小型水體蒸發(fā)中的最佳動(dòng)力分餾，本研究在安徽的官渡實(shí)驗(yàn)基地，選擇了一個(gè)小型水體（魚(yú)塘，面積約為6900 m²）和常規(guī)氣象觀測(cè)的大小兩種蒸發(fā)皿（直徑20 cm和60 cm）進(jìn)行觀測(cè)實(shí)驗(yàn)，反算各水體的ε_k（圖1）。結(jié)果表明：（1）基于IMB模型和CG模型估算的魚(yú)塘、大蒸發(fā)皿和小蒸發(fā)皿ε_k最優(yōu)值分別為10.2 ± 4.9 ‰、14.3 ± 6.8 ‰和7.0 ± 3.1 ‰（圖1）。（2）在小型水體中，動(dòng)力分餾系數(shù)的參數(shù)化方案中，OS算法要比LK算法更接近真實(shí)值（圖1）。（3）ε_k與T_s – T_a存在顯著的負(fù)相關(guān)關(guān)系（y = -3.69 x + 18.16; r = -0.67，p < 0.01），這也說(shuō)明了對(duì)流湍流強(qiáng)度越強(qiáng)，動(dòng)力分餾效應(yīng)強(qiáng)度越弱（圖2）。（4）ε_k與水體面積無(wú)顯著相關(guān)性。小型水體的ε_k（¹⁸O）平均值（9.64 ± 2.80 ‰）與大型水體的平均值（10.10 ± 2.60 ‰）沒(méi)有明顯差別（p = 0.85, t = -0.20, 雙邊T檢驗(yàn)）。基于表1中的10個(gè)數(shù)據(jù)，ε_k(¹⁸O)的總體平均值為9.74 ± 2.60 ‰，對(duì)應(yīng)的湍流參數(shù)n為0.34，ε_k（D）則為8.5 ± 2.29 ‰。以上成果將為小型水體蒸發(fā)過(guò)程同位素動(dòng)力分餾過(guò)程提供科學(xué)依據(jù)，推進(jìn)氫氧穩(wěn)定同位素在內(nèi)陸水循環(huán)中的研究。研究成果以“Isotopic kinetic fractionation of evaporation from small water bodies”為題發(fā)表在《Journal of Hydrology》雜志上。南京信息工程大學(xué)肖薇教授為通訊作者，謝成玉博士為第一作者，南京信息工程大學(xué)劉壽東教授和張彌副教授、中國(guó)科學(xué)院地理科學(xué)與資源研究所王晶苑副研究員以及南京信息工程大學(xué)博士研究生張圳、胡勇博、劉強(qiáng)、碩士研究生錢(qián)雨妃和朱昊為合作作者。該研究得到了國(guó)家重點(diǎn)研發(fā)計(jì)劃和中國(guó)國(guó)家自然科學(xué)基金的支持。

來(lái)源：網(wǎng)絡(luò)