A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.     

Turbulence decay in stratified and homogeneous marine layers

A spectral approach is applied to shear-induced turbulence in stratified layers. A system of spectral equations for stationary balance of turbulent energy and temperature variances was deduced in the vicinity of the local shear scale L_U = (ε/U_Z³)^1/2. At wavenumbers between the inertial-convective (k^−5/3) and wak turbulence (k⁻³) subranges, additional narrow spectral intervals—‘production’ subranges—may appear (E k⁻¹, E_T k⁻²). The upper boundary of these subranges is determined as L_U, and the lower boundaries as L_R (ε/U_ZN²)^1/2(χ/T_Z²). It is shown that the scale L_U is a unique spectral scale that is uniform up to a constant value for every hydrophysical field. It appears that the spectral scale L_U is equivalent to the Thorpe scale L_Th for the active turbulence model. Therefore, if turbulent patches are generated in a background of permanent mean shear, a linear relation between temperature and mass diffusivities exists. In spectral terms, the fossil turbulence model corresponds to the regime of the Boldgiano-Obukhov buoyancy subrange (E k^−11/5, E_T k^−7/5). During decay the buoyancy subrange is expanded to lower and higher wavenumbers. At lower wavenumbers the buoyancy subrange is bounded by L_** = 3(χ^1/2/N^1/2T_Z), which is equivalent to the Thorpe scale L_Th. In such a transition regime only, when the viscous dissipation rate is removed from the set of main turbulence parameters, the Thorpe scale does not correlate with the buoyancy scale L_N ε^1/2/N^3/2 and fossil turbulence is realized. Oceanic turbulence measurements in the equatorial Pacific near Baker Island confirm the main ideas of the active and fossil turbulence models.     

Prediction of the abiotic degradability of organic compounds in the troposphere

A statistically relevant correlation between the reaction rate coefficient, k _OH, for the OH radical reaction with 161 organic compounds in the gas phase at 300 K, and the corresponding vertical ionisation energies E _i,v, reveals two classes of compounds: aromatics where –log(k _OH/cm³s^-1)3/2E _i,v(eV)–2 and aliphatics where –log(k _OH/cm³s^-1)4/5E _i,v(eV)+3. The prediction of the rate coefficient, k _OH, for the reaction of OH with organic molecules from the above equations has a probability of about 90%. Assuming a global diurnal mean of the OH radical concentration of 5×10⁵ cm³, the upper limit of the tropospheric half-life of organic compounds and their persistence can be estimated.     

Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-eddy Simulations. Part 1: Free Convection

Convection in a quasi-steady, cloud-free, shear-free atmospheric boundary layer is investigated based on a large-eddy simulation model. The performed tests indicate that the characteristic (peak) values of statistical moments at the top of the mixed layer are proportional to the interfacial scales (from gradients of scalars in the interfacial layer). Based on this finding a parameterization is proposed for profiles of scalar variances. The parameterization employs two, semi-empirical similarity functions F_m(z/z_i) andF_i(z/z_i), multiplied by a combination of the mixed-layer scales and the interfacial scales.     

Three-dimensional modeling of synthetic cold fronts approaching the Alps

