短时段水面蒸发量计算方法的选择

分别采用折算数法，彭曼公式法，道尔顿公式法计算都昌蒸发站旬与日短时段水面蒸发量，比较表明，用道尔顿公式计算短时段水面蒸发量的误差最小，彭曼公式次之，折算系数法的误差最大，不宜作短时段水面蒸发量的计算。     

冬小麦农田日蒸散量的计算

本文从小气候观测资料着手，采用彭曼法、能量平衡法、波温比法和空气动力学等方法，对处于抽穗至乳熟期的冬小麦农田日蒸散量做了尝试性计算。着重考虑了彭曼公式的修正，并以水量平衡法为标准，对以上各方法的精度做了评价与误差分析。结果表明，订正后的彭曼公式可较为准确地计算各种能量、水分供应条件下有作物覆盖农田的日蒸散量，其它方法则存在较明显的不确定性误差。     

彭曼公式应用中的两个问题的探讨

通过对水面蒸发量的模拟比较,揭示了应用彭曼公式计算水面蒸发量中存在的两个问题:E_n项的权重偏小和E_a计算方式不当。提出了解决办法,建立了经改进后的修正彭曼公式,其精度有明显提高。     

我国蒸发力计算的气候学方法

目前,测定蒸发力还存在许多经济和技术上的困难,所以现在多采用气候学方法计算。彭曼综合法是其中最好的方法之一。但因彭曼公式中系数的确定是经验性的,具有局地性,直接应用于我国蒸发力的计算误差较大。本文利用我国10个蒸发实验站水面蒸发实测资料作为检验标准,逐项修正彭曼公式中各物理力量的系数,建立了适合我国实     

我国蒸发力计算的气候学方法

本文根据我国的具体情况,利用10个大型蒸发池自由水面蒸发的实测资料,对计算水面蒸发力的彭曼公式作了修订。修正后的彭曼公式可用于我国水面蒸发力的计算;乘以合适的折算系数和作物系数,还可用来估算我国的自然蒸发力和农田可能蒸散。     

蒸发力的计算及其动态变化特征分析

利用彭曼公式对广西忻城站蒸发力进行了气候学计算,并对其50a动态变化趋势进行了年和生长季等不同时段的分布规律分析.     

计算可能蒸发的一个气候学实用方法

在热量平衡基础上，综合彭曼公式和布德科公式的计算理论，推导出一个计算可能蒸发的气候学实用方法。在郑州地区的试验计算取得了令人满意的满意。     

三种计算可能蒸发方法的比较

彭曼(Penman)、桑斯维特(Thor-nthwaite)和布德科(Будыко)计算可能蒸发的方法,是当前国内外普遍采用的三种方法。本文根据这三种计算方法的公式,对它们的精确度、计算复杂性等方面得出了比较结果。     

麦田蒸散量的计算方法

作物叶茎蒸腾与株间土壤蒸发之和叫蒸散。作物各生育期的蒸散量直接反映了作物的需水情况,对于研究农田水分平衡、确定旱涝指标以及农业气候区划等都有重要意义。 五十年代初,彭曼根据能量平衡方程,导出一个蒸散量计算公式,但由于所涉及的物理量多而且难于精确测定,故其应用受到限制。布雷特莱—泰勒尔在彭曼公式的基础上,根据统计分析,略去了土壤热通量等对蒸散影响很小的项,建立了如下的计算模式;     

蒸发量的计算及其动态变化特征分析

应用彭曼公式对蒸发量进行了气候学计算 ,并对近 40a灯塔站蒸发量的动态变化趋势进行了年和生长季等不同时段的分布规律分析 ,以便为农业水分分析、农业气候影响评价、作物水分平衡等提供气候依据。     

The variation of the heat sources in East China in the early summer of 1984 and their effects on the large-scale circulation in East Asia

