NIM陆面过程模式的研究Ⅱ：青藏高原夏季陆面过程的数值模拟

采用南京气象学院（NIM）5层陆面过程模式，利用1979年5－8月“青藏高原气象科学实验”资料模拟和分析了夏季青藏高原不同地区的陆面特征和地表能量特征。并将模拟值与根据观测资料计算得到的感热和潜热以及观测得到的净辐射、土壤温度、土壤热通量进行了对比。结果表明，NIM5层陆面过程模式可以模拟青藏高原夏季不同下垫面情形下的能量交换过程。     

典型干旱区陆面模式模拟检验

利用2000年9月至2001年8月"我国西北干旱区陆-气相互作用试验(NWC-ALIEX)"敦煌站的陆面过程观测资料,基于已有的参数化结果,模拟了敦煌主要陆面特征。结果表明:典型干旱区敦煌夏季感热通量与潜热通量差异显著,感热几乎是潜热的4倍,冬季二者都很小。模式对地表温度模拟较好,但高估了浅层土壤湿度峰值;对向上辐射模拟较好,但对净辐射的模拟存在较大偏差;同时高估了地表能量的峰值。模式结果表明参数化方案对干旱区陆面过程模式具有一定的改进作用。     

CoLM模式对青藏高原中部BJ站陆面过程的数值模拟

利用公共陆面模式Common Land Model(CoLM)及"全球协调加强观测计划之亚澳季风青藏高原试验"(CAMP/Tibet)中那曲地区Bujiao(BJ)站2002—2004年的观测资料对该地区进行了单点数值模拟试验。通过比较模拟与观测的地表能量通量,表明CoLM较成功地模拟了该地区的能量分配。模式对向上的短波辐射、向上的长波辐射、净辐射及土壤热通量模拟得较好,但冬季存在偏差。进一步比较了模拟和观测的土壤温度及土壤湿度,发现浅层60 cm土壤温度模拟较好,深层存在偏差,表现为土壤温度变化滞后于实际变化。土壤湿度总体偏小,尤其是冬季冻结期,土壤冻融过程中忽略了土壤液态水在温度0℃以下仍能存在,含冰量模拟偏高。     

SiB2和SiB3对高寒草甸和茶树地表能量通量模拟的比较

运用简单生物圈模式第2版(SiB2)和第3版(SiB3),分别模拟青藏高原两个观测站(那曲、安多)和长江三角洲苏州东山观测站的近地面能量收支,并与相应观测数据进行比较研究,分析SiB2、SiB3模拟结果和观测资料产生差异的原因,以此来认识上述地区地表能量收支特点。结果表明,SiB2和SiB3模拟的近地面能量通量与观测数据有较好的一致性。对感热通量,那曲和安多站SiB3比SiB2模拟的结果更接近观测资料,但苏州站SiB2模拟的结果与观测资料更吻合;对潜热通量,SiB3比SiB2模拟的日变化与观测资料更一致,SiB3的模拟结果与观测资料(除苏州站外)相关系数都在0.8以上;对地表土壤热通量,SiB2和SiB3模拟结果与观测数据相关系数都在0.8以上;对净辐射通量,SiB2和SiB3模拟结果与观测资料相关系数接近1.0。与SiB2相比,SiB3引用通用陆面模式的土壤描述并增加对冠层空间层温度、湿度和痕量气体的预报,使其能够改善潜热通量和土壤热通量的模拟,但对复杂下垫面的感热和净辐射通量模拟能力提高不明显。     

陆面过程模式CoLM和NCAR_CLM3.0对中国典型森林生态系统陆气相互作用的模拟 I. 不同模式模拟结果的初步分析

