Development and Validation of an Evaporation Duct Model.Part Ⅰ:Model Establishment and Sensitivity Experiments

Based on the Coupled Ocean-Atmospheric Response Experiment（COARE）bulk algorithm and the Naval Postgraduate School（NPS）model,a universal evaporation duct（UED）model that can flexibly accommodate the latest improvements in component（such as stability function,velocity roughness,and scalar roughness）schemes for different stratification and wind conditions,is proposed in this paper.With the UED model,the sensitivity of the model-derived evaporation duct height（EDH）to stability function（Ψ）,ocean wave effect under moderate to high wind speeds,and scalar roughness length parameterization,is investigated,and relative contributions of these factors are compared.The results show that the stability function is a key factor influencing the simulated EDH values.Under unstable conditions,the EDH values from stability functions of Fairall et al.（1996）and Hu and Zhang（1992）are generally higher than those from others;while under stable conditions,unreasonably high EDHs can be avoided by use of the stability functions of Hu and Zhang（1992）and Grachev et al.（2007）.Under moderate to high wind speeds,the increase in velocity roughness length z₀ due to consideration of the true ocean wave effect acts to reduce modeled EDH values;this trend is more pronounced under stable conditions.Although the scalar roughness length parameterization has a minor effect on the model-derived EDH,a positive correlation is found between the scalar roughness length z_0qand the model-derived EDH.     

现场大气稳定度分类方法探讨

大气稳定度是反映湍流活动的主要指标,恰当地选择大气稳定度分类方法是大气环境影响评价的基本问题之一。目前国内最常用的P-C法(根据P-T法修改成适合我国应用的形式)有较强的经验性,分类精度不高。大气稳定度实质上是大气热力和动力过程对湍流产生、发展或抑制能力的一种客观量度,因此理查逊数和莫宁-奥布霍夫长度等指标将比P-C法能更客观地表征大气稳定度状况。本文利用野外实测资料,给出总体理查逊数划分大气稳定度级别的方法和标准,以及不同分类方法结果比较。     

Estimating crop top microclimates from weather station data

Abstract

A mathematical model (Microsim) was developed to estimate the microclimate at the top of nearby crops using inputs of weather station data and some knowledge about crop characteristics, such as height, albedo, and leaf area index. The model was tested using data measured simultaneously over a weather station and over each of two crops ‐ corn and soybean. Temperatures at the top of unstressed, uniform crops on level terrain within 1600 m of a recording weather station were estimated within 1.0° C 96% of the time for a corn crop and 92% of the time for a soybean crop. Winds at crop top were estimated within 0.4 m s^?1 92% of the time for corn and 100% of the time for soybean. Energy balance flux density estimates for the corn crop resulted in correlation coefficients of r > 0.89 for each of R_n, LE, H and G. microsim worked well under atmospheric conditions that ranged from very stable to unstable.

An enhancement was made to the model to describe wind and temperature profiles based on the complete fetch characteristics of the sites. This resulted in significantly better wind estimates, but had the disadvantage of requiring more information about the crop and weather station surroundings.     

Hybrid Methods for Computing the Streamfunction and Velocity Potential for Complex Flow Fields over Mesoscale Domains

Three types of previously used numerical methods are revisited for computing the streamfunction ψ and velocity potential χ from the horizontal velocity v in limited domains. The first type, called the SOR-based method, uses a classical successive over-relaxation (SOR) scheme to compute ψ (or χ) first with an arbitrary boundary condition (BC) and then χ (or ψ) with the BC derived from v . The second type, called the spectral method, uses spectral formulations to construct the inner part of (ψ, χ)—the inversion of (vorticity, divergence) with a homogeneous BC, and then the remaining harmonic part of (ψ, χ) with BCs from v . The third type, called the integral method, uses integral formulas to compute the internally induced (ψ, χ)—the inversion of domain-internal (vorticity, divergence) using the free-space Green?s function without BCs and then the remaining harmonic ψ (or χ) with BCs from v minus the internally-induced part. Although these methods have previously been successfully applied to flows in large-scale and synoptic-scale domains, their accuracy is compromised when applied to complex flows over mesoscale domains, as shown in this paper. To resolve this problem, two hybrid approaches, the integral-SOR method and the integral-spectral method, are developed by combining the first step of the integral method with the second step adopted from the SOR-based and spectral methods, respectively. Upon testing these methods on real-case complex flows, the integral-SOR method is significantly more accurate than the integral-spectral method, noting that the latter is still generally more accurate than the three previously-used methods. The integral-SOR method is recommended for future applications and diagnostic studies of complex flows.     