The distributions and daily variations of the apparent heat source (Q1) and the apparent moisture sink (Q2) in East China in the early summer of 1984 have been estimated with the budget calculation method. It has been found that during this time period, there occurred three significant episodes of strong heating that corresponded to the three events of heavy rainfalls prior to, during and post to the onset of mei-yu (plum rains). The peaks of Q1 were generally found at 200 hPa, with the heating rate of 6°－10°C/day observed, while the peaks of Q2 were located at about 700 hPa, with their magnitudes being 12o-20°C/day. The vertical distribution of Q1 and Q2 indicates the importance of eddy vertical flux. In other words, the convective activity plays a very important role in the processes of precipitation in East Asia in the early summer. This result is different from the finding obtained by Luo and Yanai (1984) in their calculation of the case of 1979. They pointed out that in the early summer of 1979 the continuous precipitation dominated the region of East China.Among the three terms of Q1 and Q2, the maximum contribution was made from the adiabatic term, which was caused by strong ascending motion. The adiabatic cooling produced by this term may compensate for the heating created by the condensation process.In addition, it has been revealed that the three significant heating processes were closely related to the seasonal transition from spring to summer in East China. One major synoptic event associated with it showed up in the sudden jump of the upper tropospheric, subtropical jet-stream from 30°N to 40°N. So did the plane-tary frontal zone in East China.     

Spectral Characteristics and Correction of Long-Term Eddy-Covariance Measurements Over Two Mixed Hardwood Forests in Non-Flat Terrain

We present turbulence spectra and cospectra derived from long-term eddy-covariancemeasurements (nearly 40,000 hourly data over three to four years) and the transferfunctions of closed-path infrared gas analyzers over two mixed hardwood forests inthe mid-western U.S.A. The measurement heights ranged from 1.3 to 2.1 times themean tree height, and peak vegetation area index (VAI) was 3.5 to 4.7; the topographyat both sites deviates from ideal flat terrain. The analysis follows the approach ofKaimal et al. (Quart. J. Roy. Meteorol. Soc. 98, 563–589, 1972) whose results were based upon 15 hours of measurements atthree heights in the Kansas experiment over flatter and smoother terrain. Both thespectral and cospectral constants and stability functions for normalizing and collapsingspectra and cospectra in the inertial subrange were found to be different from those ofKaimal et al. In unstable conditions, we found that an appropriate stabilityfunction for the non-dimensional dissipation of turbulent kinetic energy is of the form ₍₎ = (1 - b^-)^-1/4 - c^-, where representsthe non-dimensional stability parameter. In stable conditions, a non-linear functionG_xy() = 1 + b_xy^c _xy (c_xy < 1) was found to benecessary to collapse cospectra in the inertial subrange. The empirical cospectralmodels of Kaimal et al. were modified to fit the somewhat more (neutraland unstable) or less (stable) sharply peaked scalar cospectra observed over forestsusing the appropriate cospectral constants and non-linear stability functions. Theempirical coefficients in the stability functions and in the cospectral models varywith measurement height and seasonal changes in VAI. The seasonal differencesare generally larger at the Morgan Monroe State Forest site (greater peak VAI) andcloser to the canopy.The characteristics of transfer functions of the closed-path infrared gas analysersthrough long-tubes for CO₂ and water vapour fluxes were studied empirically. This was done by fitting the ratio between normalized cospectra of CO₂ or watervapour fluxes and those of sensible heat to the transfer function of a first-order sensor.The characteristic time constant for CO₂ is much smaller than that for water vapour. The time constant for water vapour increases greatly with aging tubes. Three methods were used to estimate the flux attenuations and corrections; from June through August, the attenuations of CO₂ fluxes are about 3–4% during the daytime and 6–10% at night on average. For the daytime latent heat flux (Q_E), the attenuations are foundto vary from less than 10% for newer tubes to over 20% for aged tubes. Correctionsto Q_E led to increases in the ratio (Q_H + Q_E)/(Q^* - Q_G) by about 0.05 to0.19 (Q_H is sensible heat flux, Q^* is net radiation and Q_G is soil heat flux),and thus are expected to have an important impact on the assessment of energy balanceclosure.     

Net air�Csea surface heat flux during 1984�C2004 over the North Pacific and North Atlantic oceans (10��N�C50��N): annual mean climatology and trend

