气象学的进展

PDF
气象学的进展/2020年/文章

研究文章|开放获取

体积 2020年 |文章的ID 6364531 | 22 页面 | https://doi.org/10.1155/2020/6364531

风速度无关两个源能量平衡模型基于一个理论梯形地表温度之间的关系和部分植被蒸散估算

学术编辑器:Hashiguchi欲之
收到了 2019年5月17日
修改后的 2019年11月09
接受 2019年12月10
发表 2020年2月14日

文摘

的准确估算陆地蒸散在异构表面使用卫星图像和几个气象观测仍然是一个具有挑战性的任务。风速(u),这是已知的具有高时空变异,是一个重要的约束在上述任务。在这项研究中,风速度无关两个源能量平衡(WiTSEB)模型的基础上提出了一个理论地表温度(Tr)部分植被覆盖率(fc)梯形空间和两阶段土壤水分蒸发蒸腾损失总量分解方法。理论上干燥的温度界限的Tr- - - - - -fc梯形迭代计算u通过使用一个假设没有明智的被水浸透的表面和大气之间的热交换在给定的大气条件下在垂直方向。WiTSEB HiWATER-MUSOEXE-12进行在黑河流域中游的八个景观利用ASTER影像。结果表明,WiTSEB提供可靠的估计在潜热通量(LE)和root-mean-square-errors (RMSE)和确定系数为68.6 W m−2和0.88,分别。的比值植被蒸腾的RMSE组件勒是5.7%。灵敏度分析表明WiTSEB不加重对气象和遥感输入的灵敏度与其他两个源模型进行比较。错误的估计Tr和观察到的土壤热通量的结果在LE高估或低估的部分风景。两级蒸发蒸腾分解方法是仔细验证了地面观测。

1。介绍

土壤水分蒸发蒸腾损失总量在建模中起着重要作用soil-vegetation-atmosphere陆地水文循环和能量交换的系统。卫星遥感技术潜在的区域和全球等提供了一个有效的方式由于常规和大空间尺度地表观测的属性(例如,反照率、地表温度、Tr,和植被指数,VI)。许多遥感模型意味着不同的理论提出了复杂性和假设1- - - - - -4]。其中,模型使用flux-profile关系(表示为巨大的温度梯度和气动阻力)的基础上Monin-Obukhov相似(MOS)理论和能量平衡理论物理计算勒(例如,SEBAL [5],度量[6],MOD16 [7,8],TSEB [9,他们…10],STSEB [11])或蒸发率(例如,实用12]和TTME [13]),主要是提出了。考虑到地表的不同待遇,这些基于物理方法可以进一步分为一个源的方法(例如,SEBAL、度量和实用和两个源的方法(例如,TSEB TTME,他们…,MOD16,和STSEB)。一个源方法治疗陆地表面作为一个“大叶”和坚持之间地表和大气湍流交换发生在一定的高度。虽然一个源算法报告可靠的热通量估计(6,14,15),他们往往需要一个精确的先验校准不能保证通过各种各样的表面条件,尤其是在稀疏的树冠地区(1- - - - - -3,16- - - - - -20.]。两个源模型对土壤和植被不同“来源”热量和水交换和模拟显热通量,H和LE组件组件的温度(即不同。,土壤和树冠temperatures), thereby representing an advancement to avoid the priori local calibration in the one-source model [9]。两个源模型被证明是强大的各种地貌和水文气象条件(20.- - - - - -24]。

一般情况下,风速u(或摩擦速度 )空气动力阻力的核心参数的计算,r一个——和两个源物理模型是基于flux-profile关系的基础上,金属氧化物半导体理论,和u(或 )大大影响H和估计。桑切斯et al。11强调10%的扰动u导致HSTSEB模型扰动和估计的17%和4%。长,辛格(13报道称,20%的扰动 带来的大约12.0%的变异勒估计TTME模型。20%的扰动u在他们…模型导致了7.4%的变异勒(10]。王等人。4增加或减少25%u导致8.6%的增加或减少11.5%H估计修改SEBAL模型。韦伯斯特et al。25)发现,uT一个那么有影响力吗Tr他们…和实用的估计。

然而,高质量的网格u(或 )不是经常。它不能被远程感知。尽管再分析数据集(例如,美国国家环境预报中心的能源部(NCEP / DOE)摘要采用国家大气研究中心(NCEP / NCAR),全球土地数据同化系统(GLDAS),和中国气象局的土地(CLDAS)提供网格数据同化系统u,精度基本上取决于地面观测网络,需要改进的匹配等估计的要求。Decker et al。26)表示,六再分析产品RMSE之间的1.5年代−1和4.5年代−1对33站在全球的观察。此外,再分析数据的时空尺度通常是太粗结合finer-scale遥感图像(27]。

有许多研究试图减少依赖u。一些研究人员改变了r一个来代替的u独立的抗性,μet al。7,8和姚明et al。28)认为,r一个是平行于辐射传输阻力和对流传输阻力和r一个定量关系T一个或者只和大气压力T一个。这种假设似乎过于简单模拟湍流通量的热量和动量的显著差异(29日]。别人介绍了一个策略来构建一个实际的像素和一个参考站点之间的桥梁。邱et al。30.假设实际的r一个等于在一个引用的裸露的土壤表面干燥。这种假设是不符合观测的刘et al。29日]在裸露的土壤和玉米树冠。太阳et al。31日随后提出了一个假设u引用的实际像素等于干裸露的土壤大气内表层修改邱的方法缺乏考虑空气动力学特征的地表。Nishida et al。27)采用类似的假设和发现的估计u是最大的误差源之一。的估计误差u推导了不稳定评估的变量(例如,Tr,Rn土壤热通量,G年代)引用裸露的土壤表面干燥。此外,相当于u假设不考虑表面粗糙度的影响u(25]。

为了克服依赖u,本文模拟气动阻力在植被覆盖的网站引用水饱和的假设缺乏合理的引用站点和上层大气之间的热交换在垂直方向,提出了u独立的两个源能量平衡(WiTSEB)模型。部分2介绍了WiTSEB模型。部分3描述了研究区和数据。部分4介绍了结果。部分5论述了灵敏度、误差来源、不确定性分析、与其他模型和比较。最后,部分6提供了本研究得出的结论。

2。材料和方法

WiTSEB模型包括三大模块,即。,两个源模型框架(部分2。1),分解Tr通过一个两阶段等方法(32,33(部分2。2),和一个迭代过程的计算c和勒年代(部分2。3)。第二个模块是主要基于一个理论u独立的Tr- - - - - -fc从一个梯形空间修改Tr王等人提出的vi梯形空间。4和王et al。34]。

2.1。两个源模型框架

两个源模型通常包含层配置(9)和补丁的配置(11,35]。层方法相比,补丁方法提供可比的估算准确率(24和有优势的要求更少的树冠结构的附加信息分配净辐射,Rn,土壤和树冠组件(11]。因此,WiTSEB采用补丁框架(11模拟能量通量: 在下标c年代分别代表了植被和土壤组件以后;Rn表示净辐射(W·m−2);H是显热通量(W·m−2);勒是潜热通量(W·m−2);GG年代土壤是土壤热通量和土壤热通量组件(W·m−2),分别;G年代可以估计分数(CG)Rns(36),CG变化从0.2到0.5取决于土壤类型和土壤水分,和一个常数值CG(= 0.35)和其他两个源模型(9- - - - - -11,13];fc是部分植被覆盖率和推荐的计算方法Choudhury et al。36),但取代了归一化差植被指数(NDVI)的增强型植被指数(以)因为NDVI在渐近结果(饱和)信号在高生物量条件[和扩展的问题7]。 在增强型植被指数马克斯和增强型植被指数最小值的以完整的植被和裸露的土壤表面,分别和n是一个经验系数,与叶树冠内定位和分布。

R数控Rns计算为 在哪里αcα年代代表植被反射率(无量纲)和裸露的土壤,分别;我们假设实际的像素α是一个加权组合的αcα年代fc,αc被固定为0.20 (11];Rd是向下短波辐射(W·m−2)和计算方法提出的艾伦et al。6];T一个空气温度(K);TcT年代树冠的温度和土壤(K),分别;ε是地表发射率(无量纲);ε一个大气辐射率(无量纲)和计算方法建议Brutsaert [37];和σ波尔兹曼常数。

HcH年代计算如下: 在哪里ρ是空气密度(公斤⋅m−3),Cp空气比热容是一个恒压(1004 J K−1公斤−1),r交流空气动力阻力(s m−1气象观测高度之间的传热z(m)和植被表面;它是表示为 在哪里公斤ydF4y2Ba是卡门常数(= 0.41);d是零点位移高度(m)和等于2/3呢h,在那里h植被的高度(m); 树冠表面的摩擦风速;lMonin-Obukhov长度(m);φhφ是热量和动量的稳定功能,根据计算的价值l;和z0hc是热的树冠粗糙度长度(m)和计算树冠动量粗糙度长度(z0mc)和无量纲参数( ),也就是说, ,在哪里 提出的方法计算Brutsaert [37]: 在哪里Cd的阻力系数是叶子元素和集到0.2 (38),Ct叶的传热系数,用户体验是一个无量纲阻力区域密度的函数。

r作为空气动力阻力(s m−1)之间的传热z裸露的土壤表面和计算 在哪里 土壤表面的摩擦风速在树冠和z0海关是热的土壤粗糙度长度(m)和土壤粗糙度长度计算的势头(z0女士), ,在我们组z0女士为0.01米; 计算(39] 在哪里是粗糙度雷诺数( , 空气的运动粘度)。

r党卫军空气动力阻力(s m−1)在树冠和土壤表面之间z0女士+d计算的方法推荐在曾庆红et al。40]: 在哪里c年代动荡的传递系数。