Flux-Variance Method for Latent Heat and Carbon Dioxide Fluxes in Unstable Conditions

Applied previously to momentum and heat fluxes, the present study extends the flux-variance method to latent heat and CO₂ fluxes in unstable conditions. Scalar similarity is also examined among temperature (θ), water vapour (q), and CO₂ (c). Temperature is adopted as the reference scalar, leading to two feasible strategies to estimate latent heat and CO₂ fluxes: the first one relies on flux-variance similarity relations for scalars, while the second is based on the parameterization of relative transport efficiency in terms of scalar correlation coefficient and a non-dimensional quantity. The relationship between the θ-to-q transport efficiency (λ _{θ q} ) and θ-q correlation coefficient (R _{θ q} ) is used to describe the intermediate hydrological conditions. We also parameterize the θ-to-c transport efficiency (λ _{θ c} ) as a function of the θ-c correlation coefficient (R _{θ c} ) by introducing a new non-dimensional ratio (α). The flux-variance method is a viable technique for flux gap-filling, when turbulence measurements of wind velocity are not available. It is worth noting that the extended method is not exempt from a correction for density effects when used for estimating water or carbon exchange.     

Temperature–Humidity Dissimilarity and Heat-to-water-vapour Transport Efficiency Above and Within a Pine Forest Canopy: the Role of the Bowen Ratio

Over the past 15 years atmospheric surface-layer experiments over heterogeneous canopies have shown that the vertical transfer of sensible heat and water vapour exhibit a strong dissimilarity. In particular, the sensible-heat-to-water-vapour transport efficiencies generally exceed unity. One of the main consequences is that evaporation (latent heat flux) computed by the flux-variance method is overestimated, as persistently demonstrated by comparisons with evaporation obtained with the eddy-correlation method. Various authors proposed to take into account the temperature–humidity dissimilarity to extend the applicability of the flux-variance method in order to compute evaporation from non-uniform surfaces. They attempted to connect the sensible-heat-to-water-vapour transport efficiency (λ) to the correlation coefficient between temperature and humidity turbulent fluctuations (R _Tq ). This approach was found to be successful over ‘wet’ surfaces for which λ can be approximated by R _Tq and ‘dry’ surfaces for which λ can be approximated by 1/R _Tq . However, no solution has been proposed until now for intermediate hydrological conditions. We investigated this question using eddy-correlation measurements above and inside a pine forest canopy. For both levels, our data present a strong likeness with previously published results over heterogeneous surfaces. In particular, they confirm that λ is R _Tq in wet conditions and 1/R _Tq in dry conditions. Moreover, we defined the range of the Bowen ratio (Bo) values for which those two approximations are valid (below 0.1 and greater than 1, respectively) and established a relationship between λ, R _Tq and Bo for the intermediate range of Bo. We are confident that this new parameterization will enlarge the applicability of the flux-variance method to all kinds of heterogeneous surfaces in various hydrological conditions     

Variation on PAR to global solar radiation ratio along altitude gradient in Naeba Mountain