CoLM(Common Land Model)和NCAR_CLM3.0(NCAR Community Lanol Model 3.0)是目前国际上广为应用的两个发展比较完善的陆面过程模式.本研究利用中国陆地生态系统通量观测研究网络(Chi-naFLUX)在长白山温带混交林和千烟洲亚热带人工针叶林观测站点的长期连续强化观测资料,对这两个模式在上述地区的模拟性能进行了初步评估.同观测资料的对比表明,两个模式均能较好地模拟出观测站点地表能量和水分平衡的基本特征,其中,CoLM对潜热通量的模拟性能更好.以对潜热通量为期1年的日均值的模拟为例,COLM和CLM3.0在长白山观测站模拟值和观测值时间序列的相关系数分别为0.80和0.65,在千烟洲站分别为0.69和0.64,均通过了0.01的信度检验;两个模式对全年平均的模拟与观测日平均值的比值在长白山分别为1.21和0.86,在千烟洲分别为0.83和0.60.研究结果表明,这两个陆面过程模式可以作为研究这两种典型森林生态系统陆气交换的基本工具.同时,对模式模拟性能差异的深入分析将有助于进一步改进陆面模式的参数化过程,为相关研究奠定更坚实的基础.     

青藏高原西部冻融期陆面过程的模拟分析

利用CAMP/Tibet中CEOP-EOP3改则站2002年10月—2003年9月的观测资料作为强迫场,运用陆面过程模式CoLM(Common Land Model),对青藏高原西部陆面特征的模拟研究表明,在高原西部地表能量平衡过程中,冬半年,感热通量占主要地位,潜热通量较小;尤其在冻结期,潜热通量几乎等于零。但在高原西部的融冻期,潜热通量有显著变化。在干季向湿季转化时段的5月中下旬,表层土壤由于融冻而引起的频繁水分相变,使得潜热通量随之变化并开始增加,Bowen比由大变小。地表有效通量的变化与降水及土壤表层频繁的冻结—消融相联系。     

改则地区气候变化对植被生理过程的影响及反馈效应的模拟研究

利用1997年10月1日至1998年9月30日设置在青藏高原西部改则的自动气象站观测资料作为强迫场,采用大气-植被相互作用模式(AVIM)对改则地区气候变化对植被生长过程的影响及反馈效应进行了模拟研究。结果表明,AVIM模式对青藏高原西部陆面过程具有一定模拟能力,能够较真实地模拟出地表特征量的变化特点。通过敏感性试验发现,青藏高原气候变化对植被生理生长过程有明显影响:降水增加有利于植被生长,尤其在雨季最为明显,其他季节无太大变化;气候变暖对植被生理过程的综合作用是植被净光合作用的变化,即春季增强,夏季减弱,秋季和冬季变化不大;"暖湿化"对高原植被生态系统的影响主要是春季和夏季植被活动增强,尤其春季最为明显。植被物理特性参数可以在相当大程度上改变陆面过程,进而导致高原热源发生变化,因此,为准确估计地表能量收支,对模式陆面参数进行深入研究是必要的。     

青藏高原单点地气交换过程的模拟试验

利用GAME-Tibet 1998 IOPs的外场观测资料,采用单点陆面过程气柱模式CLSM对青藏高原的昌都(Amdo)站和改则(Gaize)站的地气交换过程进行了模拟试验,模拟结果与观测分析的对比表明：模式对两个站陆面特征的变化和陆气之间的交换给出了比较合理的模拟结果,模拟结果与观测结果较一致。这在一定程度上表明,CLSM模式具有描述夏季高原地区两种不同下垫面陆面状况基本过程的能力。     

藏北高原地表能量和边界层结构的数值模拟

利用耦合了NCAR LSM陆面过程的中尺度模式MM5V3.7和2002年8月CAMP/Tibet加强期的观测资料,对藏北高原地区地气交换过程进行了48 h模拟研究。模式较好地模拟了该地区的山谷风环流;并将模拟的地表通量在中尺度区域上与NCEP/NCAR全球大气再分析格点资料(NNRP)获得的结果进行了比较,同时也与单站的实测值进行了比较,结果显示:模拟的地表通量与NNRP得到的结果比较吻合,同时可以得到雨季时藏北、藏东地区潜热通量大于感热通量,而高原西部感热通量大于潜热通量,这与观测试验分析结果一致;与单站试验结果比较,模拟的感热通量与实测值一致,潜热通量的模拟值和实测值有一定差别。模拟的边界层位温廓线与实测值比较,模式模拟的对流混合层和夜间残留层都与实测结果吻合,但模拟的混合层高度较实测值高。由此来看,中尺度模式MM5V3.7能够较好地模拟藏北高原的地表能量和边界层结构特征,但还需要进一步完善陆面过程和物理过程参数化方案。     

青藏高原积雪对陆面过程热量输送的影响研究