2.2。分解Tr基于梯形空间
2.2.1。建设理论Tr- - - - - -fc梯形空间没有u

王等人。34构建一个理论Tr(图六世梯形空间1(一)一个像素的基础上)。王等人。4)进一步细化Trvi梯形空间估计半干旱地区等的分水岭。从理论上讲,有四个顶点代表极端条件下在梯形空间;也就是说,Point 1 represents a well-watered vegetation in which the soil moisture of the root zone is sufficient and has a minimum canopy temperature,Tc、最小值;2代表的缺水的植物根区的土壤水分明显不足,最大冠层温度,Tc,马克斯;第3点所代表的是一个富饶的土壤有最小的土壤温度,T年代、最小值;和点4代表了裸露的土壤中,土壤蒸发干燥严重抑制,最大的土壤温度,T年代,马克斯。温度的四个顶点T年代vi梯形计算 在下标1、2、3和4代表顶点绘制在图的值1,Gf表示的比率GRn,Gf3Gf4分别设置为0.25和0.35,r厘米r残雪最小和最大树冠抗性(s m−1)和设置为12.5和625年,分别Δ饱和蒸汽压的斜率是空气温度(kPa C−1),γ是常数(KPa⋅C−1),VPD的蒸汽压赤字是空气(hPa)。

Trvi梯形方法具有强大的理论基础,但它需要网格u计算r作为r交流,这可能显示不确定性尤其在异构的风景4]。在这里,我们进一步修改Tr在两个方面构建一个vi梯形u独立的Tr- - - - - -fc梯形(图1 (b))。

首先,我们使用fc取代增强型植被指数在构建梯形空间的分解是一致的TrTr通常说的意思吗TcT年代加权的fc(9]:

第二,我们简化了计算极端温度采用一个假设显热通量在垂直方向有实力的优势在一定气象条件是微不足道的。它相当于以前的研究的假设13,41],边边缘理论忽略了水平平流时速度平衡等。这个假设意味着是一个富饶的中性大气条件优势金属氧化物半导体理论基础。然后,我们可以计算r交流2r作为4中性大气条件(例如,r交流0r作为0根据方程()16一个)和(16 c),而不使用u作为输入时Tc、最小值T年代、最小值是已知的。在实践中,我们采用平均T一个有实力的风景的Tc、最小值T年代、最小值。的计算r作为4(或r交流2在方程()13)(或方程(11)可以分为两个部分:r作为0(或r交流0)和大气稳定校正功能,f(φ女士4,φ海关4),或f(φmc2,φhc2)[29日]:

计算的过程Tc,马克斯T年代,马克斯一个像素的基础上总结了左边虚线框的图2,它分为两个主要步骤:(1)计算初始r交流2r作为4(例如,ther交流0r作为0通过方程(16一个)和(16 c),分别)f(φmc2,φhc2)= 0和f(φ女士4,φ海关4)= 0。(2)迭代计算Tc,马克斯T年代,马克斯(方程(16 b)和(16 d);R数控2(方程(8)),Rns4(方程(8 b)),Hc2(= 0.9R数控),H年代4(=Rns4(1−Gf4));lc2lc4;φmc2,φhc2,φ女士4,φ海关4; (结合方程(9)和方程(11)), (结合方程(10)和方程(13)), (方程(12)), (方程(14));和f(φmc2,φhc2)(类似于方程(方程18)),f(φ女士4,φ海关4)(方程(18)),r交流2(方程(11)),r作为4(方程(13)),直到r交流2r作为4是稳定的,即,d我fference between two adjacent calculations are smaller than five percent. Normally, the stability can be satisfied within 10 times.

应该指出的是,Hc2设置为0.9R数控2因为存在表皮角质层蒸腾在极端干旱条件下(42]。

2.2.2。计算组件的温度

我们采用两级等方法Tr- - - - - -fc梯形估计T年代Tc。两级等方法等与土壤水分的动态变化分为两个连续的阶段在一个给定的大气和植被条件下,即,土壤蒸发压力和植被蒸腾强调阶段阶段。

在土壤蒸发压力阶段,地表土壤水分(0 - 5厘米)逐渐减少,勒年代不同的极大值为0。相比之下,根区的土壤水分保留在绑定,可以维持一个树冠斯皮德这样近的潜在增长率。这些变化反映在Tr- - - - - -fc梯形(图1 (b)温度变化对fc下的三角形。T年代从最小值增加(T年代、最小值)到最大值(T年代,马克斯),Tc依然不变,给出的

在植被蒸腾强调阶段,根区土壤水分的大幅减少和抑制cTc从最小值不同(Tc、最小值)到最大值(Tc,马克斯)。此外,T年代维护T年代,马克斯在极端缺乏表面土壤水分。温度的变化与fc在这个阶段可以反映出上三角图中描述1 (b)T年代Tc可以计算为

2.3。计算勒年代,勒c,勒

上述假设也在这部分计算HcH年代,即,u独立的r交流0r作为0结合TcT年代被用来估计HcH年代分别由一个迭代过程。

计算的迭代过程HcH年代(图包含两个主要部分2,对虚线框表示):(1)计算初始r交流r作为,即,u独立的r交流0r作为0(2)迭代计算Hc(方程(9)),H年代(方程(10));lcl年代;φmc,φhc,φ女士,φ海关; (结合方程(9)和方程(11)), (结合方程(10)和方程(13)); (方程(12)), (方程(14)),f(φmc,φhc),f(φ女士,φ海关),r交流(方程(11)),r作为(方程(13));和r党卫军(方程(15)),直到r交流,r作为,r党卫军是稳定的

c和勒年代计算每个林冠下残余能量(方程(5))或土壤组件(方程(6))。勒估计为乐的意思c和勒年代加权的fc通过方程(3)。

3所示。研究区域与数据

3.1。研究区和地面观测

多尺度观测实验蒸散的异质性黑河流域2012年盟军遥测实验研究称为HiWATER-MUSOEXE-12 [43)是2012年5月和9月之间进行的沙漠绿洲过渡带在黑河流域的中游(图3)。年平均(1961 - 2010)的空气温度和降水的区域是7.4°C和128.7毫米44]。

HiWATER-MUSOEXE-12包含两个嵌套的观测矩阵,即大面积的30公里×30公里和核心区域位于盈科5.5公里×5.5公里,蹄兔灌区。HiWATER-MUSOEXE-12装备21自动气象站(AWS)和22涡度相关(EC)站在DM特大型发电厂(两个)和1水同位素站(DM特大型发电厂)在不同土地覆盖(包括玉米、沙漠、戈壁、荒漠草原,果园、蔬菜、住宅、湿地)。气象数据,包括T一个相对humdityrh向上/向下短波辐射,向上/向下长波辐射、净辐射、多层土壤湿度,土壤温度,观察21 AWS 10分钟的间隔。T一个rh是空间插值研究区与反距离加权法(IDW) [25]。土壤热通量G通过计算板卡尔法(45]。这个方法包括两个部分:吸热通量,计算使用fc加权平均的三个heat-plates(6厘米以下每个通量塔)周围的地面测量,以及在土壤蓄热的变化,计算了土壤温度变化速率,土壤湿度,土壤孔隙度在2和4厘米的深度。G不计算在S4, HZZ ZY网站由于缺乏土壤孔隙度测量。质量控制上G是主要的基础上进行理性的土壤温度的变化。植被高度,h测量,通常在每个AWS站的间隔h两个测量由线性插值计算。

EC的预处理通量平均在30分钟,通常分为三个质量等级如下:0级(平稳性的质量评估方法,Δst < 30,和湍流积分特性测试,ITC < 30),一级(Δst < 100和ITC < 100),和2级(Δst < 100和ITC > 100)。为了确保质量,EC通量从2级,从怀疑仪器漂移,从异常G(主要是指土壤温度测量异常)被否决。EC通量通常被认为是作为能量不平衡。波文比(BR)和残余能量(重新)是两种常用的方法用于解决能量不平衡问题在电子商务系统46]。在目前的研究中使用的BR方法。DM站,比勒(LE蒸腾的组件树冠/ LE,勒树冠是树冠蒸腾的组分在混合像素,=勒c×fc勒(LE)和蒸发组件土壤/ LE,勒土壤是土壤蒸发组件在混合像素,=勒年代×(1−fc)是使用腔衰荡光谱测量的(crd)水蒸气同位素(模型L1102-I、Picarro Inc .)和代表平均13:00-15:00(当地时间)44];等及其蒸发组件通过梯度扩散计算方法和。戈登模型,分别。

3.2。遥感数据

ASTER(先进星载热发射和反射辐射仪)一级b图像八万里无云的日期(6月24日6月15日,7月10日,8月2日,8月11日,8月18日,8月27日和9月3日)获得来自美国地质调查局网站(https://earthexplorer.usgs.gov/)。ASTER为可见光和近红外波段的观测提供了空间分辨率15 m 12:10 - 12:20(当地时间)。ASTER L1B可见红外数据重新取样90通过最近邻方法。α通过拟合计算窄带(相对于红光47]: 在哪里ρ1,ρ2,ρ3代表的相对于红光band1、二级和band3覆盖了可见光和近红外光谱,分别。通过双波段的增强型植被指数计算方法提出的江et al。48)认为ASTER传感器没有提供一个蓝色乐队:

土地覆盖图像是通过视觉获得解释。Trε产品提供的寒冷和干旱地区科学数据中心;这些产品估计temperature-emissivity分离(te)算法的基础上(49)使用ASTER热红外波段的空间分辨率90米。

4所示。结果

4.1。验证WiTSEB输出

WiTSEB输出的散点图(Rn,G,H,LE)与地面观测显示在图中4。WiTSEB输出的准确性进一步用均值偏移误差(MBE),均方根误差(RMSE),和确定系数(r2)与地面观测相比,表中列出1


通量 LUCC 观察到的平均(W m−2) 估计平均(W m−2) (W m MBE−2) RMSE (W m−2) r2

Rn 住宅 550.0 581.8 31.8 37.1 0.82
戈壁 474.7 496.1 21.4 28.6
蔬菜 640.3 614.4 −25.9 34.5
湿地 617.7 682.3 64.6 65.1
沙漠 512.3 517.2 4所示。9 35.1
沙漠草原 507.6 522.8 15.2 20.1
果园 730.2 673.9 −56.3 58.1
玉米 659.5 648.3 −11.2 23.6
整体 631.3 626.0 −5.3 31.9

G 住宅 - - - - - - - - - - - - - - - - - - - - - - - - 0.22
戈壁 135.4 149.7 14.3 19.4
蔬菜 84.4 66.7 −17.7 38.9
湿地 NA NA NA NA
沙漠 129.1 173.9 44.8 50.8
沙漠草原 - - - - - - - - - - - - - - - - - - - - - - - -
果园 108.7 110.1 1。4 38.3
玉米 86.9 73.6 −12.7 45.9
整体 94.3 87.9 −6.4 43.1

H 住宅 230.3 142.7 −87.6 98.5 0.73
戈壁 263.1 327.1 64.0 86.8
蔬菜 121.3 131.1 9.8 49.9
湿地 79.6 28.5 −51.1 62.5
沙漠 292.2 359.1 67.1 84.3
沙漠草原 318.0 357.0 39 69.9
果园 122.2 49.3 −72.9 85.4
玉米 85.3 69.4 −15.9 50.7
整体 131.2 119.5 −11.7 63.5

住宅 246.4 281.3 34.9 62.5 0.88
戈壁 75.8 33.7 −42.1 66.7
蔬菜 438.1 399.7 −38.4 115.1
湿地 458.1 528.8 70.7 76.6
沙漠 90.9 17.0 −73.9 95.0
沙漠草原 136.4 28.1 −108.3 118.7
果园 509.9 517.6 7.7 35.6
玉米 484.3 502.6 18.3 56.4
整体 405.4 408.1 2。7 68.6

Rn最好的建模是基于方程(1),(8)和(8 b)的总体RMSE 31.9 W m−2在四个组件的地表能量平衡方程。有智商有点低估了5.3 W m−2与一个r20.82。估计的性能Rn随风景。它低估了玉米,蔬菜,和果园景观,但被高估了在沙漠草原,戈壁、住宅和湿地网站。最大的错误是发生在湿地景观与MBE 64.6 W m−2和65.1 W的RMSE m−2,这将是可能的严重低估α

由于缺乏土壤孔隙度测量、估计的准确性G只有在玉米、蔬菜、果园、沙漠和戈壁的网站。G估计常数Rns(= 0.35),整体−6.4 W m的英帝国勋章−2,43.1 W的RMSE m−2和一个r20.22,相比之下tower-based测量。固定比例策略可能的低性能的主要原因G估计。

合理的估计之间的协议是获得HMBE和tower-based观测和RMSE−11.7和63.5 W m−2,分别。蔬菜和玉米网站相对较小的偏差H估计在对比其他的风景。负H测量是在几个人玉米网站(图4 (c)),而他们的H估计是积极的,因为假设是缺乏会计水平对流的影响。

估计2.7 W m WiTSEB勒的英帝国勋章−2,68.6 W的RMSE m−2和一个r20.88(表1)。玉米和果园网站显示比干更紧密的协议与地面观测站点,即。沙漠、戈壁、沙漠和草原网站,使用更小的偏差和RMSE较低。

此外,我们也估计勒(我们叫它TSEB模型u为了方便WiTSEB)相同,即,p一个tchfr一个米e and the two-stage decomposing method, except that TSEBu员工u作为输入来计算空气动力学抗性(例如,r交流2,r作为4,r交流,r作为,r党卫军)。结果在表2LE-estimated TSEB的准确性u利用ASTER影像从八个日期HiWATER-MUSOEXE-12 21塔观测。尽管TSEBu在蔬菜方面优于WiTSEB RMSE,湿地,沙漠,沙漠和草原,它显示了稍低的整体性能,产生一个MBE 8.0 W m−2,70.7 W的RMSE m−2和一个r20.87。


精度 住宅 戈壁 蔬菜 湿地 沙漠 沙漠草原 果园 玉米 整体

观察到的 246.4 75.8 438.1 458.1 90.9 136.4 509.9 484.3 405.4
估计 258.7 36.0 445.6 499.5 26.4 71.8 528.3 506.5 413.4
MBE 12.3 −39.8 7.5 41.4 −64.5 −64.6 18.4 22.2 8.0
RMSE 66.6 74.7 85.9 47.0 83.6 102.0 38.0 68.2 70.7

4.2。与先前的研究相比

HiWATER-MUSOEXE-12数据集被用来评估不同评估方法在之前的研究(表3)。通过不同的方法之间的不同的RMSE LE 17.9和133 W−2,中位数为74.8 W m−2。在那里,TD-TSEB [28)和非参数模型(51)移除依赖的优势u提供一个RMSE 89.8和133 W m−2,分别。WiTSEB使用相似的输入和两个源框架作为TD-TSEB但RMSE较低(即。41.9 W m−2DM站)。虽然这精度比较缺乏绝对的意义,因为假设不同,它可以说明可接受的精度。


作者 模型/方法 输入 网站 图像(日期) 精度
地面 遥感产品/同化 RMSE MBE r2

姚明et al。(28] TD-TSEB T一个,P副总裁(蒸汽压) α,ε,Tr、归一化植被指数 DM 陆地卫星(5 6月一日至9月30日) 89.8 −63.6 0.81
李等人。50] 实用 T一个,P,u特定的湿度,Rd,Rld α,ε,Tr,fc,hc赖, S1-S17、ZY GB,量HZZ ASTER(5月30日,6月15日,6月24日,10 7月,8月2日,8月8日,27日,8月和9月3日) 64年b
63年r
24b
−8r
锅等。51] 非参数 Tr,α,ε归一化植被指数,T一个,P,rh,理查德·道金斯 ZY, S1, S4,肌力,GB, f MODIS(6月25日至9月9日) 133年 −59 0.48
徐et al。52] TVDA u,T一个,rh,P,Rd,Rld,Tr赖, DM (6月7日到9月4日) 116年 46
TSEB 132.9 82.2
歌等。53] 两个来源+贝叶斯 Rd,Rld,u,T一个,rh,P,hc Tr,α,ε归一化植被指数,赖 S1-S17 ASTER(10 7月和8月11日) 17.9 0.7 0.98
歌等。47] TSEBCT T一个,u,P副总裁,Rd,ε 归一化植被指数,Tr,α S1-S17 ASTER(6月15日、6月24日10 7月,8月2日,8月11日,8月18日,8月27日,9月3日,9月12日) 52b
58r
3b
−21r
TSEB T一个,u,P副总裁,Rd,ε 归一化植被指数,Tr,α 87年b
71年r
50b
25r
杨et al。54] 他们… T一个赖,hc,u,fc,rh、磷、VPD 归一化植被指数,α,εTr S1-S17、ZY GB,量HZI ASTER(10 7月,8月2日,8月11日,8月18日,8月27日,9月3日) 42.3r −15.7r 0.88
TSEB 49.8r 4所示。8r 0.85
MOD16 95.3r −34.5r 0.46
黄等。55] SEBS-NDWI T一个,u,P,Rd Tr,α,ε归一化植被指数,NDWI S1-S17、GB、量、ZY HZZ MODIS(6月15日到9月5日) 84.1 37.6
实用 T一个,u,P,Rd副总裁 Tr,α,ε、归一化植被指数 117.8 69.4
马等。56] 实用 T一个,u,P,rh Tr,α,ε归一化植被指数,赖 S1-S17、ZY GB,量HZZ ASTER(6月15日、6月24日10 7月,8月2日,8月11日,8月18日,8月27日,9月3日,9月12日) 102.9b 68年b 0.77
修改后的实用 74.8b 23b 0.81
壮族和吴57] 改善TSEB u,T一个,rh,P,Rd,Rld Tr,α,ε、归一化植被指数 S6-S14 ASTER(8月11日,18 8月、9月3日和9月12日) 35.1b 4所示。1b
歌等。58] TSEB 赖、植物宽度、株高、叶片宽度、叶片长度,T一个,u,rh,fc,Tr DM (6月18日,6月21日,6月24日,7月10日,7月13日,8月18日,9月8日,9月12日) 61年 31日
朱et al。59] Shuttleworth-Wallace +贝叶斯 赖,T一个,u,rh,Rd,Rld、土壤水分 DM (5月26日至9月13日) 80.7 24.2 0.83

bc代表观察到的勒通过BR和纠正方法,分别。
4.3。LE分区