Summary We investigated the ratio of photosynthetically active photon flux (Q _p ) to global solar radiation (R _s ) at three sites along different altitudes in Naeba Mountain, Japan at various temporal scales based on 3 years measurement data (1999–2001). The lowest values of the ratio ever reported were found for all sites on both an hourly and a daily scale. A similar slight diurnal pattern was found for all sites based on the monthly mean hourly values of the ratio. However, different sites exhibited different seasonal courses. Statistically significant altitude dependencies were found for the hourly Q _p /R _s under both clear and cloudy weather conditions (t-test, P < 0.001). For clear weather conditions, the hourly Q _p /R _s exhibited an increasing trend with altitude at an average rate of 3.6% per km. The increasing rate was more noticeable below 900 m (8.8% per km) than above (0.7% per km). The inverse trend was found for hourly Q _p /R _s with the altitude under cloudy weather conditions. The hourly Q _p /R _s decreased from 550 m to 1500 m at a rate of 1.8% per km. Again, a major decrease occurred below 900 m, which had the rate of 4.2% per km, compared with 0.2% per km over 900 m. Although the same tendencies were noted for daily Q _p /R _s , under clear sky conditions, they were not as statistically significant as the hourly counterpart (t-test, P < 0.021). The increasing rate of Q _p /R _s at this scale under clear weather conditions was near that of the hourly rate, but below the 900 m rate was reduced to near half of the hourly rate (4.7% per km). And the rate over 900 m increased to 2.7% per km. On the other hand, statistically significant altitude effect was noted for the daily under cloudy weather conditions (t-test, P < 0.002). A rapidly decreasing rate was found for it along the altitude gradient. The reason was due to the large decreasing rate below 900 m (12.4% per km). But again similar to the hourly Q _p /R _s , a very small decreasing rate of daily Q _p /R _s was found over 900 m under cloudy weather conditions. These results suggest the necessity of considering the altitude dependency of Q _p /R _s in future studies.     

Regional shear stress of broken forest from radiosonde wind profiles in the unstable surface layer

Mean wind speed profiles were measured by tracking radiosondes in the unstable atmospheric boundary layer (ABL) over the forested Landes region in southwestern France. New Monin-Obukhov stability correction functions, recently proposed following an, analysis by Kader and Yaglom, as well as the Businger-Dyer stability formulation were tested, with wind speeds in the surface sublayer to calculate the regional shear stress. These profile-derived shear stresses were compared with eddy correlation measurements gathered above a mature forest stand, at a location roughly, 4.5 km from the radiosonde launch site. The shear stress values obtained by means of the newly proposed stability function were in slightly better agreement with the eddy correlation values than those obtained by means of a Businger-Dyer type stability function. The general robustness of the profile method can be attributed in part to prior knowledge of the regional surface roughness (z ₀=1.2 m) and the momentum displacement height (d ₀=6.0 m), which were determined from neutral wind profile analysis. The 100 m drag coefficient for the unstable conditions above this broken forest surface was found to beu _* ² /V ₁₀₀ ² =0.0173.     

The role of raindrop coalescence and breakup in rainfall modeling

In a one dimensional model, the microphysical properties of rainwater content (M), rainfall rate (R) and radar reflectivity factor (Z) are examined in five numerical experiments. In general, neglecting both coalescence and breakup produces less than 10% of the errors in M and R within the top 2 km of a rainshaft (below cloud base). This is true regardless of the rain intensity and downdraft strength for Mashall-Palmer distributions at the top boundary. However, as the distance from the top boundary increases, the errors of M, R and Z calculated without including the collisional processes become larger. These errors tend to increase as the rain intensity increases or the downdraft decreases. The errors due to neglecting breakup are generally larger than those resulting from the exclusion of both coalescence and breakup. At lower levels, the errors for Z are usually very large unless the rainfall is very light (less than 10 mm h⁻¹ ). For heavy rainfall with strong downdraft, coalescence and breakup can be neglected for fall distances greater than 2 km without appreciable errors. On the other hand, for heavy rainfall in a weak downdraft, coalescence and breakup must be included to ensure accuracy.     

Dependence of geostrophic drag on intensity of convection,baroclinicity, and acceleration