Using the Objectively Analyzed air?Csea Fluxes dataset (and also the National Oceanography Centre Southampton Flux Dataset v2.0), we examined both the annual mean climatology and trend of net air?Csea surface heat flux (Q _net) for 1984?C2004 over the North Pacific and North Atlantic oceans (10°N?C50°N). The annual mean Q _net climatology shows that oceans obtain the positive Q _net over much of the North Pacific and North Atlantic oceans. Exceptions are the regions of western boundary currents (WBCs) including the Kuroshio and its extension off Japan and the Gulf Stream off the USA and its extension, where oceans release lots of heat into the atmosphere, mainly ascribed to the large surface turbulent heat loss. The statistically significant negative Q _net trends occurred in the WBCs, while the statistically significant positive Q _net trends appeared in the central basins of Northern Subtropical Oceans (CNSOs) including the central basin of Northern Subtropical Pacific and the central basin of Northern Subtropical Atlantic. These indentified Q _net trends, which are independent of both El Ni?o-Southern Oscillation (ENSO) and ENSO Modoki but closely related to global warming forcing, are predominately due to the statistically significant surface latent heat (LH) trends. Over the WBCs, the positive LH trends are mainly induced by the sea surface temperature increasing, indicating the ocean forcing upon overlying atmosphere. In contrast, over the CNSOs, the negative LH trends are mainly caused by the near-surface air specific humidity increasing, indicative of an oceanic response to overlying atmospheric forcing.     

Structural variation of an atmospheric heat source over the Qinghai-Xizang Plateau and its influence on precipitation in northwest China

NCEP/NCAR reanalysis data and a 47-year precipitation dataset are utilized to analyze the relationship between an atmospheric heat source (hereafter called < Q₁ >) over the Qinghai-Xizang Plateau (QXP) and its surrounding area and precipitation in northwest China. Our main conclusions are as follows: (1) The horizontal distribution of < Q₁ > and its changing trend are dramatic over QXP in the summer. There are three strong centers of < Q₁ > over the south side of QXP with obvious differences in the amount of yearly precipitation and the number of heat sinks predominate in the arid and semi-arid regions of northwest China (NWC), beside the northern QXP with an obvious higher intensity in years with less precipitation. (2) In the summer, the variation of the heat source's vertical structure is obviously different between greater and lesser precipitation years in eastern northwest China (ENWC). The narrow heat sink belt forms between the northeast QXP and the southwestern part of Lake Baikal. In July and August of greater precipitation years, the heating center of the eastern QXP stays nearly over 35^oN, and at 400 hPa of the eastern QXP, the strong upward motion of the heating center constructs a closed secondary vertical circulation cell over the northeast QXP (40^o--46^oN), which is propitious to add precipitation over the ENWC. Otherwise, the heating center shifts to the south of 30^oN and disappears in July and August of lesser precipitation years, an opposite secondary circulation cell forms over the northeast QXP, which is a disadvantage for precipitation. Meanwhile, the secondary circulation cell in years with more or less precipitation over the ENWC is also related to the heat source over the Lake Baikal. (3) The vertical structure of the heat source over the western QXP has obvious differences between greater and lesser precipitation years in western northwest China in June and July. The strong/weak heat source over the western QXP produces relatively strong/weak ascending motion and correspondingly constructs a secondary circulation cell in lesser/greater precipitation years.     

A CLOUD-RESOLVING MODELING STUDY OF SURFACE RAINFALL PROCESSES ASSOCIATED WITH LANDFALLING TYPHOON KAEMI(2006)