五天(6月24日,7月10日8月2日,8月11日和8月27日)在八个日期的ASTER crd水蒸气同位素观察DM站。为期五天的树冠蒸腾组件的平均比率(LE树冠/ LE)估计为84.6%,低于3.2%的观察(表4)。一项研究表明,crd略有高估了勒的比率树冠/ LE [44]。因此,勒的准确性树冠/ LE估计通过WiTSEB有望高。相比之下,歌曲等。58)采用TSEB模型使用一个eight-date观察Tr从一个侥幸Ti55热红外成像仪和其他观察结果(表3)来估计勒树冠和RMSE MBE / LE和获得了1%和2%,分别。杨et al。54)评估他们…的表演,TSEB和MOD16模型使用ASTER影像从六个日期和MBE决定了勒树冠/ LE−1.1%、4.1%分别和−20.5%。姚明et al。(28)使用陆地卫星图像从五个日期评估TD-TSEB MBE模型并获得−11.1%。


项目 观察到的平均 估计平均 MBE RMSE

树冠/ LE (%) 87.8 84.6 −3.2 4所示。6 5。7
土壤/ LE (%) 12.2 15.4 3.2 4所示。6 5。7
树冠(W m−2) 469.3 468.4 −0.9 43.9 49.3
土壤(W m−2) 63.8 85.7 21.9 25.3 33.7
勒(W m−2) 533.1 553.9 20.8 35.3 41.9

我们评估的性能WiTSEB(表4在勒土壤(土壤蒸发组件在混合像素,=勒年代×(1−fc))和勒树冠估计在DM火车站。结果表明,MBE (RMSE)的勒土壤,勒树冠,勒是21.9(33.7),−0.9(49.3),和20.8 (41.9)W m−2,分别。相比之下,歌曲等。58)获得的RMSE勒在61 W m−2通过使用TSEB模型和地面观测。勒的值土壤,勒树冠和LE TD-TSEB模型低估了的28)分别为4.2,60.1和63.6 W m−2的RMSE 16.5, 91.3,和89.8 W m−2,分别。

4.4。时空变化

估计LE地图展现在图8日期5。这些地图显示显著的空间差异主要由土地覆盖的空间分布。一般来说,勒大高的地方fc(例如,wetland, corn, vegetable, and orchard) and low at places with sparse vegetation (i.e., Gobi, desert, and desert steppe). Wetland has the largest LE given sufficient water supply, with an average of 485.4 W m−2在八个日期。虽然降水稀少(2012年6月和9月之间的累积沉淀∼100毫米),果园,玉米,和vegetable-covered区域产生一个高勒因为灌溉,平均值为452.1,422.1和393.3 W m−2,分别。此外,区域土壤水分灌溉可以放大空间差异增加空间变化在LE尽管相同的土地覆盖类型;例如,勒在6月24日东北地区较高的核心区域(蓝色框)比在其他领域(图5(左列)。土壤水分在4厘米深度明显高于在S3, S9, S10(36 - 39%)比其他corn-covered网站(17 - 29%)。居民区的LE平均为262.6 W m−2,主要来源于树冠蒸腾(174.6 W m−2)。由于水压力(即。,lack of precipitation and no irrigation), desert, Gobi, and desert steppe provide a small LE, with an average of 44.5, 66.3, and 69.0 W m−2,分别。8月18日是唯一和为期两天的前期有效降水(即约会。,每日3 - 5毫米的降水大于0.5毫米);因此,戈壁、沙漠和荒漠草原始终有一个明显高于勒比其他日期8月18日。

LE地图的比较在不同的日期表示时态的变化勒估计在高的地区fc由于物候变化。勒在玉米和蔬菜区域的平均值明显高于6月7月和8月和9月。果园和湿地地区也显示时间变化,但他们的范围较小比玉米和蔬菜。

在图6,勒的比率树冠勒(LE树冠/ LE)在沙漠和荒漠草原的表面的最高平均0.95和0.92,分别在8个日期。这个结果是由于土壤蒸发是极其强调考虑到严重的干燥土壤表面,而树冠可以使用深层土壤水分蒸腾通过维护一个轻微的蒸腾。勒树冠/ LE沙漠和荒漠草原的表面低于在其他日期8月18日给土壤蒸发与降水的增加。作物(包括玉米和蔬菜)和果园地区类似的勒树冠/ LE平均约为0.86。住宅最低勒树冠/ LE平均为0.72。

5。讨论

5.1。敏感性分析

WiTSEB模型输入包括以下细节:(1)地面气象观测:T一个,rh,h;(2)遥感数据:Tr,α,ε,增强型植被指数;和(3)的参数:Gf3,Gf4,α1,α2,α3,α4,z0mc,r厘米,r残雪。输入模型灵敏度分析理解误差的来源和误差传播机制WiTSEB。勒的敏感性th输入,年代,被定义为 在哪里勒0代表平均勒估计实际输入,勒±代表的平均LE 8个研究区域估计的日期th输入增加(+)或减少(−)的比例值或对实际的勒,和变化范围和步骤2 K和0.5 K温度输入(T一个Tr其他输入)和20%和5%,分别。

在表5温度变量,WiTSEB敏感。TrT一个分别是消极和积极与勒。2 (1)K减少或者增加Tr可能会导致一个15.9(8.1)17.2(8.5)% %增加或减少勒估计,分别。−2的扰动(−1)K或2 (1)KT一个将导致LE扰动−22.7(−10.8)%或18.9(9.9)%,分别。WiTSEB显示类似的灵敏度大小其他遥感变量(α,ε,增强型植被指数)和低于敏感Tr。勒估计是积极的增强型植被指数,以增加或减少20%导致增加或减少4.5% 1.9%勒估计,分别。估计是负面的α,ε,rh。增加或减少20%α将导致减少了5.1%或5.3%的增长。WiTSEB构成的参数。LE小于3%时的摄动参数扰乱20%。勒估计显示出积极的关系Gf3,α1,α3,z0mc,r残雪但负相关Gf4,α4,h,r厘米。−二十(20)%的扰动α2,α4,h会导致LE扰动0.1 (−0.1)%,1.9 (−2.3)%,0.2 (−0.4)%。


变化(K / %) 20 (−−2) 15 (−−1.5) 10 (−−1) −5 (−0.5) 5 (0.5) 10 (1) 15 (1.5) 20 (2)

Tr 15.9 12.1 8.1 4所示。1 −4.2 −8.5 −12.8 −17.2
T一个 −22.7 −16.6 −10.8 −5.3 5。1 9.9 14.5 18.9
α 5。3 3.9 2。6 1。3 −1.3 −2.6 −3.9 −5.1
rh 3.7 2。9 2。0 1。0 −1.1 −2.4 −3.7 −5.2
增强型植被指数 −4.5 −3.2 −2.0 −0.9 0.8 1。3 1。7 1。9
ε 6.2 4所示。6 3.0 1。5 −1.5 −3.0 −4.4 −5.8
Gf3 −2.5 −1.9 −1.3 −0.6 0.7 1。3 2。0 2。7
Gf4 2。7 2。1 1。4 0.7 −0.8 −1.7 −2.6 −3.5
α1 −0.3 −0.2 −0.1 −0.1 0.1 0.1 0.2 0.3
α2 0.1 0.0 0.0 0.0 0.0 0.0 0.0 −0.1
α3 −1.9 −1.4 −1.0 −0.5 0.5 1。0 1。5 2。0
α4 1。9 1。5 1。0 0.5 −0.5 −1.1 −1.7 −2.3
z0mc −0.6 −0.4 −0.3 −0.1 0.1 0.3 0.4 0.5
h 0.2 0.1 0.1 0.1 −0.1 −0.3 −0.4 −0.4
r厘米 0.2 0.2 0.1 0.1 −0.1 −0.1 −0.2 −0.3
r残雪 −0.2 −0.1 −0.1 0.0 0.0 0.0 0.1 0.1

的变化T rT 一个在K和其他变量的百分比(%)。

类似的敏感性分析进行对许多trapezoid-based两个源模型(表6),比如他们…[10],TTME [13],ESVEP [33]。与这些模型相比,WiTSEB避免使用u作为输入,不增加其他投入的依赖。


模型 Tr(K) T一个(K) α α2 α4 h
−2 (−1) 2 (1) −2 (−1) 2 (1) −20% 20% −20% 20% −20% 20% −20% 20%

TTME 28.9 (14.4) −28.6 (−14.3) −27.7 (−13.9) 27.6 (23.8) 1。4 −1.4 4所示。9 −5.7 1。5 −1.5
他们… 17.3 (9.5) −23.2 (−11) 21 (−−10) 15.3 (8.5) 7.1 −7.0 0.73 −0.75 2.02 −2.04 2。4 −2.1
ESVEP (6.5) (−7.7) (−6.9) (5.7) 4所示。7 −5.0 −0.5 0.5

此外,如果地面RdRld观察是足够的,我们可以使用它们作为输入WiTSEB。这种方式,WiTSEB可以相对的性能改善,即。的RMSE LE减少从68.6 W m−262.9 W m−2r2从0.88增加到0.89。我们进行了灵敏度分析。这两个RdRld勒估计是积极的。5%的扰动RdRld勒估计变异可能导致6.6%和3.1%,这几乎是一样的ESVEP [33]和略低于STSEB [11]。

5.2。误差源分析
5.2.1。错误Trα估计

Tr观察是计算 在哪里RRld是观察到的向上和向下长波辐射(W m−2),分别。陆地表面发射率ε来源于Tr产品,是估计的算法。

在图7(一),Tr产品与观察,显示一致性,MBE美,RMSE,r22.9 2.0 1.1 K, K, K,分别和0.89。的Tr产品显示在玉米和果园面积精度高。同样,LE估计在这两个领域产生高精度(图4 (d)和表1)。Tr高估了在戈壁、荒漠草原和沙漠地区,MBE 2.1, 4.6,和5.4 K,分别。勒估计稀疏植被区域的负面偏见42.1,108.3和73.9 W m−2,分别。6月24日6月15日,8月2日,8月11日,8月27日和9月3日Tr在沙漠草原网站是高估了4.5,4.6,6.0,3.4,5.6,和5.6 K,分别,所以可能会呈现点(fc,Tr)信封显示在图1 (b)和导致勒估计约为0。在住宅的网站,Tr大大低估了,这是符合34.9 W m的LE高估−2呈现在图4 (d)和表1

类似于Tr,α也估计从ASTER影像。估计α有0.021的RMSE和一个吗r20.65。α一般显示轻微高估与AWS(图的观察7 (b)0.01)与英帝国勋章。然而,α给大大低估了在湿地,−0.05与英帝国勋章,这创造了相对高估的Rn(图4(一)和表1)和勒(即。,一个nMBEof 70.7 W m−2),尽管Tr是高估了(即。0.64 K),英帝国勋章。

5.2.2。错误T一个插值

我们计算出标准差(SD) 90×90分辨率的网格量化空间可变性(表7),采用分析crossvalidation (LOOCV)方法分析插值精度(图8)。0.1 K和0.3 K之间的SD不同平均为0.2 K在八个日期,从而表明较低的空间变异性。内插的T一个显示了观察的良好的关系r20.97。内插的T一个低于观察吗T一个0.2 K和有一个美和RMSE 0.4 K和0.5 K,分别。这两个指数暗示T一个可能不是的主要误差源WiTSEB模型在这一研究领域。


日期 最小值 马克斯 的意思是 SD

6月15 297.7 299.0 298.4 0.1
24日6月 297.8 300.2 299.1 0.3
7月10 299.1 300.9 300.0 0.1
8月2 300.5 302.8 301.8 0.3
8月11日 297.9 299.2 298.6 0.1
8月18日 296.3 297.5 297.0 0.2
8月27日 298.6 300.7 299.8 0.3
9月3日 291.4 293.0 292.3 0.2

然而,研究人员时需要谨慎使用方法复杂的地形和罕见站区域因为插值方法的性能敏感的气象数据的密度和可变性。韦伯斯特et al。25)评估七个插值方法在不同的地形和景观异质性高澳大利亚东南部地区,发现没有一个空间插值方法可以提供一个可靠的性能在各种条件和气象数据。