This paper examines the practical importance of stability, baroclinicity, and acceleration in the bulk ABL similarity formulations, in light of the random errors inherent in field measurements. This is done by propagating the measurement uncertainties through a theoretical model for the bulk ABL similarity functionsA ₀ andB ₀, under a range of assumed (but always unstable) conditions. It is shown that random measurement errors and acceleration effects may overwhelm most effects of baroclinicity and stability, once conditions are at least slightly unstable. Because of this, it is hard to discern a clear functional dependence ofA ₀ andB ₀ on degree of instability. Thus, for a given value ofh _i/z₀, whereh _i is the inversion height andz ₀ is the surface roughness, the geostrophic drag coefficient, which depends onA ₀ andB ₀, and weakly onh _i/z₀, may also be taken to be nearly independent of degree of instability.     

Measurements and estimations of photosynthetically active radiation in Beijing

The measurements of the photosynthetic photon flux density (Q_p) and other solar components have been in Beijing for 2-year period. The Q_p, broadband solar radiation (R_s) and the PAR fraction (Q_p / R_s) showed similar seasonal features that peaked in value during the Summer and reached their lowest value during the Winter. The PAR fraction ranged from 1.68 E M J^− 1 (Winter) to 1.98 E M J^− 1 (Summer) with an annual mean value of 1.83 E M J^− 1. The analysis of the hourly values also revealed a diurnal pattern, with higher values of Q_p and R_s being observed around noon. The PAR fraction increased from 1.78 to 1.89 μE J^− 1 (hourly values), as the sky conditions changed from clear to cloudy. The monthly mean hourly PAR fraction also revealed a diurnal variation, however, with lower values being observed around noon during most months. In November, the diurnal variations showed an opposite feature in comparison with other months. This is mainly attributed to the diurnal variations in the water vapor concentration.Two models were developed to estimate Q_p from R_s. The models consisted of atmospheric parameters that were found to cause substantial changes to the PAR fraction, such as sky clearness, brightness and path length. The estimated Q_p obtained via different equations was much closer to the observed values, with relative errors below 20% in Beijing. The Q_p and R_s data collected at three stations with featuring different climate types from within Beijing were used for verifying the transferability of the models. The correlation coefficients between the measured and estimated Q_p values decreased at these stations, and the relative error increased. This indicates that the estimation models need to be modified accordingly for the local climatic conditions.     

Derivation of the relationship between the Obukhov stability parameter and the bulk Richardson number for flux-profile studies

Using the relationship between the bulk Richardson numberR _z and the Obukhov stability parameterz/L (L is the Obukhov length), formally obtained from the flux-profile relationships, methods to estimatez/L are discussed. Generally,z/L can not be uniquely solved analytically from flux-profile relationships, and it may be defined using routine observations only by iteration. In this paper, relationships ofz/L in terms ofR _z obtained semianalytically were corrected for variable aerodynamic roughnessz ₀ and for aerodynamic-to-temperature roughness ratiosz ₀/z _T, using the flux-profile iteration procedure. Assuming the so-called log-linear profiles to be valid for the nearneutral and moderately stable region (z/L<1), a simple relationship is obtained. For the extension to strong stability, a simple series expansion, based on utilisation of specified universal functions, is derived.For the unstable region, a simple form based on utilisation of the Businger-Dyer type universal functions, is derived. The formulae yield good estimates for surfaces having an aerodynamic roughness of 10^–5 to 10^–1 m, and an aerodynamic-to-temperature roughness ratio ofz ₀/z _T=0.5 to 7.3. When applied to the universal functions, the formulae yield transfer coefficients and fluxes which are almost identical with those from the iteration procedure.     

Current Observations of Internal Tides and Parametric Subharmonic Instability in Luzon Strait