The detailed surface rainfall processes associated with landfalling typhoon Kaemi(2006) are investigated based on hourly data from a two-dimensional cloud-resolving model simulation. The model is integrated for 6 days with imposed large-scale vertical velocity, zonal wind, horizontal temperature and vapor advection from National Center for Environmental Prediction (NCEP) / Global Data Assimilation System (GDAS) data. The simulation data are validated with observations in terms of surface rain rate. The Root-Mean-Squared (RMS) difference in surface rain rate between the simulation and the gauge observations is 0.660 mm h^-1, which is smaller than the standard deviations of both the simulated rain rate (0.753 mm h^-1) and the observed rain rate (0.833 mm h^-1). The simulation data are then used to study the physical causes associated with the detailed surface rainfall processes during the landfall. The results show that time averaged and model domain-mean P_s mainly comes from large-scale convergence (Q_WVF) and local vapor loss (positive Q_WVT). Large underestimation (about 15%) of P_s will occur if Q_WVT and Q_CM (cloud source/sink) are not considered as contributors to P_s. Q_WVF accounts for the variation of P_s during most of the integration time, while it is not always a contributor to P_s. Sometimes surface rainfall could occur when divergence is dominant with local vapor loss to be a contributor to P_s. Surface rainfall is a result of multi-timescale interactions. Q_WVE possesses the longest time scale and the lowest frequency of variation with time and may exert impact on P_s in longer time scales. Q_WVF possesses the second longest time scale and lowest frequency and can explain most of the variation of P_s. Q_WVT and Q_CM possess shorter time scales and higher frequencies, which can explain more detailed variations in P_s. Partitioning analysis shows that stratiform rainfall is dominant from the morning of 26 July till the late night of 27 July. After that, convective rainfall dominates till about 1000 LST 28 July. Before 28 July, the variations of in rainfall-free regions contribute less to that of the domain-mean Q_WVT while after that they contribute much, which is consistent to the corresponding variations in their fractional coverage. The variations of Q_WVF in rainfall regions are the main contributors to that of the domain-mean Q_WVF, then the main contributors to the surface rain rate before the afternoon of 28 July.     

Measuring C T 2 , C Q 2 , and C TQ in the unstable surface layer,and relations to the vertical fluxes of heat and moisture

The structure parameters of temperature (C _T ² ), humidity (C _Q ² ) and temperature-humidity (C _TQ ) were observed at a height of 4 m in the unstable surface layer using thin platinum wires and two Ly- hygrometers. Two ways of measuring structure parameters were employed: one using spaced sensors, the other using time-delayed observations at one location. It is found that the three structure parameters follow free-convection scaling down to -z/L 0.02. The scaling functions % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₁ (of C _T ² ), % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₂ (of C _TO ) and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₃ (of C _Q ² ) are found to be related through % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₂/% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₁ 0.69 and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₃/% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₁ 0.84. The usefulness of the structure parameters for inferring the fluxes of heat and water vapor, as well as the Bowen ratio, is demonstrated. The scatter is about 30% on either side of the mean.This work was done while the author was a visiting scientist at the Wave Propagation Laboratory, NOAA, ERL, Boulder, U.S.A. He received support from the Netherlands Minister for Science Policy and the U.S. Army Research Office.     

Vertical motion diagnostics related to Q-vectors,and the reorganisation of a cold frontal system

Summary Estimates of vertical motion obtained by explicit solution of the quasi-geostrophic omega equation expressed in terms of Q-vectors (Hoskins et al., 1978) were compared with qualitative estimates based on the divergence of Q-vectors. The horizontal distributions of each were broadly similar, but differed in their detail. At any level, the relative magnitudes of centres of vertical motion varied between the parameters. The precise locations of the main centres also differed. Greater differences existed in the vertical distributions of the parameters. The divergence ofQ andF _Q of Hoskins et al. (1978) both reached maximum values in upper levels, whereas peak values of occurred in middle levels. A vertical distribution similar to was achieved ifF _Q was divided by a stability parameter. The diagnosed vertical velocities appeared consistent with the restructuring of a summertime cold front as it crossed waters south of Australia.With 6 Figures     

Thermodynamic Characteristics of Tropical Cyclones with Rapid Intensity Change over the Coastal Waters of China

In order to investigate the different thermodynamic mechanisms between rapid intensifying (RI) and rapid weakening (RW) tropical cyclones (TCs), the thermodynamic structures of two sets of composite TCs are analyzed based on the complete-form vertical vorticity tendency equation and the NCEP/NCAR reanalysis data. Each composite is composed of five TCs, whose intensities change rapidly over the coastal waters of China. The results show that the maximum apparent heating source Q ₁ exists in both the upper and lower troposphere near the RI TC center, and Q ₁ gets stronger at the lower level during the TC intensification period. But for the RW TC, the maximum Q ₁ exists at the middle level near the TC center, and Q ₁ gets weaker while the TC weakens. The maximum apparent moisture sink Q ₂ lies in the mid troposphere. Q ₂ becomes stronger and its peak-value height rises while TC intensifies, and vice versa. The increase of diabatic heating with height near the TC center in the mid-upper troposphere and the increase of vertical inhomogeneous heating near the TC center in the lower troposphere are both favorable to the TCs’ rapid intensification; otherwise, the intensity of the TC decreases rapidly.     