5.2.3。引入的错误假设

规避的依赖u一个假设,即。,water-saturated surface is under neutral conditions, is used to construct theTr- - - - - -fc梯形(图1 (b)),通过平均计算LE组件T一个蔬菜、玉米和湿地景观。使用一个特定的条件T一个作为Tc、最小值T年代、最小值是一种常见的实践构建一个理论吗Trvi或Tr- - - - - -fc梯形空间(10,13,27,32]。唐,李33发现估计Tc、最小值T年代、最小值通过莫兰的方法(60略有不同T一个在禹城火车站被农田覆盖在58晴空天从2010年3月到2011年10月。王等人。4发现计算Tc、最小值通过修订莫兰低于像素的方法T一个2.9 K,然而,这种差异与半干旱植被条件。然而,王子和Goward [61年)认为,两个温度之间的差异的大小大约是2 K。

在这里,我们分析的误差传播WiTSEB变量(例如,R数控1,Rns3,r作为0,r交流0,r作为,r交流,T年代,马克斯,Tc,马克斯,T年代,Tc,Hc,H年代,H和LE)通过引入的减少Tc、最小值,T年代、最小值在3 K,通过考虑王等人的研究。4]。在图93 K下降T年代、最小值导致变量的大小增加或减少不到4.5%,大小是0R数控1,r交流0,r交流,Tc,马克斯估计。3 K下降Tc、最小值而在Tc、最小值T年代、最小值结果在增加Hc估计分别为15.1%和16.3%,而对其他变量的影响并不明显,5.5%和6.9%H估计在勒估计和−−3.3%和4.3%,分别。

5.3。的不确定性G

EC通量测量广泛被称为一个能量不平衡,即。观察通量之和(H+ LE)小于可用的能源(RnG),从而产生了大约16%的不确定性在HiWATER-MUSOEXE-12勒(62年]。原因导致能源unclosure包括测量不匹配的足迹63年)、水平平流和忽视蓄热测量G(64年]。这里,我们着重分析G对勒观察测量误差及其影响。板卡尔的方法(45用于观察G。网站(包括S4住宅、HZZ沙漠草原和ZY湿地)不用于计算蓄热变化缺乏土壤孔隙度的测量,从而导致G观察不到实际值。因此,纠正勒通过BR方法将高于实际价值。HZZ沙漠草原有一个系统性的错误更少Tr估计(4.6 K和5.4 K MBE),但有一个更大的LE低估(−108.3 W m−2与−73.9 W m−2比沙漠站点)。这种对比可以解释观察到的错误造成的勒G观察,即。,underobservedG增强BR-corrected勒和加剧勒估计的低估。此外,underobservedG在剩余站点减轻了过高的估计严重低估造成的Tr

图所示的结果4 (b)和表1表明估计之间的相对较大的分歧G和测量G。这个分歧是可能由于固定值CG在整个研究区域的所有8天。CG据报道,随时间变化,土壤类型和土壤水分(36]。具体地说,CG显示昼夜变化(65年),但通常是恒定在当地中午太阳周围的时间范围(66年)和时间范围内的ASTER天桥时间;因此,设置的策略CG在给定的常数土壤类型和土壤水分条件一天是可以接受的。有系统地高估了G在戈壁和沙漠网站,对比蔬菜站点的性能。这种差异可以解释为土壤性质差异(54),即。,dry rock and sandy soil has lower soil heat capacity than vegetable soil due to the lower soil moisture and higher soil porosity [67年]。土壤水分的影响CG在这项研究中,很难定量评估因为我们缺少地面数据隔离土壤特性的影响玉米唯一景观显示统计学意义。重要的是,当地的校准CG研究区域的网站将提高协议;然而,本研究的目的是建立一个广泛的应用u独立的两个源能量平衡模型,而不是那些特定的网站优化。

5.4。的比较Tr分解方法

的分解TrTcT年代两个源模型是一个关键过程。四种直接方法,即multiangle方法(68年),经验法(69年),湿润等值线法(70年),两级等方法(32,33),采用;此外,一个间接法,即Priestley-Taylor迭代法(9),使用。鉴于缺乏两个视图角度的图片,在这里,除了multiangle方法,我们主要讨论三个直接的方法,即。经验,湿润等值线,两级等方法。

实证方法使用之间的经验关系T年代TcTrT一个(例如,T年代Tc=C一个(TrT一个),在那里C一个= 0.1,= 2)[69年)段Tr。我们之间的关系TcTr在一个给定的气象条件下(T一个基于Lhomme = 293 K)的关系。在图10,Tc近似线性增加Tr在高植被(fc> 0.7),而Tc随的增加而减小Trfc< 0.7条件。这种现象是自然状态,相反的TrTc增加与水压力和植被蒸腾约束。此外,Tc甚至小于0的情况下Tr> 315 K(例如,TrT一个在稀疏地区> 20 K)。事实上,平均TrT一个在ASTER天桥八日期时间在沙漠,沙漠草原和戈壁地区为20.9,18.2,和18.1 K,分别。因此,Tc错误地估计了在这些地区。这些现象验证实证关系的限制,特别是在稀疏区域(71年]。此外,Lhomme经验系数随地区而异。詹et al。71年)报道,系数C一个在0.07和0.27之间波动设置为2基于测量数据从90年的87年,季风的横笛和沃希托河的92个站点。

湿润的等值线法最初提出PCACA模型(70年)的假设几乎直等压湿润与等价的T年代Tc在一个统一的大气环境和土壤表面均匀。这种方法有相似的两级使用梯形空间段等方法Tr。然而,两种方法之间的分解组件温度大大不同。在这里,TcT年代分解的湿润等值线法表示的Tr- - - - - -fc梯形空间Tc,张T年代,张做一个区分。在图11,Tc几乎等于T一个平均为297.9 K DM站在八个日期。相比之下,Tc,张大于Tc3.5 K。卡尔森(41)宣称,Tc主要是由根区土壤水分。吴et al。72年)的深度报道,玉米根区是20 - 80厘米生长期。这里,我们计算了在根区土壤水分的平均土壤湿度20岁,40岁,80厘米,发现玉米根区土壤水分是略大于所有日期的田间持水量,从而表明植被蒸腾并不强调。因此,Tc理论上约等于T一个。的平均值T年代307.7 K和明显高于吗T年代,张5.6 K。T年代的主要原因是土壤表面水分(∼2厘米),考虑到吗Tr在实际是“皮肤温度”深处,给定的波长的电磁辐射可以穿透(73年]。

此外,显著差异也体现在T年代Tc。的平均值T年代Tc在八个日期是9.7 K,T年代,张Tc,张只有0.8 K。6月24日重要水分胁迫发生与土壤表面的水分含量只有10.5%。T年代理论上必须明显大于Tc;然而,T年代,张Tc,张只有2.2 K(图11),这可能是相对较小而观察。先前的研究已经报道,观察到T年代Tc可能超过10 K。例如,Colaizzi et al。74年)发现,T年代Tc中午结束10 K,赛季中期初在陆地棉区域,原始林区,德克萨斯州,美国根据观测的红外热成像(模型红外热成像/ c, Exergen Corp .),水城,麻萨诸塞州,美国)。田et al。75年)报道,大豆T年代Tc有超过20 K的极值中午利用热红外成像仪(福禄克红外FlexCam Ti55,侥幸作物。美国)。一些研究人员观察到T年代Tc经常超过20 K在半干旱的环境中(例如,Chehbouni et al。76年),休谟et al。77年],Kustas et al。78年])。Kustas和诺曼22)强调,平均水平T年代Tc有一个极端的25 K的情况下。艾伦et al。79年)表明,T年代Tc可以在公开的20 - 30 K和干燥表面。