Acoustic Doppler Current Profiler (ADCP) measurements of the velocity structure in the thermocline in Luzon Strait are presented. The statistics for current, vertical variation of the inertia-gravity waveband, parametric subharmonic instability (PSI), and current shear are analyzed. It was found that 1) barotropic flow primarily consists of a nearly circular mixed tide. Diurnal tides are strongest and show smooth variance with a fortnightly spring-neap cycle, indicative of the astronomical tide-generating force. However, the semi-diurnal band power exhibits a high-frequency oscillation as a result of non-linear interactions. The high-frequency band power with high values during the spring tide oscillates with the tidal cycle. Near-inertial wave motions showing random variance may be caused by changes in the wind forcing at the sea surface or by random forcing. 2) Baroclinic velocities exhibit strong shear structure. The observed large changes in the amplitude of the baroclinic velocity and the limited vertical extent of the high-velocity cores may be interpreted as internal wave beams that pass through the observed water column. Semi-diurnal tides are dominant in the baroclinic velocity. Kinetic energy spectra also revealed that additional peaks were centred at sum-tidal-inertial interaction frequencies (such as M₂ ?+?f) and difference-interaction frequencies (such as M₂ ???f). The spectral exponent of the baroclinic velocity is ω^?α (1?<?α?<?3). 3) Strong non-linear interactions among internal waves exist, and the semi-diurnal (M₂) component plays a key role in these interactions. Bicoherence analysis showed that M₂/2 waves were non-linearly coupled with the dominant M ₂ internal tide. 4) The polarization relations were used to diagnose observational internal tidal motions. Diurnal waves propagate to the east-northeast, and the semi-diurnal–diurnal waves propagate westward. In the case of diurnal tides, the minor to major axis ratio is different from the expected value of f/ω_{K 1} because of the deviation of inclinations, whereas, for semi-diurnal tides, it is close to the expected value of f/ω _M2 at depths from 30 to 150 m.

RÉSUMÉ?[Traduit par la rédaction] Nous présentons les mesures de la structure des vitesses dans la thermocline faites par profileur de courant à effet Doppler (ADCP) dans le détroit de Luçon. Nous analysons les statistiques sur le courant, la variation verticale de la gamme d'ondes d'inertie-gravité, l'instabilité subharmonique paramétrique (PSI) et le cisaillement du courant. Il ressort que 1) l’écoulement barotrope consiste principalement en une marée mixte presque circulaire. Les marées diurnes sont les plus fortes et présentent une variance régulière dans un cycle vives-eaux mortes-eaux de deux semaines qui révèle la nature astronomique de la force qui produit les marées. Cependant, la puissance de bande semi-diurne affiche une oscillation de haute fréquence causée par des interactions non linéaires. La puissance de bande de haute fréquence avec des valeurs élevées durant les vives-eaux oscille avec le cycle des marées. Les mouvements des ondes quasi inertielles affichant une variance aléatoire peuvent être dus à des variations dans le forçage par le vent à la surface de la mer ou à un forçage aléatoire. 2) Les vitesses baroclines exhibent une structure fortement cisaillée. Les changements marqués observés dans l'amplitude des vitesses baroclines et l'extension verticale limitée des c?urs de vitesse élevée peuvent être interprétés comme des faisceaux d'ondes internes qui traversent la colonne d'eau observée. Les marées semi-diurnes sont dominantes dans les vitesses baroclines. Les spectres d’énergie cinétique ont aussi révélé que des crêtes additionnelles étaient centrées aux fréquences d'interaction cumulative marée-inertie (comme à M₂ ?+?f) et aux fréquences d'interaction différentielle (comme à M₂ ???f). L'exposant spectral de la vitesse barocline est ω^?α (1?<?α?<?3). 3) Il existe de fortes interactions non linéaires entre les ondes internes et la semi-diurne (M₂) joue un rôle clé dans ces interactions. L'analyse de bicohérence a montré que les ondes M₂/2 étaient non-linéairement couplées avec la marée interne dominante M₂ . 4) Nous avons utilisé les relations de polarisation pour diagnostiquer les mouvements de marée internes déduits des observations. Les ondes diurnes se propagent vers l'est-nord-est et les ondes semi-diurnes se propagent vers l'ouest. Dans le cas des marées diurnes, le rapport de l'axe secondaire à l'axe principal diffère de la valeur attendue (f/ω_{K 1} ) à cause de l’écart des inclinaisons alors que pour les marées semi-diurnes, il est proche de la valeur attendue (f/ω _M2 ) à des profondeurs de 30 à 150 m.     