利用1979-2016年中国区域长时间序列逐日雪深资料,分析了青藏高原积雪深度与积雪日数的分布及变化特征,并将积雪期划分为三个阶段(积累期、鼎盛期和消融期),结合ERA-Interim月平均再分析资料,分析了积雪与地表热状况(气温、地表和土壤温度)和能量输送量(地表净短波辐射、地表净长波辐射、感热通量、潜热通量、地表热通量和土壤热通量)的相关关系,初步探讨了积雪在高原陆面过程中的作用。结果表明:研究时间范围内青藏高原积雪(深度和日数)主要呈减少趋势,仅在黄河源区及高原边缘地区为增加趋势,积雪鼎盛阶段(1-2月)的减少趋势最显著;高原积雪对地表主要起降温作用,深层土壤温度对积雪的响应存在滞后性,积雪的减少抑制了土壤向上的热量输送进而不利于冻土的发育;高原积雪与地表感热和地表热通量主要呈现负相关关系,潜热通量与积雪也呈负相关特征但比感热通量的相关性小。由于ERA-Interim资料对高原积雪深度的描述与本研究使用的卫星遥感积雪深度存在较大偏差(包括空间分布、气候倾向率、年际变化以及绝对大小等),导致本研究中积雪与地表热状况和热通量的相关度不高,需要通过陆面模式模拟做进一步探讨。     

The solubility and behaviour of acid gases in the marine aerosol

The following Henry's law constants (K _H/mol²kg^-2atm^-1) for HNO₃ and the hydrohalic acids have been evaluated from available partial pressure and other thermodynamic data from 0°–40°C, 1 atm total pressure: HNO ₃ , 40°C–5.85×10⁵; 30°C–1.50×10⁶; 25°C–2.45×10⁶; 20°C–4.04×10⁶; 10°C–1.15×10⁷; 0°C–3.41×10⁷. HF, 40°C–3.2; 30°C–6.6; 25°C–9.61; 20°C–14.0; 10°C–32.0; 0°C–76. HCl, 40°C–4.66×10⁵; 30°C–1.23×10⁶; 25°C–2.04×10⁶; 20°C–3.37×10⁶; 10°C–9.71×10⁶; 0°C–2.95×10⁷. HBr, 40°C–2.5×10⁸; 30°C–7.5×10⁸; 25°C–1.32×10⁹; 20°C–2.37×10⁹; 10°C–8.10×10⁹; 0°C–3.0×10¹⁰. HI, 40°C–5.2×10⁸; 30°C–1.5×10⁹; 25°C–2.5×10⁹; 20°C–4.5×10⁹; 10°C–1.5×10¹⁰; 0°C–5.0×10¹⁰. Simple equilibrium models suggest that HNO₃, CH₃SO₃H and other acids up to 10x less soluble than HCl displace it from marine seasalt aerosols. HF is displaced preferentially to HCl by dissolved acidity at all relative humidities greater than about 80%, and should be entirely depleted in aged marine aerosols.     

Rate constants for the reactions of the NO3 radical with HCOOH/HCOO− and CH3COOH/CH3COO− in aqueous solution between 278 and 328 K

The kinetics of the aqueous phase reactions of NO₃ radicals with HCOOH/HCOO^– and CH₃COOH/CH₃COO^– have been investigated using a laser photolysis/long-path laser absorption technique. NO₃ was produced via excimer laser photolysis of peroxodisulfate anions (S₂O ₈ ^2– ) at 351 nm followed by the reactions of sulfate radicals (SO ₄ ^– ) with excess nitrate. The time-resolved detection of NO₃ was achieved by long-path laser absorption at 632.8 nm. For the reactions of NO₃ with formic acid (1) and formate (2) rate coefficients ofk ₁=(3.3±1.0)×10⁵ l mol^–1 s^–1 andk ₂=(5.0±0.4)×10⁷ l mol^–1 s^–1 were found atT=298 K andI=0.19 mol/l. The following Arrhenius expressions were derived:k ₁(T)=(3.4±0.3)×10¹⁰ exp[–(3400±600)/T] l mol^–1 s^–1 andk ₂(T)=(8.2±0.8)×10¹⁰ exp[–(2200±700)/T] l mol^–1 s^–1. The rate coefficients for the reactions of NO₃ with acetic acid (3) and acetate (4) atT=298 K andI=0.19 mol/l were determined as:k ₃=(1.3±0.3)×10⁴ l mol^–1 s^–1 andk ₄=(2.3±0.4)×10⁶ l mol^–1 s^–1. The temperature dependences for these reactions are described by:k ₃(T)=(4.9±0.5)×10⁹ exp[–(3800±700)/T] l mol^–1 s^–1 andk ₄(T)=(1.0±0.2)×10¹² exp[–(3800±1200)/T] l mol^–1 s^–1. The differences in reactivity of the anions HCOO^– and CH₃COO^– compared to their corresponding acids HCOOH and CH₃COOH are explained by the higher reactivity of NO₃ in charge transfer processes compared to H atom abstraction. From a comparison of NO₃ reactions with various droplets constituents it is concluded that the reaction of NO₃ with HCOO^– may present a dominant loss reaction of NO₃ in atmospheric droplets.     