Anthropogenic heat in the city of S?o Paulo, Brazil

The main goal of this work is to describe the anthropogenic energy flux (Q _F) in the city of S?o Paulo, Brazil. The hourly, monthly, and annual values of the anthropogenic energy flux are estimated using the inventory method, and the contributions of vehicular, stationary, and human metabolism sources from 2004 to 2007 are considered. The vehicular and stationary sources are evaluated using the primary consumption of energy based on fossil fuel, bio fuel, and electricity usage by the population. The diurnal evolution of the anthropogenic energy flux shows three relative maxima, with the largest maxima occurring early in the morning (??19.9 Wm^?2) and in the late afternoon (??20.3 Wm^?2). The relative maximum that occurs around noontime (??19.6 Wm^?2) reflects the diurnal pattern of vehicle traffic that seems to be specific to S?o Paulo. With respect to diurnal evolution, the energy flux released by vehicular sources (Q _FV) contributes approximately 50% of the total anthropogenic energy flux. Stationary sources (Q _FS) and human metabolism (Q _FM) represent about 41% and 9% of the anthropogenic energy flux, respectively. For 2007, the monthly values of Q _FV, Q _FS, Q _FM, and Q _F are, respectively, 16.8?±?0.25, 14.3?±?0.16, 3.5?±?0.03, and 34.6?±?0.41?MJ?m^?2?month^?1. The seasonal evolution monthly values of Q _FV, Q _FS, Q _FM, and Q _F show a relative minimum during the summer and winter vacations and a systematic and progressive increase associated with the seasonal evolution of the economic activity in S?o Paulo. The annual evolution of Q _F indicates that the city of S?o Paulo released 355.2?MJ?m^?2?year^?1 in 2004 and 415.5?MJ?m^?2?year^?1 in 2007 in association with an annual rate of increase of 19.6?MJ?m^?2?year^?1 (from 2004 to 2006) and 30.5?MJ?m^?2?year^?1 (from 2006 to 2007). The anthropogenic energy flux corresponds to about 9% of the net radiation at the surface in the summer and 15% in the winter. The amplitude of seasonal variation of the maximum hourly value of the diurnal variation increases exponentially with latitude.     

The Influence of NaBr/NaCl Ratio on the Br--Catalysed Production of Halogenated Radicals

The activation of Br^- and Cl^- to atomic Br and Cl in sea-spray aerosol was investigated in smog-chamber experiments. In the presence of O₃, hydrocarbons and NaCl aerosol alone no activation was observed. By adding Br^- to the aerosol, the chain reaction: Br + O₃ BrO, BrO + HO₂ HOBr, HOBr HOBr(aq), HOBr(aq) + H⁺ + Br^- Br₂ (6), HOBr(aq) + H⁺ + Cl^- BrCl (7) was verified. The step from reaction (6) to (7) is accompanied by a decrease of the Br^-/Cl^- ratio from 1/600 to less than 1/2000. In the absence of sulphate, the chain is initiated by the reaction of OH(aq) with Br^-. The pH value decreases to less than 2 during the first minutes of the experiment and later on to almost 1 (in the absence of NO_x or SO₂). This is caused by the formation of oxalic acid from alkanes and toluene. In stopped flow experiments, the reduction of Br₂ by oxalic acid was observed to occur through a two-step mechanism: HC₂O₄ ^- + Br₂ Br^- + BrC₂O₄H (k₂₂, k_-22), BrC₂O₄H Br^- + H⁺ + 2 CO₂ (23) with the following rate constants and ratios of rate constants, k ± 2: k₂₂k_-23 / k_-22 = (2.9 ± 0.3) · 10^-4 s^-1, k_-22 / k_-23 = 7000 ± ₃₀₀₀ ¹³⁰⁰⁰ M^-1, k₂₂ = 2 ±_-1 ⁴ M^-1 s^-1, and k_-23 > 0.1 s^-1, k_-22 > 600 M^-1 s^-1. Oxalic acid may be responsible for the inhibition of the chain reaction observed at the end of the experiments.