等方法是基于湿度等值线估计土壤蒸发分数(LE年代/ (RnsG年代(LE)和植被蒸腾分数c/R数控假设他们拟线性变化T年代Tc分别在梯形空间;这些方法包括PCACA [70年],TTME [13,估计潜在蒸发分数(LE年代/勒sp,勒sp是潜在的土壤表面蒸发)和植被潜在蒸腾分数(LE吗c/勒cp)拟线性T年代Tc分别(例如,ETEML(杨et al。80年])。根据我们的假设,我们可以推断出粗糙的蒸发和蒸腾分数之间的平等33)和土壤之间的平等和植被潜在蒸发比率。在这里,我们计算了使用观察LE潜在蒸发和植物蒸腾作用的分数土壤/乐、乐、和气象数据(例如,T一个,rh,uDM站评估扣除),观察到勒年代计算了勒土壤/勒×勒/ (1−fc);fc通过计算安德森等人提出的方法。81年),勒sp和勒cP使用公式计算推荐的μet al。8和关和威尔逊82年]: 在哪里r年代,r一个,r合计土壤表面的气动阻力,树冠表面,和总,分别和r年代、最小值是最低表面电阻。

在图12,潜在蒸发部分显然是低,波动幅度大于潜在的蒸腾分数。潜在蒸发分数范围在0.09和1.32之间,平均和标准偏差为0.47和0.24,分别。值得注意的是,潜在蒸发分数大于1 7月30日主要是由于小EsP从一个高了rh值(79%)。潜在蒸散分数在0.81和1.37之间波动,平均和标准偏差为0.95和0.1,分别。此外,潜在的蒸腾分数在灌浆期期(约在8月)平均为0.98,高于在成熟时期(大约8月和后),为0.91。因此,湿润等值线法有可能高估了勒土壤和低估了勒树冠。的MBE勒土壤和勒树冠通过湿润等值线法是53.4−11.1 W m−2分别,从而证实了结论。

两个先决条件是隐含在两级等方法。第一个前提是在根区土壤水分减少落后,表面上(∼2厘米)。第二个前提是土壤水分变化的司机辐射温度等变化。图12说明表面土壤湿度显示了戏剧性的减少漫灌的两个实例之间的间隔(DM展出漫灌7月2日,7月28日和8月28日44])。在根区土壤水分减少在间隔和维护的范围大于田间持水量,从而表明没有水的压力。Tc被认可和广泛使用作为水的可用性指标的计算等作物水分胁迫指数(CWSI) [83年]。莫兰et al。60]发达CWSI裸露土壤的概念并提出了水分亏缺指数,表明T年代大约是线性相关的表面土壤水分和土壤蒸发在给定的T一个,Rn、VPD和r一个验证,比达尔和毕雷矿泉水84年]。因此,应用两阶段适合分段等方法Tr在一个给定的气象条件。

6。结论

减少依赖地面气象数据,特别是数据称为高空间和时间变化,即。、风速u,显著减轻不确定性的遥感等大规模和异构表面。在这项研究中,我们提出一个WiTSEB模型通过使用一个简化的理论表面温度(Tr)植被覆盖率(fc两级等)梯形Tr分解方法,和两个源补丁框架。

WiTSEB在HiWATER-MUSOEXE-12网站开展的沙漠绿洲过渡带Zhangye城市黑河流域中游的使用ASTER影像在八个日期。Rn,H和估计,从21通量塔观测结果吻合较好。RMSE和r2勒估计68.6 W m−2和0.88,分别。的精度Rn,G,H,勒估计显著不同土地覆盖类型(包括玉米、沙漠、戈壁、荒漠草原,果园、蔬菜、住宅、湿地)。MBE和RMSE勒树冠/ LE估计−3.2%和5.7%,分别比crd水蒸气同位素测量。勒估计显示高跨景观空间变化。勒估计通常是大型高植被覆盖区域和低稀疏植被区域。此外,灌溉策略放大空间的差异。勒估计在农田(玉米和蔬菜)高7月和8月,6月和9月低。勒估计在果园和湿地区域显示6月和8月之间的时序变化更小,但是明显大于9月的价值。降水的主要原因是勒在低植被时空变化的区域。

WiTSEB是最敏感的T一个T年代但对其他气象、遥感和其他参数。增加了2 KTrT一个导致降低17.2%和勒估计,增加了18.9%。其他投入引起的扰动敏感性20%,达到6.2%。与其他两个源模型相比,WiTSEB不加重气象和遥感输入的灵敏度。误差分析表明,估计误差Tr有助于LE高估或低估在某些土地覆盖类型。观察到的误差G会加剧LE低估和减轻LE高估。理性的评估表明,假设,这是用来避免u在构建Tr- - - - - -fc梯形空间和迭代计算H核心组件,收益率最小影响变量的估计。

一个两阶段的合理性等模式分解Tr是全面验证了使用地面气象观测的DM站。相比经验和湿润等值线,两级等模式提供了更大的灵活性的自然环境,给勒估计精度高。

与其他卫星模型,WiTSEB有一定的限制,因为一些不经常使用输入,如树冠高度(h)[1),这是一个必要的因素对模拟的表面粗糙度和空气动力阻力。参数的不确定性(例如,r残雪,r厘米,Gf,α)在极端的条件下构建Tr- - - - - -fc梯形空间也需要进一步研究。

数据可用性