Drag and heat transfer relations for the planetary boundary layer

A theory is offered for the drag and heat transfer relations in the statistically steady, horizontally homogeneous, diabatic, barotropic planetary boundary layer. The boundary layer is divided into three regionsR ₁,R ₂, andR ₃, in which the heights are of the order of magnitude ofz ₀,L, andh, respectively, wherez ₀ is the roughness length for either momentum or temperature,L is the Obukhov length, andh is the height of the planetary boundary layer. A matching procedure is used in the overlap zones of regionsR ₁ andR ₂ and of regionsR ₂ andR ₃, assuming thatz ₀ L h. The analysis yields the three similarity functionsA(),B(), andC() of the stability parameter, = u _*/fL, where is von Kármán's constant,u _* is the friction velocity at the ground andf is the Coriolis parameter. The results are in agreement with those previously found by Zilitinkevich (1975) for the unstable case, and differ from his results only by the addition of a universal constant for the stable case. Some recent data from atmospheric measurements lend support to the theory and permit the approximate evaluation of universal constants.     

Intensity Evolution of Zonal Shear Line over the Tibetan Plateau in Summer: A Perspective of Divergent and Rotational Kinetic Energies

Based on the ERA5 reanalysis datasets during 1980–2019, a total of eleven zonal shear lines (ZSLs) that caused heavy precipitation and lasted more than 60 hours over the Tibetan Plateau in summer are selected for composite analysis. By decomposing the kinetic energy (K) near the ZSL into divergent and rotational kinetic energies (K_D and K_R) and the kinetic energy of interaction between the divergent wind and the rotational wind (K_RD), the influence of the rotational and divergent winds on the evolution of the ZSL intensity is investigated from the perspective of K_D and K_R. The main results are as follows. The ZSL is a comprehensive reflection of rotation and convergence. The intensity evolution of ZSL is essentially synchronized with those of K, K_R, and K_RD but lags behind K_D by about three hours. The enhancement of K is mainly contributed by K_R, which is governed by the conversion from K_D to K_R. Furthermore, the increase in the conversion from K_D to K_R is controlled by the geostrophic effect term Af, which is determined by the joint enhancement of the zonal rotational and meridional divergent wind components (u_R and v_D). Therefore, the joint enhancement of u_R and v_D controls the increase of the ZSL intensity, leading to increased precipitation.     

珠海凤凰山森林下垫面干季和湿季气象要素的对比分析与动量和感热交换系数的参数化研究

利用珠海凤凰山陆气相互作用观测塔站2014年11月至2016年5月的观测数据,对比分析了干湿季森林下垫面能量通量和气象要素的变化特征,分析了在不同稳定度下3个风向范围（315°～45°、45°～135°和135°～225°）的动量和感热交换系数随冠层表面风速的变化特征,并对动量和感热交换系数进行了参数化研究。结果表明:干季感热和潜热通量值相当,湿季潜热远大于感热。干季和湿季的夜晚都出现负感热现象,感热从大气向森林输送。相对湿度的变化幅度大,与该地气象状况密切相关,相对湿度的垂直梯度夜晚较大,白天较小。干季的气温垂直梯度比湿季的明显。风速在冬季变化平缓,夏季变化剧烈,低层风速随高度变化梯度明显,高层较紊乱。各高度风向差异不大。中性和近中性状态下,在风向为315°～45°、45°～135°和135°～225°时,动量交换系数C_dn分别为0.05、0.0055和0.022,感热交换系数C_hn分别为0.0055、0.003和0.004。在稳定和不稳定状态下,动量交换系数C_d、感热交换系数C_h随冠层表面风速v明显发生变化,稳定条件下,C_d、C_h随v的增大而增大;不稳定条件下,C_d、C_h随v的增大而减小。分不同风向对森林冠层C_d、C_h在稳定和不稳定条件下与v的关系进行了拟合,得到了参数化公式。     

Anomalies in Flux-Gradient Relationships Over Forest