Summary A numerical model was used to study the behaviour of prototype cold fronts as they approach the Alps. Two fronts with different orientations relative to the Alpine range have been considered. One front approaches from west, a second one from northwest. The first front is connected with southwesterly large-scale air-flow producing pre-frontal foehn, whereas the second front is associated with westerly largescale flow leading to weak blocking north of the Alps.Model simulations with fully represented orography and parameterized water phase conversions have been compared with control runs where either the orography was cut off or the phase conversions were omitted. The results show a strong orographic influence in case of pre-frontal foehn which warms the pre-frontal air and increases the cross-frontal temperature contrast leading to an acceleration of the front along the northern Alpine rim. The latent heat effect was found to depend much on the position of precipitation relative to the surface front line. In case of pre-frontal foehn precipitation only falls behind the surface front line into the intruding cold air where it partly evaporates. In contrary, precipitation already appears ahead of the front in the case of blocking. Thus, the cooling effect of evaporating rain increases the cross-frontal temperature difference only in the first case causing an additional acceleration of the front.List of symbols C _pd specific heat capacity of dry air at constant pressure (C _pd =1004.71 J kg^–1 K^–1) - C _pv specific heat capacity of water vapour at constant pressure (C _pv =1845.96 J kg^–1 K^–1) - C _f propagation speed of a front - x, y horizontal grid spacing (cartesian system) - , horizontal grid spacing (geographic system) - t time step - E turbulent kinetic energy - f Coriolis parameter - g gravity acceleration (g=9.81 ms^–1) - h terrain elevation - H height of model lid (H=9000 m) - k Karman constant (k=0.4) - K _Mh horizontal exchange coefficient of momentum - K _Hh horizontal exchange coefficient of heat and moisture - K _Mz vertical exchange coefficient of momentum - K _Hz vertical exchange coefficient of heat and moisture - l mixing length - l _c specific condensation heat (l _c =2500.61 kJ kg^–1) - l _f specific freezing heat (l _f =333.56 kJ kg^–1) - l _s specific sublimation heat (l _s =2834.17 kJ kg^–1) - longitude - m ₁,m ₂,m ₃ metric coefficients - p pressure - Exner function - Pr Prandtl number - latitude - _M profile function - q _v specific humidity - q _c specific content of cloud droplets - q _i specific content of cloud ice particles - q _R specific content of rain drops - q _S specific content of snow - R _d gas constant of dry air (R _d =287.06 J kg^–1 K^–1) - R _v gas constant of water vapour (R _v =461.51 J kg^–1 K^–1) - r _E radius of earth (r _E =6371 km) - Ri _F flux Richardson number - density of dry air - t time - T temperature - _dia period of diastrophy - potential temperature - _v virtual potential temperature - _e equivalent potential temperature - U relative humidity - u, v, w cartesian wind components - u _F ,v _F front-normal and front-parallel wind components - x, y, z cartesian coordinates - w ^* transformed vertical wind component - W _R speed of falling rain - W _S speed of falling snow - z ^* transformed vertical coordinate Abbreviations GND (above) ground level - MSL (above) mean sea level With 12 Figures     

Decaying stratified turbulence: comparison between a two-point closure EDQNM model and direct numerical simulations

A new statistical two-point closure EDQNM (Eddy Damped Quasi Normal Markovian) model for axisymmetric stratified turbulence is presented. This model takes into account the detailed anisotropic structure of the flow, including angular dependence, and its closure assumption may involve explicit effects of the stratification. In the present paper, a simplified closure assumption is used, in which these explicit effects are not taken into account. The model is tested against three-dimensional direct numerical simulations (DNS) of decaying homogeneous stratified turbulence. A very good agreement is found between DNS results and EDQNM predictions, when energy transfers between the different modes of motion are considered. However, transfers among wavenumbers are insufficiently damped by the EDQNM model, because of the simplified closure assumption, thus yielding a somewhat less satisfactory agreement. An interesting prediction of the EDQNM model is the existence of a ‘directional’ anisotropy at all scales of motion, to the smallest: the level of the total (kinetic + potential) energy density spectrum E(k,θ_k) varies with the angle θ_k that the wavenumber vector k makes with the vertical.     

On the Scaling Laws of the Velocity-Scalar Cospectra in the Canopy Sublayer Above Tall Forests

The scaling laws of the vertical (F _wc ) and longitudinal (F _uc ) velocity-scalar cospectra within the inertial subrange are explored using dimensional arguments and a simplified cospectral budget in the canopy sublayer above three distinct forested ecosystems. The cospectral budget was shown to be consistent with plausible scaling laws originating from dimensional considerations. Using the analytical solution to the novel cospectral budget, it was shown that F _wc (k) and F _uc (k) are governed by the linear superposition of two terms that scale as k ^−2/3−α and k ^−β , where k is the wavenumber, −α is the exponent of the velocity spectrum, and β( ≥ 7/3) depends on the ratio of the similarity constants for the pressure-scalar covariance and the flux transport terms. It was also demonstrated that, when the magnitude of the mean scalar concentration gradient is large, the k ^−2/3−α term dominates the velocity-scalar cospectral budget. For such a case, correcting for biases emerging from high frequency losses in eddy-covariance scalar flux measurements can be readily formulated by using the measured velocity spectral exponent in the inertial subrange.     

土壤数据集对全球陆面过程模拟的影响