地面观测数据用于支持本研究的结果公布在应用到寒冷和干旱地区科学数据中心在兰州(http://card.westgis.ac.cn/)。L1B ASTER影像是在美国地质调查局网站免费下载(https://earthexplorer.usgs.gov/)。遥感产品,包括反照率和增强型植被指数,可从相应的作者。

的利益冲突

作者宣称没有利益冲突。

确认

作者感谢黑河计划科学数据中心,中国国家自然科学基金,提供地面和近地表测量获得的黑河流域盟军遥测实验研究(HiWATER, http://westdc.westgis.ac.cn/hiwater)。这项工作得到了中国国家重点研究和发展计划(批准号2017 yfc0405801-02)的技术创新基础珠江液压研究所(批准号[2018]ky015), CRSRI研究项目(不开放。CKWV2017529 /肯塔基州),中央大学(没有基础研究基金。2017 b614x14),江苏省研究生创新研究与实践项目(没有。KYCX17_0419)。

引用

  1. Z.-L。李,r . Tang z Wan et al .,“回顾当前区域蒸散估算方法从遥感数据,”传感器,9卷,不。5,3801 - 3853年,2009页。视图:出版商的网站|谷歌学术搜索
  2. j·d·卡尔马、t·r·Mcvicar和m·f·麦凯布”估计地表蒸发:回顾方法使用遥感地表温度数据,”在地球物理调查卷,29号4 - 5,421 - 469年,2008页。视图:出版商的网站|谷歌学术搜索
  3. 王k和r·e·迪金森,“回顾全球陆地蒸散:观察、建模、气候学、和气候变化,”地球物理评论,50卷,不。2,1-54,2012页。视图:出版商的网站|谷歌学术搜索
  4. X.-G。王,黄d . w . Wang Yong, x陈,“修改SEBAL模型基于梯形地表温度和植被指数之间的关系,实际蒸散估算,“遥感》第六卷,没有。7,5909 - 5937年,2014页。视图:出版商的网站|谷歌学术搜索
  5. w·g·m·Bastiaanssen m . Menenti r . A . Feddes和A . A . m . Holtslag”土地遥感地表能量平衡算法(SEBAL)。1:制定。”《水文卷,212 - 213,198 - 212年,1998页。视图:出版商的网站|谷歌学术搜索
  6. r·g·艾伦·m·Tasumi, r . Trezza”卫星能量平衡映射蒸散与内化校准(公制)模型,”灌溉和排水工程》杂志上,卷133,不。4、395 - 406年,2007页。视图:出版商的网站|谷歌学术搜索
  7. 问:μ,f . a . Heinsch m .赵和s . w .跑步,“发展全球蒸散算法基于MODIS和全球气象数据,”环境遥感,卷111,不。4、519 - 536年,2007页。视图:出版商的网站|谷歌学术搜索
  8. 问:μm .赵,s . w .运行“MODIS全球陆地蒸散算法的改进,”环境遥感,卷115,不。8,1781 - 1800年,2011页。视图:出版商的网站|谷歌学术搜索
  9. j·m·诺曼·w·p·Kustas, k . s .休谟”方法估算土壤和植被能源通量观测的定向辐射表面温度,“农业和森林气象学,卷77,不。3 - 4、263 - 293年,1995页。视图:出版商的网站|谷歌学术搜索
  10. y杨和美国商”,混合双源方案和梯形基于框架蒸散模型(他们…)使用卫星图片:算法和模型试验,”地球物理学研究杂志:atm,卷118,不。5,2284 - 2300年,2013页。视图:出版商的网站|谷歌学术搜索
  11. j·m·桑切斯w . p . Kustas诉Caselles和m . c·安德森,“造型表面能量通量在玉米使用一块两个源模型和辐射土壤和冠层温度观测,”环境遥感,卷112,不。3、1130 - 1143年,2008页。视图:出版商的网站|谷歌学术搜索
  12. z苏”,表面能量平衡系统(实用)估计湍流热通量,”水文和地球系统科学》第六卷,没有。1,第100 - 85页,2002。视图:出版商的网站|谷歌学术搜索
  13. d . v . p·辛格,“蒸散的两个源梯形模型从卫星图像(TTME),“环境遥感,卷121,不。2、370 - 388年,2012页。视图:出版商的网站|谷歌学术搜索
  14. w·g·m·Bastiaanssen e . j . m . Noordman h . Pelgrum g .大卫b . p . Thoreson和r·g·艾伦,”SEBAL模型与遥感数据在实际现场条件下改善水资源管理,“灌溉和排水工程》杂志上,卷131,不。1,第93 - 85页,2005。视图:出版商的网站|谷歌学术搜索
  15. h·苏·m·f·麦凯布·e·f·伍德,z苏,和j·h·Prueger“蒸散在SMACEX建模:比较两种方法对当地对地区级的核预测,“水文气象学杂志》第六卷,没有。6,910 - 922年,2005页。视图:出版商的网站|谷歌学术搜索
  16. r·g·艾伦·m·Tasumi, a·莫尔斯”卫星能量平衡映射834年蒸散与内化校准(公制)应用程序,”835年灌溉和排水工程》杂志上,卷133,不。4、395 - 406年,2007页。视图:出版商的网站|谷歌学术搜索
  17. j .诺曼·m·安德森,w . Kustas”是单源遥感surface-flux模型太简单?”航会议论文集,卷852,不。170年,2006年。视图:出版商的网站|谷歌学术搜索
  18. w·p·Kustas”,估计蒸散,和两层模型的传热部分树冠覆盖,“应用气象学杂志卷,29号8,704 - 715年,1990页。视图:出版商的网站|谷歌学术搜索
  19. d . Troufleau j . p . Lhomme b . Monteny和比达尔,“显热通量和辐射表面温度在稀疏的萨赫勒地区的植被。我:知识库的实验分析−1参数”,《水文卷,188 - 189。1,第838 - 815页,1997。视图:出版商的网站|谷歌学术搜索
  20. w·j·蒂默曼w . p . Kustas m·c·安德森et al。”一个相互比对表面能量平衡算法的土地(SEBAL)和两个源能量平衡(TSEB)建模方案,“环境遥感,卷108,不。4、369 - 384年,2007页。视图:谷歌学术搜索
  21. m·安德森,j·m·诺曼·g·r·Diak et al .,“一个两个源time-integrated模型使用热红外遥感估算地表通量,”环境遥感,60卷,不。2、195 - 216年,1997页。视图:出版商的网站|谷歌学术搜索
  22. w . p . Kustas和j·m·诺曼”评价土壤和植被热通量与辐射温度预测模型使用一个简单的两个源部分树冠覆盖,“农业和森林气象学,卷94,不。1、13、1999页。视图:出版商的网站|谷歌学术搜索
  23. w . Kustas m·安德森,“热红外遥感地表建模的发展,“农业和森林气象学,卷149,不。12日,第2081 - 2071页,2009年。视图:出版商的网站|谷歌学术搜索
  24. f·李,w . p . Kustas j . h . Prueger c·m·尼尔,t·j·杰克逊,“实用程序的远程sensing-based两个源能量平衡模型在低收入和high-vegetation覆盖条件下,“水文气象学杂志》第六卷,没有。6,878 - 891年,2005页。视图:出版商的网站|谷歌学术搜索
  25. e·韦伯斯特,d .坡道,r·t·金斯福德”空间表面能量平衡算法的灵敏度气象数据在异构环境中,“环境遥感卷,187年,第319 - 294页,2016年。视图:出版商的网站|谷歌学术搜索
  26. m·德克尔·m·a . Brunke k .坂口,曾x, z . Wang和m . g . Bosilovich”评估从戈达德宇航中心再分析产品,摘要,并使用通量塔ECMWF的观察,“杂志的气候,25卷,不。6,1916 - 1944年,2012页。视图:出版商的网站|谷歌学术搜索
  27. k . Nishida r . r . Nemani s . w . et al。”一个操作遥感算法的地表蒸发,”地球物理研究杂志》的atm,卷108,不。D9, 2003年。视图:出版商的网站|谷歌学术搜索
  28. y姚明,梁,j . Yu et al .,”一个简单的温度域两个源模型估算农田地表能量通量从陆地卫星图像,”地球物理研究杂志》,122卷,2017年。视图:出版商的网站|谷歌学术搜索
  29. 毛l . s . Liu, d, l .贾”评估参数化的空气动力阻力使用实地测量传热,”水文和地球系统科学,11卷,不。2、769 - 783年,2007页。视图:出版商的网站|谷歌学术搜索
  30. g . y .邱t .矢野,k . Momii”一种改进的方法来测量蒸发从土壤裸露的基于比较干燥的土壤表面的表面温度,“《水文,卷210,不。1 - 4、93 - 105年,1998页。视图:出版商的网站|谷歌学术搜索
  31. z太阳,小王,b .松下et al .,“开发一个简单的遥感蒸散模型(Sim-ReSET):算法和模型试验,”《水文,卷376,不。3 - 4、476 - 485年,2009页。视图:出版商的网站|谷歌学术搜索
  32. h .太阳,“两级梯形:地表温度的一个新的解释和部分植被覆盖空间,”IEEE选择杂志的主题应用地球观测和遥感,9卷,不。1,第346 - 336页,2016。视图:出版商的网站|谷歌学术搜索
  33. r . Tang和z l .李“end-member-based两个源方法估算地表蒸散遥感数据,”IEEE地球科学和遥感,55卷,不。10日,5818 - 5832年,2017页。视图:出版商的网站|谷歌学术搜索
  34. d . w . Wang黄,X.-G。王,Y.-R。Liu和f .周”,估计使用梯形的土壤水分遥感地表温度和植被指数之间的关系,“水文和地球系统科学,15卷,不。5,1699 - 1712年,2011页。视图:出版商的网站|谷歌学术搜索
  35. e·m·切·r·j·哈丁”,应用聚合模型表面热通量的萨赫勒地区的老虎布什”农业和森林气象学,卷72,不。3 - 4、213 - 235年,1995页。视图:出版商的网站|谷歌学术搜索
  36. b . Choudhury s Idso, r . Reginato”分析的实证模型下的土壤热通量估算蒸发种植小麦作物的红外温度建立能量平衡方程,”农业和森林气象学,39卷,不。4、283 - 297年,1987页。视图:出版商的网站|谷歌学术搜索
  37. w·Brutsaert蒸发到大气中:理论、历史和应用程序施普林格科学与商业媒体,柏林,德国,2013年。
  38. z苏,t . Schmugge w . p . Kustas和w·j . Kustas”两个模型估计粗糙度高的地表和大气之间的传热,”应用气象学杂志,40卷,不。11日,第1951 - 1933页,2001年。视图:出版商的网站|谷歌学术搜索
  39. w·j·Massman”、“一个模型的研究kBH−1植物的表面使用拉格朗日理论“局部近场”,“《水文,卷223,不。1 - 2,27-43,1999页。视图:出版商的网站|谷歌学术搜索
  40. x曾庆红,m . Barlage r·e·迪金森y戴,g . Wang和k·奥尔森,“土地undercanopy湍流模型的治疗。”杂志的气候,18卷,不。23日,第5094 - 5086页,2005年。视图:出版商的网站|谷歌学术搜索
  41. t·卡尔森“概述”三角形法”估算表面蒸发蒸腾和土壤水分从卫星图像,”传感器,7卷,不。8,1612 - 1629年,2007页。视图:出版商的网站|谷歌学术搜索
  42. j·莱维特植物对环境胁迫的响应、学术出版社、剑桥、马、美国、1980。
  43. x, s .刘问:肖et al .,“多尺度数据集对理解复杂的eco-hydrological过程在异构绿洲系统中,“科学数据,4卷,不。1,文章ID 170083, 2017。视图:出版商的网站|谷歌学术搜索
  44. 温x、b .杨太阳x, x李,“蒸散分区通过原位氧同位素测量在绿洲农田,”农业和森林气象学卷,230 - 231,89 - 96年,2016页。视图:出版商的网站|谷歌学术搜索
  45. c . Liebethal b Huwe, t . Foken“敏感性分析两个地面热通量计算方法,”农业和森林气象学,卷132,不。3 - 4、253 - 262年,2005页。视图:出版商的网站|谷歌学术搜索
  46. t . e .线w . p . Kustas j·m·诺曼et al .,“纠正涡度相关通量在草原,低估了”农业和森林气象学,卷103,不。3、279 - 300年,2000页。视图:出版商的网站|谷歌学术搜索
  47. w·p·l .歌曲,美国Liu Kustas et al .,”应用程序的远程sensing-based两个源能量平衡模型映射领域表面通量和复合组件表面温度,“农业和森林气象学卷,230 - 231年8-19,2016页。视图:出版商的网站|谷歌学术搜索
  48. 江z, a . Huete k Didan, t .三浦”的发展没有蓝色的双波段增强型植被指数乐队”环境遥感,卷112,不。10日,3833 - 3845年,2008页。视图:出版商的网站|谷歌学术搜索
  49. d·h . Li太阳,y Yu et al .,“评估VIIRS和MODIS LST产品在中国西北干旱区,”环境遥感,卷142,不。1,第121 - 111页,2014。视图:出版商的网站|谷歌学术搜索
  50. c . y . Li黄,j .侯j .顾g .朱和李x”映射每天蒸散ASTER和MODIS图像的基于时空融合在黑河流域灌溉农业地区,中国西北,”农业和森林气象学卷,244 - 245,82 - 97年,2017页。视图:出版商的网站|谷歌学术搜索
  51. x x, y,和粉丝,“卫星表面蒸散与非参数方法:检索准确性评估在半干旱地区,”气象学的进展ID 1584316条,卷。2016年,14页,2016。视图:出版商的网站|谷歌学术搜索
  52. t, s m . Bateni s . a .马古利斯l .歌曲,s . Liu和z徐,“分区蒸散成土壤蒸发和树冠蒸腾通过两个源变分资料同化系统”水文气象学杂志,17卷,不。9日,第2370 - 2353页,2016年。视图:出版商的网站|谷歌学术搜索
  53. y的歌,l·金·g·朱,m·马”为一个简单的两个源参数估计使用贝叶斯推理蒸散模型及其应用遥感估算的潜热通量在区域范围内,“农业和森林气象学卷,230 - 231年20-32,2016页。视图:出版商的网站|谷歌学术搜索
  54. 长y, d, h .关,w .梁,c·西蒙斯和o . Batelaan”比较三种双源遥感蒸散模型在MUSOEXE-12运动:重新审视模型物理,“水资源研究,51卷,不。5,3145 - 3165年,2015页。视图:出版商的网站|谷歌学术搜索
  55. 黄c . y . Li j .顾l . Lu和李x”提高估算蒸散water-limited条件下基于实用和MODIS数据在干旱地区,”遥感,7卷,不。12日,第16814 - 16795页,2015年。视图:出版商的网站|谷歌学术搜索
  56. y, s . Liu f . Zhang et al .,“估计区域地表能量通量在异构oasis-desert表面在黑河中游HiWATER-MUSOEXE期间,“IEEE地球科学&遥感信,12卷,不。3、671 - 675年,2015页。视图:出版商的网站|谷歌学术搜索
  57. 问:壮族和b .吴”,从一种改进的两个源能量平衡模型估算蒸散使用ASTER卫星图像,”,7卷,不。12日,第6688 - 6673页,2015年。视图:出版商的网站|谷歌学术搜索
  58. l .歌曲,s . Liu x Zhang et al .,“评估和验证土壤蒸发和作物蒸腾HiWATER-MUSOEXE期间,“IEEE地球科学和遥感信,12卷,不。2、334 - 338年,2015页。视图:出版商的网站|谷歌学术搜索
  59. 苏g·f·朱,李x, y . h . et al .,”同时吸收多元数据集的两个源蒸散模型,贝叶斯方法:应用程序来春玉米在中国西北的干旱地区,”立模型开发,7卷,不。4、1467 - 1482年,2014页。视图:出版商的网站|谷歌学术搜索
  60. 井上m . s .莫兰,t·r·克拉克y,和a·维达尔,“估计作物水分亏缺使用表面空气温度和光谱植被指数之间的关系,“环境遥感卷,49号3、246 - 263年,1994页。视图:出版商的网站|谷歌学术搜索
  61. s . d .王子和s . n . Goward“全球主要生产:遥感的方法,”生物地理学杂志,22卷,不。4 - 5,815 - 835年,1995页。视图:出版商的网站|谷歌学术搜索
  62. s . j . j . Wang壮族,w . Wang Liu和z徐,“评估不确定性的涡度相关通量测量HiWATER-MUSOEXE密集的流量矩阵的基础上,“IEEE地球科学&遥感信,12卷,不。2、259 - 263年,2015页。视图:出版商的网站|谷歌学术搜索
  63. d . Baldocchi j . Finnigan k·威尔逊et al .,“测量净生态系统碳交换在高大的植被在复杂地形,”边界层气象学,卷96,不。1 - 2、257 - 291年,2000页。视图:出版商的网站|谷歌学术搜索
  64. c . Liebethal和t . Foken评价六个地面热通量参数化方法,”理论与应用气候学,卷88,不。1 - 2,43-56,2007页。视图:出版商的网站|谷歌学术搜索
  65. m·a·fiedl的用于检查电子邮件地址,“遥感数据之间的关系,表面能量平衡,植被生长区和area-averaged通量部分表面,”应用气象学杂志,35卷,不。11日,第2103 - 2091页,1996年。视图:出版商的网站|谷歌学术搜索
  66. w . p . Kustas x詹,t . j . Schmugge”结合光学和微波遥感映射在半干旱流域能量通量,”环境遥感,卷64,不。2、116 - 131年,1998页。视图:出版商的网站|谷歌学术搜索
  67. 杨k和j .王”预示校正温度的表面土壤热通量估算方法从土壤温度和湿度数据,”科学在中国系列D:地球科学,51卷,不。5,721 - 729年,2008页。视图:出版商的网站|谷歌学术搜索
  68. w . p . Kustas和j·m·诺曼”两个源的方法估算湍流通量使用多个角度热红外观测,”水资源研究,33卷,不。6,1495 - 1508年,1997页。视图:出版商的网站|谷歌学术搜索
  69. j.p. Lhomm、b . Montenya和m . Amadoua”估计显热通量辐射温度在稀疏的小米,“农业和森林气象学,卷68,不。1 - 2、77 - 91年,1994页。视图:出版商的网站|谷歌学术搜索
  70. j . r . Zhang田、h·苏x太阳,s . Chen和j .夏”两个操作两层模型的改进地面热通量检索”传感器,8卷,不。10日,6165 - 6187年,2008页。视图:出版商的网站|谷歌学术搜索
  71. x詹、w . p . Kustas和k . s .休谟”一个相互比较研究模型的显热通量部分树冠表面与遥感地表温度,“环境遥感,卷。58岁的没有。3、242 - 256年,1996页。视图:出版商的网站|谷歌学术搜索
  72. y, t·杜·李,李,r·丁和l .通“量化的玉米水吸收来自不同层次和根区交替沟灌下使用稳定的氧同位素,”农业水管理卷。168年,35-44,2016页。视图:出版商的网站|谷歌学术搜索
  73. j·m·诺曼·f·贝克尔,“术语在热红外遥感的自然表面,”遥感的评论,12卷,不。3 - 4、153 - 166年,1995页。视图:出版商的网站|谷歌学术搜索
  74. p·d·Colaizzi w . p . Kustas m·c·安德森et al。”两个源能量平衡模型的估算蒸散使用组件和复合表面温度,“水资源的进步,50卷,不。6,134 - 151年,2012页。视图:出版商的网站|谷歌学术搜索
  75. f, g .秋y, y, y Xiong,“蒸散估算及其分区的基础上延长三温度模型和MODIS产品,”《水文,卷498,不。1,第220 - 210页,2013。视图:出版商的网站|谷歌学术搜索
  76. a . Chehbouni y Nouvellon, j。Lhomme et al .,“估计表面显热通量使用辐射表面温度的双重角度观察,“农业和森林气象学,卷108,不。1,55 - 65、2001页。视图:出版商的网站|谷歌学术搜索
  77. k . s .休谟w . p . Kustas m . s . Moran w·d·尼科尔斯和m . a . Weltz”发射率和表面温度的变化对表面植被稀疏,”水资源研究,30卷,不。5,1299 - 1310年,1994页。视图:出版商的网站|谷歌学术搜索
  78. w . p . Kustas j·m·诺曼·t·j·Schmugge et al .,“映射与辐射温度、地表能量通量”热遥感地表过程CRC出版社,页205 - 253年,纽约,纽约,美国,2004年。视图:出版商的网站|谷歌学术搜索
  79. r·艾伦河,r . Trezza j·m·h·亨德里克斯,w . Bastiaanssen和j·克亚斯高,“卫星等估计农业利用SEBAL和度规,“水文过程,25卷,不。26日,第4027 - 4011页,2011年。视图:出版商的网站|谷歌学术搜索
  80. y, h·苏·r·张,j .田l·李,“一个增强两个源土地(ETEML):蒸散模型算法和评价,“环境遥感,卷168,不。168年,54 - 65年,2015页。视图:出版商的网站|谷歌学术搜索
  81. m·c·安德森,j·m·诺曼·w·p·Kustas et al .,“植被影响凝结在两个源模型的估计地表能量通量从农业景观SMACEX期间,“水文气象学杂志》第六卷,没有。6,892年,页2005。视图:出版商的网站|谷歌学术搜索
  82. h .关和j·l·威尔逊,“混合双源模型潜在蒸发和植物蒸腾作用的分区,”《水文,卷377,不。3 - 4、405 - 416年,2009页。视图:出版商的网站|谷歌学术搜索
  83. r·d·杰克逊,s . b . Idso r . j . Reginato和p . j .品特“作为作物冠层温度水压力指标,”水资源研究,17卷,不。4、1133 - 1138年,1981页。视图:出版商的网站|谷歌学术搜索
  84. a·维达尔和a .毕雷矿泉水”技术报告分析的简化关系估计每天蒸散从卫星热红外数据,”国际遥感杂志》上,10卷,不。8,1327 - 1337年,1989页。视图:出版商的网站|谷歌学术搜索

版权©2020王小刚et al。这是一个开放分布式下文章知识共享归属许可,它允许无限制的使用、分配和复制在任何媒介,提供最初的工作是正确引用。

74年 的观点 | 171年 下载 | 0 引用
PDF 下载引用 引用
下载其他格式更多的
订单打印副本订单

我们致力于分享发现相关COVID-19尽可能快速和安全。任何作者提交COVID-19纸应该通知我们help@hindawi.com以确保他们的研究顺利和尽快预印本服务器上可用。我们将提供无限的豁免的出版费用接受COVID-19相关文章。注册在这里作为一个评论家,帮助快速新提交。