Simultaneous profile and eddy correlation flux data gathered over Thetford Forest, U.K., have been analysed to find values of the vertical turbulent diffusivities K _M, K _H and K _E (for momentum, heat and water vapour transfer, respectively) at a reference height z _R, nine roughness lengths above the zero-plane displacement d. The results show: (i), that values of K _M over the forest are not significantly different from these predicted by semiempirical diabatic influence functions appropriate to much smoother surfaces such as short grass; and (ii), that K _H and K _E exceed their values predicted from the semiempirical functions by an average factor of 2 or more in unstable, near neutral and slightly stable conditions. These conclusions are strongly dependent on the assumed behaviour of d, here taken as 0.76 tree heights, independent of both property and stability. Consideration is given to an alternative analysis procedure, in which values of the zero-plane displacements d _H and d _E for heat and water vapour respectively, are obtained from the data by assuming K _H and K _E to be given by semiempirical diabatic influence functions; this procedure is shown to be unacceptable on both practical and physical grounds. To account for the anomalies in K _H and K _E, a mechanism is proposed in which the horizontally inhomogeneous temperature structure of the canopy causes free convection to be maintained by discrete; localized heat sources and/or sinks, effectively enhancing turbulent transport processes even in near neutral conditions.     

Precipitation drives interannual variation in summer soil respiration in a Mediterranean-climate, mixed-conifer forest

Predictions of future climate change rely on models of how both environmental conditions and disturbance impact carbon cycling at various temporal and spatial scales. Few multi-year studies, however, have examined how carbon efflux is affected by the interaction of disturbance and interannual climate variation. We measured daytime soil respiration (R _s) over five summers (June–September) in a Sierra Nevada mixed-conifer forest on undisturbed plots and plots manipulated with thinning, burning and their combination. We compared mean summer R _s by year with seasonal precipitation. On undisturbed plots we found that winter precipitation (PPT_w) explained between 77–96% of interannual variability in summer R _s. In contrast, spring and summer precipitation had no significant effect on summer R _s. PPT_w is an important influence on summer R _s in the Sierra Nevada because over 80% of annual precipitation falls as snow between October and April, thus greatly influencing the soil water conditions during the following growing season. Thinning and burning disrupted the relationship between PPT_w and R_s, possibly because of significant increases in soil moisture and temperature as tree density and canopy cover decreased. Our findings suggest that R _s in some moisture-limited ecosystems may be significantly influenced by annual snowpack and that management practices which reduce tree densities and soil moisture stress may offset, at least temporarily, the effect of predicted decreases in Sierran snowpack on R _s.     

Phillips模式的斜压不稳定的非线性饱和问题:能量情形

A conservation law for the Phillips model is derived. Using this law, the nonlinear saturation of purely baroclinic instability caused by the vertical velocity shear of the basic flow in the Phillips model-the case of energy-is studied within the context of Arnold's second stability theorem. Analytic upper bounds on the energy of wavy disturbances are obtained. For one unstable region in the parameter plane, the result here is a second-order correction in ε to Shepherd's; For another unstable region, the analytic upper bound on the energy of wavy disturbances offers an effective constraint on wavy (nonzonal) disturbances φ'i at any time.     

An empirical expression for aerodynamic resistance in the unstable boundary layer

In unstable conditions, the set of equations defining the aerodynamic resistance to sensible heat transfer, r _a , cannot be solved analytically. An iterative technique must be used to obtain r _a exactly, but this is cumbersome and time consuming. In this paper, a new, empirical equation is presented relating the ratio, Q, of the aerodynamic resistances in neutral and unstable conditions, to the bulk Richardson number, Ri _B . The equation takes the form Q = a + b(–Ri) ^c , where a, b and c are empirical functions of (z – d)/z _om . This model is shown to predict r _awith a mean absolute error of 0.06 s m^–1 over the ranges -15 < Ri _B < 0 and 10 < (z – d)/z _om < 2300. Statistical comparison with other equations that have been proposed for r _a in unstable conditions indicates the superior precision of the model presented here.