Characteristics of ionic components in precipitation in Kitakyushu City,Japan

Precipitation samples were collected by filtrating bulk sampler in Kitakyushu City, Japan, from January 1988 to December 1990. Volume weighted annual mean of pH was 4.93, but the pH distribution indicated that most probable value lay in the range pH 6.0–6.4. Volume weighted annual mean concentrations of major ionic components were as follows; SO ₄ ^2– : 84.2, NO ₃ ^– : 28.1, Cl^–: 86.3, NH ₄ ⁺ : 45.5, Ca²⁺: 63.3, Mg²⁺: 27.0, K⁺: 3.4, Na⁺: 69.0 µ eq l^–1. The highest concentrations of these ionic components were observed in winter and the lowest occurred in the rainy season. The ratio of ex-SO ₄ ^2– /NO ₃ ^– exhibited the lowest ratio in summer, and the highest ratio in winter. Good correlations were obtained between Cl^– and Na⁺, ex-SO ₄ ²⁺ and ex-Ca²⁺, NO ₃ ^– and ex-Ca²⁺, and NH ₄ ⁺ and ex-SO ₄ ^2– , respectively. However, no correlation between Cl^– and Na⁺ with Ca²⁺ was observed. The relationship of H⁺ with (ex-SO ₄ ^2– + NO ₃ ^– ) - (ex-Ca²⁺ + NH ₄ ⁺ ) indicated positive correlation.     

A K-Theory of Dispersion,Settling and Deposition in the Atmospheric Surface Layer

Dispersion estimates with a Gaussian plume model are often incorrect because of particle settling (β), deposition (γ) or the vertical gradient in diffusivity (K _v (z) = K ₀ + μz). These “non-Gaussian” effects, and the interaction between them, can be evaluated with a new Hankel/Fourier method. Due to the deepening of the plume downwind and reduced vertical concentration gradients, these effects become more important at greater distance from the source. They dominate when distance from the source exceeds L _β = K ₀ U/β ², L _γ  = K ₀ U/γ ² and L _μ = K ₀ U/μ ² respectively. In this case, the ratio β/μ plays a central role and when β/μ = 1/2 the effects of settling and K gradient exactly cancel. A general computational method and several specific closed form solutions are given, including a new dispersion relation for the case when all three non-Gaussian effects are strong. A more general result is that surface concentration scales as C(x) ~ γ ⁻² whenever deposition is strong. Categorization of dispersion problems using β/μ, L _γ and L _μ is proposed.     

Nitrogen and sulfur species in Antarctic aerosols at Mawson,Palmer Station,and Marsh (King George Island)

