主要计算通量方法式子汇总.docx

ChineseEcosystemResearchNetwork

海洋、陆地生态系统

碳循环和通量观测研究方法

一、 涡度相关方法(EddyCorrelationMethod)

分别有下列通量的涡度:

感热通量：H=p?Cp?①0'

潜热通量：E=p?IV-^q

CO2热通量:

E=p?wc’

w=1jwdt乂1ILw

TOC\o"1-5"\h\z0 t=1

0=1jodtr1尤0

o t=1

q=1jqdtr1艺q

0 t=1

c=1jcdtr1￡c

0 t=1

二、波文比（Bowen’sRatio）方法

（1）

A=H+E=R-G-S

n

（2）

（3

（3）

（4）

基于地表能量平衡，有效能量。有

E=1AH=PA

p+1 P+1

o CA0

p=H/E=—XAq

由上式利用剖面资料，即可计算出波文比B,然后利用方程(3)计算出H和E。

空气动力学方法

空气动力学方法是根据近地面层空气动力学特性，计算能量和物质通量的输送过程。风速、温度、湿度、co2或氧化亚氮输送梯度表达式为

(5)

(6)⑺(8)TOC\o"1-5"\h\zd$_ -H

(6)

⑺

(8)

dzpCku(z一d)"h

dq_ -LE

dzpLku(z-d)3可

v*

d -F

N2O,CO2=N2O,CO2、dT^ku*(z-d)

由(5)—(8)可得

t=

t=pk2(z一d)2

(du￥

du

dudB~

H=-pCk2(z一d)2__(<p?J-i

(10)

dudq,

dudq, 、

LE=-pLk2(z一d)2__(^q)-i

― dudC ,

F=—pk2(z—d)2 %O,CQ(qq

n2O，CO2 dz dz

(11)

)-1

m N2O，CO2

(12)

式中，k为Karman常数；y为湿度表常数，Y=*=0.67hapC—i；d为位

移长度(d=，h为植被高)；q,q,q,q 分别为风速、温度、

mhwN2O，CO2

湿度和二氧化碳及氧化亚氮的稳定度通用函数，

TOC\o"1-5"\h\zZ r

中J)=1+。广当—>0 (13)

mL mL L

zx Z、 r

中(Z)=(1+Y—)4,当三<0 (14)

mL mL—l

Zz一

9(z)=9(z)=9 (z)=1+。,当—>0 (15)

hLwL N2O,CO2L hL L ()

7 Z一

9(z)=9(z)=9 (z)=(1+Y_)-L当了<0 (16)

hLwL N2O,CO2L hL L

系数Pm，5,Yh见下表。

表风、温、湿稳定度函数表达式系数

来源

Pm

Ym

Ph

Yh

k

Businger(1971)

Panison(1970)

/

Weeb(1970)

Dyer&hicks(1970)

/

/

f状<0]二,R>0

t1-5Ri)

I

式中，；L—Moni-Obukhov长度;

R.—Richardson数。有

调

方 pC0"u3

R=g°Z;L=L七

i0f"2 kgH

3z/

热量平衡方法

这是一种以能量守恒定律为基础的计算方法。实际上是一种余项法。在环境生态研究中如有辐射观察资料时，使用此法较好。地表面热量平衡方程为

Rn=H+LE+G

或 Rn=-pCpK兽-pLK藉+G （）

其中，Rn为辐射平衡，其余符号都是已知的。如果以差分代替微分，并从中解出K,可得

v （R-G）A

K=

pCAT+pLAq

8厂 而 △厂 Aq

将~8_和~c_与成差分形式~A_和~A_，则

（R—G）=-pCK半-pLKAL

或直接求出湍流热通量和蒸发耗热项，则有

dTH=-pCpK__=pCpR-Gn1+匕华C

dT

H=-pCpK__=pCp

R-G

n

1+匕华

CAT

p

R-G

n

n

LE=-pLK四=pLd_vpCAT+pLAqAz CAT

LE=-pLK四=pLd_

vpCAT+pLAqAz CAT+1

LAq

蒸发散量计算的微气象法

1.涡动相关法

地球表面能量、物质的输送是由大大小小的极不规则的湍流涡旋完成的。通常把这种湍流运动分成两部分：即平