基于通用陆面模式（Common Land Model, CoLM）,首次评估了两套最新的全球土壤数据集GSDE（Global Soil Dataset for Earth System Model）和SG（SoilGrids）对全球陆面过程模拟的影响。比较分析了两套数据中砂粒、粘粒、砾石、有机碳的含量和容重这五个土壤属性在全球分布上的差异以及这种差异造成的对模式估计的土壤特性参数、水力热力变量的影响。结果表明,土壤特性参数在全球的空间分布主要受土壤粒径分布（砂粒、粉粒和粘粒）影响,同时也受砾石、有机碳和容重的影响。土壤资料对全球模拟结果影响主要体现在区域差异,对水文学变量的影响（R_e最大达到±100%）大于对土壤热力学变量的影响（R_e<±10%）,对地表辐射变量的影响较小（R_e<±5%）。其中,土壤体积含水量在加拿大中部及西北部、俄罗斯东南部及中西部和澳大利亚中部地区模拟结果相差较大,总径流在低纬地区模拟结果出现较大的差异,热力学变量在非洲北部、加拿大西北部以及俄罗斯中北部的差异稍大。将模拟的土壤体积含水量与站点观测相比,两套数据的表现接近,与站点观测相比都存在一定的偏差,但SG更接近观测,其中在Molly Caren站点（39°57′N,83°27′W）上SG相比GSDE整体提高约0.01～0.02。本研究表明,模式模拟结果受不同土壤数据集的影响显著,可优先考虑诸如SG较准确的土壤数据。土壤属性对陆面模拟的影响需进一步研究。     

Kinetics of the Gas-Phase Reactions of OH Radicals, NO3 Radicals and O3 with Three C7-Carbonyls Formed From The Atmospheric Reactions of Myrcene, Ocimene and Terpinolene

Rate constants for the gas-phase reactions of OH radicals, NO₃ radicals and O₃ with the C₇-carbonyl compounds 4-methylenehex-5-enal [CH₂=CHC(=CH₂)CH₂CH₂CHO], (3Z)- and (3E)-4-methylhexa-3,5-dienal [CH₂=CHC(CH₃)=CHCH₂CHO] and 4-methylcyclohex-3-en-1-one, which are products of the atmospheric degradations of myrcene, Z- and E-ocimene and terpinolene, respectively, have been measured at 296 ± 2 K and atmospheric pressure of air using relative rate methods. The rate constants obtained (in cm³ molecule^–1 s^–1 units) were: for 4-methylenehex-5-enal, (1.55 ± 0.15) × 10^–10, (4.75 ± 0.35) × 10^–13 and (1.46 ± 0.12) × 10^–17 for the OH radical, NO₃ radical and O₃ reactions, respectively; for (3Z)-4-methylhexa-3,5-dienal: (1.61 ± 0.35) × 10^–10, (2.17 ± 0.30) × 10^–12, and (4.13 ± 0.81) × 10^–17 for the OH radical, NO₃ radical and O₃ reactions, respectively; for (3E)-4-methylhexa-3,5-dienal: (2.52 ± 0.65) × 10^–10, (1.75 ± 0.27) × 10^–12, and (5.36 ± 0.28) × 10^–17 for the OH radical, NO₃ radical and O₃ reactions, respectively; and for 4-methylcyclohex-3-en-1-one: (1.10 ± 0.19) × 10^–10, (1.81 ± 0.35) × 10^–12, and (6.98 ± 0.40) × 10^–17 for the OH radical, NO₃ radical and O₃ reactions, respectively. These carbonyl compounds are all reactive in the troposphere, with daytime reaction with the OH radical and nighttime reaction with the NO₃ radical being predicted to dominate as loss processes and with estimated lifetimes of about an hour or less.     

Rapid vertical sampling in stably stratified turbulent shear flows

A new method for obtaining instantaneous vertical profiles of two components of velocity and temperature in thermally stratified turbulent shear flows is presented. In this report, the design and construction of the traversing system will be discussed and results to date will be presented. The method is based on rapid vertical sampling whereby probe sensors are moved vertically at a high speed such that the measurement is approximately instantaneous. The system is designed to collect many measurements for the calculation of statistics such as vertical wave number spectra, mean square vertical gradients, and Thorpe scales. Results are presented for vertical profiles of temperature and compared to vertical profiles measured by single-point Eulerian time averages. The quality of the vertical profiles is found to be good over many profiles. Some comparisons are made between vertical measurements and standard single-point Eulerian measurements for three cases of stably stratified turbulent shear flow in which the initial microscale Reynolds number, Re_λ≈30. In case 1, the mean conditions are characterized by a gradient Richardson number, Ri_g=0.015, for which the flow is “unstable”, meaning the spatially evolving turbulent kinetic energy (E_k) grows. In case 2, Ri_g=0.095, for which the evolving turbulent kinetic energy is almost constant. In case 3, the flow is highly stable, where Ri_g=0.25 and E_k decays with spatial evolution. The measurements indicate anisotropy in the small scales for all cases. In particular, it is found that the ratio grows initially to a maximum and then decays with further evolution. Maximum Thorpe displacements are measured and compared to single-point measures of the vertical scales. It is found that vertical length scales derived from single-point measurements, such as the Ozmidov scale, L_O=(ε/N³)^1/2 and the overturn scale, L_t=θ′/(dT/dz), do not represent well the wide range of overturning scales which are actually present in the turbulence.     

WRF4.0模式陆气耦合对不同物理参数化方案的敏感性初探

模式的不同参数化方案会对模拟的天气和气候产生影响,但如何影响陆气耦合还不是很清楚。利用WRF 4.0模式设计了一组16种不同参数化方案组合的集合试验,对华北地区2013年夏季气候进行模拟。结果表明,不同积云对流参数化方案对模拟的降水平均态影响较大,其次是短波辐射方案,而微物理方案影响较小。进一步利用潜热通量(LH,Latent Heat flux)和土壤湿度(SM,Soil Moisture)、抬升凝结高度(LCL,Lifting Condensation Level)的相关系数R_(SM,LH)和R_(LH,LCL)表征陆气耦合过程,评估了陆气耦合对不同参数化方案的敏感性。模式模拟的R_(SM,LH)在不同区域随着平均降水、相对湿度的增加以及平均短波辐射的减少而减小,而R_(LH,LCL)的变化正好相反。它们内在的物理机制类似,与气候态的干湿变化密切相关,即随着土壤湿度趋于饱和,蒸散发过程逐渐受太阳辐射的限制,陆面对大气的影响减弱。另外,陆气耦合强度受积云对流方案、微物理方案和短波辐射方案的影响而产生差别,这与不同参数化方案模拟的气候态的不同密切相关,因此本研究对于WRF 4.0模式中参数化方案特别是积云对流方案的选择同样有着重要的指导意义。     

Nonlinear interactions of spherical Rossby modes

Weakly nonlinear triad interactions between spherical Rossby harmonics are studied. Each of the harmonics has the form AP_n^m(sin θ)exp[i(mλ − σt)], where A is an amplitude and P_n^m is the associated Legendre function. Equations for the amplitudes are derived and resonance conditions are analysed. The resonance conditions differ substantially from the usual resonance conditions on a β-plane and include a Diophantine equation and a few inequalities for the integer wavenumbers n and m of the interacting modes. Particular analytical series of solutions to the resonance conditions are constructed, which show that modes with arbitrary large wavenumbers can participate in the interactions. A numerical study of the resonance conditions reveals that no more than 21% of the Rossby harmonics can participate in the triad interactions and that chains of the interacting triads soon break. Thus precise interactions (for which the resonance conditions hold exactly) do not result in any significant redistribution of energy over the spectrum. On the other hand, approximate interactions (for which the resonance conditions hold approximately) generate an intensive energy redistribution among short Rossby modes with typical scales much smaller than the Earth's radius. Thus the energy cascade is concentrated mainly in the short-wave part of the spectrum and is very weak in the large-scale domain. The efficiency of the triad interaction of Rossby modes with scales much smaller than the Earth's radius depends strongly on the existence of the so-called interaction latitude at which the local wave-vectors and frequencies of the interacting modes satisfy resonance conditions for plane Rossby waves on the β-plane approximating the neighbourhood of the latitude. If the interaction latitude exists, the interaction is intensive; in the opposite case, it is weak.     

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.     

A Parameterization of Bowen Ratio with Respect to Soil Moisture Availability

The Bowen ratio (B) is impacted by 5 environmental elements: soil moisture availability, m, the ratio of resist-ances between atmosphere and soil pores, ra/rd, atmospheric relative humidity, h, atmospheric stability, ΔT, and environment temperature. These impacts have been investigated over diverse surfaces, including bare soil, free water surface, and vegetation covered land, using an analytical approach. It was concluded that: (a) B is not a continuous function. The singularity exists at the condition αhcb=h, occurring preferably in the following conditions: weak turbulence, stable stratified stability, dry soil, and humid air, where hcb, defined by Eq.(11) is a critical variable. The existence of a singularity makes the dependence of B on the five variables very complicated. The value of B approaches being inversely proportional to m under the conditions m≥mfc (the soil capacity) and / or ra/rd→0. The proportional coefficient changes with season and latitude with relatively high values in winter and over the poles; (b) B is nearly independent of ra/rd during the day. The impact of m on B is much larger as compared to that of ra/rd on B, (c) when h increases, the absolute value of B also increases; (d) over bare soil, when the absolute surface net radiation increases, the absolute value of B will increase. The impact of RN on B is larger at night than during the day, and (e) over plant canopy, the singularity and the dependcies of B on m, ra , and h are modified as compared to that over bare soil. Also (i) during the daytime unstable condition, m exerts an even stronger impact on B, at night, however, B changes are weak in response to the change in m; (ii) the value of B is much more sensitive in response to the changes of turbulent intensity; (iii) the B response to the variation of h over a vegetation covered area is weaker; and (iv) the singularity exists at the condition hcp=h instead of αhcb=h as over bare soil, where hcp is defined by Eq.(49). The formulas derived over bare soil also hold the same when applied to free water bodies as long as they are visualized as a special soil in which the volumetric fraction of soil pore is equal to one and are fully filled with water. Finally, the above discussions, are used to briefly study the impact on the thermally induced mesoscale circulations.     

Free convection frequency spectra of atmospheric turbulence over the sea

Frequency spectra of atmospheric turbulenceS (f) in the inertial subrange are considered in the free convection regime over the sea surface in a case of motionless instrument measurements (Eulerian frequency spectra). The frequency spectra formulaef _* S (f)/ ² =c (f _*/f)^5/3 for wind velocity (=1–3), temperature (=t) and humidity (=e) fluctuations are derived on the basis of similarity theory and the –5/3 law. These relations also can be derived from a consideration of convective large-scale advection of small eddies. The frequency scalef _* = (N ^{1 2}/)^1/2 (H/z ²)^1/3 is the lower bound of the inertial subrange and it is of order 10^–2 Hz.The spectra formulae are compared with direct measurements of atmospheric turbulence from the fixed research tower in the coastal zone of the Black Sea in calm weather. It is shown that these formulae are realized at least over two to three decades of the frequency range (approximately from 10^–2 to 10 Hz) and values of the numerical coefficients are found. The derived formulae can be used for calculations of sensible and latent heat fluxes by measuring the high-frequency range of spectra at a fixed point at low wind speeds when the conventional inertial dissipation method is not applicable.     

Determination of rate constants for reactions of some hydrohaloalkanes with OH radicals and their atmospheric lifetimes

Rate constants have been measured for the gas-phase reactions of hydroxyl radical with partly halogenated alkanes using the discharge-flow-EPR technique over the temperature range 298–460 K. The following Arrhenius expressions have been derived (units 10^–13 cm³ molecule^–1 s^–1): (8.1 _–1.2 ^+1.5 ) exp{–(1516±53)/T} for CHF₂Cl (HCFC-22); (10.3 _–1.5 ^+1.8 ) exp{–(1588±52)/T} for CH₂FCF₃ (HFC-134a); (11.3 _–1.6 ^+2.1 ) exp{–(918±52)/T} for CHCl₂CF₂Cl (HCFC-122); (9.2 _–2.0 ^+2.5 ) exp{–(1281±85)/T} for CHFClCF₂Cl (HCFC-123a).The atmospheric lifetimes for the substances have been estimated to be 12.6, 12.9, 1.05, and 4.8 years, respectively, and the accuracy of the estimates is discussed.     

Internal tide generation at a continental shelf/slope junction: A comparison between theory and a laboratory experiment

Observations of internal tide generation over continental slopes in a laboratory experiment have been carried out, with the objectives of making comparisons with linear generation theory and investigating its limitations. Both continuous and layered stratification have been considered. A measure of the amplitude of the barotropic tidal forcing (and hence of non-linearity) is given by the Froude Number F = u_sb/c_w, where u_sb is the maximum barotropic tidal velocity at the shelf break, and c_w is the long-wave speed of the lowest internal mode.For continuous stratification, good agreement was obtained for “steep” slopes (α/c > 1), where α is the slope at the continental slope and c is the slope of the internal wave rays of tidal frequency), even for quite large amplitude motions (F < 1.6), and the upper limit of its quantitative usefulness was not reached. For “flat” slopes (α/c < 1) reasonable agreement was also obtained, even up to quite amplitudes (F < 3.1), although some departure from linear theory was apparent.For two-layer flows the applicability of linear theory was much more restricted. For F 0.5 there was only qualitative agreement and for larger F (>1) significant differences were observed. The latter were due to the substantial advection and associated hydraulic jumps which occured seaward of the shelf-break during the ebb-phase of the barotropic tide. Shelf-break values of F > 1 are common in the ocean.     

A theory for the scalar roughness and the scalar transfer coefficients over snow and sea ice

Although the bulk aerodynamic transfer coefficients for sensible (C _H ) and latent (C _E ) heat over snow and sea ice surfaces are necessary for accurately modeling the surface energy budget, they have been measured rarely. This paper, therefore, presents a theoretical model that predicts neutral-stability values of C _H and C _E as functions of the wind speed and a surface roughness parameter. The crux of the model is establishing the interfacial sublayer profiles of the scalars, temperature and water vapor, over aerodynamically smooth and rough surfaces on the basis of a surface-renewal model in which turbulent eddies continually scour the surface, transferring scalar contaminants across the interface by molecular diffusion. Matching these interfacial sublayer profiles with the semi-logarithmic inertial sublayer profiles yields the roughness lengths for temperature and water vapor. When coupled with a model for the drag coefficient over snow and sea ice based on actual measurements, these roughness lengths lead to the transfer coefficients. C _E is always a few percent larger than CH. Both decrease monotonically with increasing wind speed for speeds above 1 m s^–1, and both increase at all wind speeds as the surface gets rougher. Both, nevertheless, are almost always between 1.0 × 10^–3 and 1.5 × 10^–3.     

Latitudinal variation of measured NO2 photolysis frequencies over the Atlantic Ocean between 50° N and 30° S

Using a filter radiometer, the meridional profile of the NO₂ photolysis frequency, J(NO₂), was measured between 50° N and 30° S during the cruise ANTVII/1 September/October 1988 of the research vessel Polarstern on the Atlantic Ocean. Simultaneously, global broadband irradiance and acrosol were monitored. Clean marine background air with low aerosol loads (b _sp=(1–2)×10^-5 m^-1) was encountered at the latitudes 25° N–30° N and 18° S–27° S, respectively. Under these conditions and an almost cloudless sky J(NO₂) reached 7.3×10^-3 s^-1 (2 sr) for a zenith angle of 30°. Between 30° N and 30° S, the latitudinal variation of the J(NO₂) noontime maxima was less than ± 10%, while the mean value at noon was 7.8×10^-3 s^-1. For the set of all data between 50° N and 30° S, a nearly linear correlation of J(NO₂) vs. global broadland irradiance was found. The slope of (8.24±0.03)×10^-5 s^-1/mW cm^-2 agrees within 10% with observations in Jülich (51° N, 6.2° E).     

Sensitivity Tests of an Energy Balance Model to Choice of Stability Functions and Measurement Accuracy

Sensible (H) and latent (L_vE) heat fluxes are obtained by a combined energy budget – similarity model applied to observations from Melle in Belgium and Cabauw in The Netherlands. The sensitivity to both the stability functions and the accuracy of input data is investigated. In a first step, fluxes are calculated for a selection of stability functions and compared to values obtained with pre-defined (reference) functions. For the diurnal fluxes higher than 10 W m⁻² in 1996 at Melle, the root-mean-square rmsreaches 9 W m⁻² for H and 6 W m⁻² for L_vE, depending on the chosen functions. A lesser sensitivity is obtained at Cabauw and can be explained by lower absolute values of the stability parameter ζ (L involving the Obukhov length) mainly induced by higher mean wind speeds. Different stability bins are also considered. It is concluded that a more accurate assessment of the stability functions is already desirable for absolute values of L above a few metres. These values are not so scarce at Melle and should be captured in the future by an increasing number of new developing long-term measurement stations. In a second step, a statistical approach is proposed with errors depicted by both systematic biases and random fluctuations represented by means of Gaussian distributions. The results show that very accurate measurements are needed in order to maintain the mean annual value of the bias and rms below 5–10 W m⁻², and thus to allow the discrimination between the sensitivity to errors on input data and to the stability functions selection.