High volume bulk aerosol samples were collected continuously at three Antarctic sites: Mawson (67.60° S, 62.50° E) from 20 February 1987 to 6 January 1992; Palmer Station (64.77° S, 64.06° W) from 3 April 1990 to 15 June 1991; and Marsh (62.18° S, 58.30° W) from 28 March 1990, to 1 May 1991. All samples were analyzed for Na⁺, SO ₄ ^2– , NO ₃ ^– , methanesulfonate (MSA), NH ₄ ⁺ ,²¹⁰Pb, and⁷Be. At Mawson for which we have a multiple year data set, the annual mean concentration of each species sometimes vary significantly from one year to the next: Na⁺, 68–151 ng m^–3; NO ₃ ^– , 25–30 ng m^–3; nss SO ₄ ^2– , 81–97 ng m^–3; MSA, 19–28 ng m^–3; NH ₄ ⁺ , 16–21 ng m^–3;²¹⁰Pb, 0.75–0.86 fCi m^–3. Results from multiple variable regression of non-sea-salt (nss) SO ₄ ^2– with MSA and NO ₃ ^– as the independent variables indicates that, at Mawson, the nss SO ₄ ^2– /MSA ratio resulting from the oxidation of dimethylsulfide (DMS) is 2.80±0.13, about 13% lower than our earlier estimate (3.22) that was based on 2.5 years of data. A similar analysis indicates that the ratio at Palmer is about 40% lower, 1.71±0.10, and more comparable to previous results over the southern oceans. These results when combined with previously published data suggest that the differences in the ratio may reflect a more rapid loss of MSA relative to nss SO ₄ ^2– during transport over Antarctica from the oceanic source region. The mean²¹⁰Pb concentrations at Palmer and Marsh and the mean NO ₃ ^– concentration at Palmer are about a factor of two lower than those at Mawson. The²¹⁰Pb distributions are consistent with a²¹⁰Pb minimum in the marine boundary layer in the region of 40°–60° S. These features and the similar seasonalities of NO ₃ ^– and²¹⁰Pb at Mawson support the conclusion that the primary source regions for NO ₃ ^– are continental. In contrast, the mean concentrations of MSA, nss SO ₄ ^2– , and NH ₄ ⁺ at Palmer are all higher than those at Mawson: MSA by a factor of 2; nss SO ₄ ^2– by 10%; and NH ₄ ⁺ by more than 50%. However, the factor differences exhibit substantial seasonal variability; the largest differences generally occur during the austral summer when the concentrations of most of the species are highest. NH ₄ ⁺ /(nss SO ₄ ^2– +MSA) equivalent ratios indicate that NH₃ neutralizes about 60% of the sulfur acids during December at both Mawson and Palmer, but only about 30% at Mawson during February and March.     

Further studies of the temperature dependence of the rate coefficients for the reactions of OH with a series of ethers under simulated atmospheric conditions

The following temperature-dependent rate coefficients (k/cm³ molecule^–1 s^–1) of the reactions of hydroxyl radicals with aliphatic ethers have been determined over the temperature range 247–373 K by a competitive flow technique: diethyl ether,k _OH=5.2×10^–12 exp[(262±150)/T]; methyln-butyl ether,k _OH=5.4×10^–12 exp[(309±150)/T]; ethyln-butyl ether,k _OH=7.3×10^–12 exp[(335±150)/T]; di-n-butyl ether,k _OH=5.5×10^–12 exp[(502±150)/T] and di-n-pentyl ether,k _OH=8.5×10^–12 exp[(417±150)/T]. The data have been measured relative to the rate coefficientk(OH + 2,3-dimethylbutane)=6.2×10^–12 cm³ molecule^–1 s^–1 independent of temperature.Previous discrepancies in the room-temperature rate coefficients for the OH reactions with ethyln-butyl ether and di-n-butyl ether, obtained in the flow and static experiments of Bennett and Kerr (J. Atmos. Chem. 8, 87–94, 1989;10, 29–38, 1990) compared with those of Wallingtonet al. (Int. J. Chem. Kinet. 20, 541–547, 1988;21, 993–1001, 1989) and of Nelsonet al. (Int. J. Chem. Kinet. 22, 1111–1126, 1990) have been resolved. The results are considered in relation to the available literature data and evaluated rate expressions are deduced where possible. The data are also discussed in terms of structure-activity relationships.     

Kinetics of the reactions of hydroxyl radicals with aliphatic ethers studied under simulated atmospheric conditions: Temperature dependences of the rate coefficients

Rate coefficients have been measured for the reactions of hydroxyl radicals with five aliphatic ethers over the temperature range 242–328 K. Competitive studies were carried out in an atmospheric flow reactor in which the hydroxyl radicals were generated by the photolysis of methyl nitrite in the presence of air containing nitric oxide. The reaction of OH with 2,3-dimethyl-butane was used as the reference reaction and the following Arrhenius parameters have been obtained for the reactions: OH+RORproducts: