Validation of parametrisations for the meridional energy and moisture transport used in simple climate models

Parametrisations of meridional energy and moisture transport used in zonally averaged climate models are validated using reanalysis data and results from a doubling CO₂-experiment from a general circulation model. Global meridional fluxes of moisture and sensible heat are calculated by integrating surface and top-of-the-atmosphere vertical fluxes from one pole to the other. The parametrisations include an eddy-diffusion term, representing down-gradient transport of specific humidity and temperature due to the transient atmospheric eddies at mid- and high latitudes, and simple representations of the mean meridional circulation. Qualitative and quantitative agreement between the increased hydrological cycle in the 2×CO₂-run from the GCM and the parametrisation is found. The performance for the sensible heat flux shows larger differences to the GCM results, particularly at low latitudes. Seasonal variations of the moisture and sensible heat transport are well captured by parametrisations including the influence of the mean meridional circulation. Interannual variability cannot be simulated. An examination of the parametrisations on different spatial scales suggests that they should not be used for small scales. Furthermore, two closures for the zonal distribution of precipitation were examined. They are used in zonally averaged atmosphere models coupled to an ocean model with different ocean basins at one latitudinal belt. An assessment of both the reanalysis data and the GCM results shows that both closures exhibit very similar behaviour and are valid in the long-term mean and seasonal cycle. Interannual variability is not captured well. They become invalid for spatial scales smaller than 10^∘. Received: 30 November 1998 / Accepted: 4 July 1999     

An analytic formulation for NO and NO2 flux profiles in the atmospheric surface layer

The chemical reactivity of NO and NO₂ is so rapid that their fluxes and concentrations can be considerably modified from that expected for conserved variables in the atmospheric surface layer, even as low as a meter above the surface. Fitzjarrald and Lenschow (1983) have calculated flux and mean concentration profiles for NO, NO₂ and O₃ in the surface layer using numerical techniques. However, their solutions do not approach the photostationary state at large heights. Here we solve a simpler set of equations analytically (i.e. we assume a constant O₃ concentration and neutral hydrodynamic stability), and are able to show how the flux profiles behave at large heights assuming that the concentrations approach their photostationary values. We find, for example, that at large heights the ratio of the flux of NO to that of NO₂ is equal to the ratio of their concentrations. These results are relevant to estimating surface fluxes of NO and NO₂, and are most applicable to nonurban environments where NO and NO₂ concentrations are usually much less than O₃ concentration.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Determination of line averages of sensible heat flux using an optical method

An optical method is discussed for the determination of the line average of the vertical temperature gradient, which is subsequently used for the calculation of the scaling temperature and of the sensible heat flux. The method requires conventional surveying instrumentation only, and the measurements are carried out easily. The attainable accuracy is estimated as 10% or better. Observations are also possible during night or light-wind conditions. In a field experiment, results from this optical method were compared with measurements taken with a conventional covariance eddy-flux instrument. Excellent agreement between both methods was found.     

Comparison of Two Air Pollution Episodes over Northeast China in Winter 2016/17 Using Ground-Based Lidar

This study analyzes and compares aerosol properties and meteorological conditions during two air pollution episodes in 19–22 (E1) and 25–26 (E2) December 2016 in Northeast China. The visibility, particulate matter (PM) mass concentration, and surface meteorological observations were examined, together with the planetary boundary layer (PBL) properties and vertical profiles of aerosol extinction coefficient and volume depolarization ratio that were measured by a ground-based lidar in Shenyang of Liaoning Province, China during December 2016–January 2017. Results suggest that the low PBL height led to poor pollution dilution in E1, while the high PBL accompanied by low visibility in E2 might have been due to cross-regional and vertical air transmission. The PM mass concentration decreased as the PBL height increased in E1 while these two variables were positively correlated in E2. The enhanced winds in E2 diffused the pollutants and contributed largely to the aerosol transport. Strong temperature inversion in E1 resulted in increased PM_2.5 and PM₁₀ concentrations, and the winds in E2 favoured the southwesterly transport of aerosols from the North China Plain into the region surrounding Shenyang. The large extinction coefficient was partially attributed to the local pollution under the low PBL with high ground-surface PM mass concentrations in E1, whereas the cross-regional transport of aerosols within a high PBL and the low PM mass concentration near the ground in E2 were associated with severe aerosol extinction at high altitudes. These results may facilitate better understanding of the vertical distribution of aerosol properties during winter pollution events in Northeast China.     

Moving Beyond Monin–Obukhov Similarity Theory in Modelling Wind-Speed Profiles in the Lower Atmospheric Boundary Layer under Stable Stratification

Monin–Obukhov similarity theory (MOST) is commonly used to model the wind-speed profile at altitudes relevant to wind-power production (e.g. 10–200 m). Though reasonably accurate for unstable to weakly stable stratification, this approach becomes less accurate under increasingly stable stratification, largely due to the constant-flux surface layer assumed by MOST becoming shallower than the altitude range of interest. Furthermore, above the surface layer, the Coriolis force has a considerable influence on the wind-speed profile (in particular in the formation of low-level jets) that cannot be modelled using similarity theory. Our goal is to compare the accuracy of alternative extrapolation models that are more physically appropriate above the surface layer. Using data from the 213-m Cabauw meteorological tower in the Netherlands between July 2007 and June 2008, it is shown that MOST is accurate only at low altitudes and low stability, and breaks down at high altitudes and high stability. Local similarity is generally more accurate than MOST across all altitudes and stabilities, though the model requires turbulent flux data at multiple altitudes that is generally impractical. In contrast, a two-layer MOST–Ekman model is found to be comparable to the other models at low stability ranges and considerably more accurate in the high stability range, while requiring only a measure of surface stability and the geostrophic wind.     

Daily sensible heat flux estimation from a single measurement of surface temperature and maximum air temperature

Surface energy budget investigations of a range of agricultural surfaces in France and the African Sahel demonstrate consistent linear relationships between daily totals of sensible heat flux (H _d) and the difference between a once-a-day radiative measurement of surface temperature and the maximum air temperature at a height of 2 m. Surface temperature was measured with nadir-viewing radiothermometers near 1400 h (LST). The average residual standard error in the estimate ofH _d was 0.6 mm of equivalent evaporation.An equation for the daily sensible heat flux (H _d) having a form analogous to Dalton's evaporation formula was derived from surface energy budget considerations. This equation discriminates well between relatively homogeneous, low-cover surfaces where surface exchange characteristics can be assumed to be simple fractions of the height of the roughness elements. By contrast, data from two other crops with discontinuous plant cover suggest a much reduced sensitivity to canopy architecture. This result is not unreasonable if scalar transport were controlled by the thermal conductivity of a layer of still air close to ground level which is sheltered by the plant canopy. There is scope for further experimental and theoretical work on this matter.     

Airborne measurements of turbulent trace gas fluxes and analysis of eddy structure in the convective boundary layer over complex terrain

Using the new high-frequency measurement equipment of the research aircraft DO 128, which is described in detail, turbulent vertical fluxes of ozone and nitric oxide have been calculated from data sampled during the ESCOMPTE program in the south of France. Based on airborne turbulence measurements, radiosonde data and surface energy balance measurements, the convective boundary layer (CBL) is examined under two different aspects. The analysis covers boundary-layer convection with respect to (i) the control of CBL depth by surface heating and synoptic scale influences, and (ii) the structure of convective plumes and their vertical transport of ozone and nitric oxides. The orographic structure of the terrain causes significant differences between planetary boundary layer (PBL) heights, which are found to exceed those of terrain height variations on average. A comparison of boundary-layer flux profiles as well as mean quantities over flat and complex terrain and also under different pollution situations and weather conditions shows relationships between vertical gradients and corresponding turbulent fluxes. Generally, NO_x transports are directed upward independent of the terrain, since primary emission sources are located near the ground. For ozone, negative fluxes are common in the lower CBL in accordance with the deposition of O₃ at the surface.The detailed structure of thermals, which largely carry out vertical transports in the boundary layer, are examined with a conditional sampling technique. Updrafts mostly contain warm, moist and NO_x loaded air, while the ozone transport by thermals alternates with the background ozone gradient. Evidence for handover processes of trace gases to the free atmosphere can be found in the case of existing gradients across the boundary-layer top. An analysis of the size of eddies suggests the possibility of some influence of the heterogeneous terrain in mountainous area on the length scales of eddies.     

H2释放扰动电离层的数值模拟

由于电离层中分子性离子与电子的复合要比氧离子与电子的复合快得多,因此H₂在电离层高度释放可有效地引起电子的消耗。本文基于一个包括中性气体扩散方程和离子化学反应方程的二维动力学模型,对H₂在电离层高度释放过程进行了数值模拟研究,并分析了不同释放条件下的电子扰动特性。结果表明:1)500 mol H₂释放后,迅速向周围空间扩散,释放中心处的电子密度30 s内下降了近4%,F₂层临界频率下降了1%左右;2)在不同高度处释放H₂时,最大的电子密度相对变化率并不是在峰值高度附近处释放时出现的;3)释放化学物质的量越多,电子密度的扰动幅度也越大,但两者之间并不存在线性关系;4)相同量的H₂在电离层峰值高度处释放,白天的电子密度扰动幅度要大于夜间的扰动幅度。     

A complete dynamical ozone budget measured in the tropical marine boundary layer during PASE

The Pacific Atmospheric Sulfur Experiment (PASE) was a field mission that took place aboard the NCAR C-130 airborne laboratory over the equatorial Pacific Ocean near Christmas Island (Kirimati, Republic of Kiribati) during August?CSeptember, 2007. Eddy covariance measurements of the ozone fluxes at various altitudes above the ocean surface, along with simultaneous mapping of the horizontal gradients provided a unique opportunity to observe all of the dynamical components of the ozone budget in this remote marine environment. The results of six daytime and two sunrise flights indicate that vertical transport into the marine boundary layer from above and horizontal advection by the tradewinds are both important source terms, while photochemical destruction consisting of 82% photolysis (leading to OH production), 11% reaction with HO₂, and 7% reaction with OH provides the main sink. The overall photochemical lifetime of ozone in the marine boundary layer was found to be 6.5 days. Ocean uptake of ozone was observed to be fairly slow (mean deposition velocity of 0.024?±?0.014 cm s^?1) accounting for a diurnally averaged loss rate that was ??30% as large as the net photochemical destruction. From the measurement of deposition velocity an ozone reactivity of ??50 s^?1 in seawater is inferred. Due to the unprecedented measurement accuracy of the dynamical budget terms, unobserved photochemistry was able to be deduced, leading to the conclusion that 3.9?±?3.0 ppt (parts per trillion by volume) of NO is present on average in the daytime tropical marine boundary layer, broadly consistent with several previous studies in similar environments. It is estimated, however, that each ppt of BrO hypothetically present would counter each ppt of NO above the requisite 3.9 ppt needed for budget closure. The long-term budget of ozone is further analyzed in the buffer layer, between the boundary layer and free troposphere, and used to derive an entrainment velocity across the trade wind inversion of 0.51 ± 0.38 cm s^?1.     

Fluxes of latent heat over the oceans: climatological studies and application of satellite observations

A short review is given of the different methods by which latent heat fluxes (or evaporation) over oceans are determined. In more detail, the applicability of the bulk aerodynamical formula is discussed. This formula is mainly used for climatological studies of heat fluxes and for the application of satellite data. As an example, for climatological studies we selected the work of Isemer and Hasse, who did a re-processing of the so-called Bunker data set to determine heat fluxes over the North Atlantic Ocean. In order to check their results, Isemer and Hasse calculated the annual mean heat budget for each latitudinal delt and derived from it the required ocean heat transport. With the aid of inverse modelling, the derived ocean transport was compared with the observed ocean transport and some of the used coefficients (e.g. bulk coefficients for latent and sensible heat flux) were altered. Though the ocean heat transport is changed by a large amount (at the equator 0.3 PW, original Bunker data; 0 PW, Isemer and Hasse; 0.76 PW, after inverse modelling (all northwards)) the overall patterns of the fields of the energy fluxes remain almost unchanged. The bulk coefficient for latent heat flux for example is altered by 5.6%.The geophysical parameters necessary for the bulk aerodynamic method can be determined from satellite observations: SST, q₀, u₀. Studies are described which used data from a microwave radiometer on SEASAT and NIMBUS7 to determine latent heat flux. An error calculation shows that the obtained accuracy is between 26 and 35 W m⁻². This accuracy is adequate enough to allow reasonable estimates to be made of these fluxes. More satellites are planned for launch with microwave radiometers and scatterometers which will increase the possibility of determining geophysical parameters more accurately for use in the bulk aerodynamic formula. They will provide the database from which large-scale fieldsof latent heat flux (for time scales shorter than a month or even for actual situations) can be derived.     

Coherent Structures near the Surface in a Strongly Sheared Convective Boundary Layer Generated by Large-Eddy Simulation

Well-developed low speed and high temperature streaks in association with the alignment of convection cells are observed in a large-eddy-simulation (LES) generated strongly sheared convective boundary-layer flow, which is driven by a geostrophic wind speed of 15 m s^-1 and a surface kinematic heat flux of 0.05 K m s^-1. Vortices that drive streaky structures are identified through an eigenvalue method (lambda;₂method) near the surface. These vortices are highly elongated along the quasi-streamwise direction alternating sign of the x-component of vorticity (_x). By conditional sampling of fully developed vortices, a statistically significant coherent structure is educed. The educed vortex is elongated to the streamwise direction with the elevation angle of about 17° above the horizontal surface. However, the horizontal tilting is not clearly demonstrated in the present simulation. Fluctuation fields in the domain of the educed vortex show the existence of a low speed and high temperature streak as a direct consequence of momentum and heat transport by vortical motions. The strong ejection(upward transport of low momentum or high temperature)occurring at the higher level than that of the strong sweep (downward transport of high momentum and low temperature) can be explained from the spatial distribution of the fluctuationfields of velocity and temperature. The contribution of ejection to the Reynolds stress at z/h₁ = 0.18 is about 75%, which is slightly greater than that (70% at z/h₁ = 0.173) for the neutrally stratified atmospheric boundary layer. Ejection is also found to be dominant for the turbulent heat flux.     

Trends in middle- and upper-level tropospheric humidity from NCEP reanalysis data

The National Centers for Environmental Prediction (NCEP) reanalysis data on tropospheric humidity are examined for the period 1973 to 2007. It is accepted that radiosonde-derived humidity data must be treated with great caution, particularly at altitudes above the 500 hPa pressure level. With that caveat, the face-value 35-year trend in zonal-average annual-average specific humidity q is significantly negative at all altitudes above 850 hPa (roughly the top of the convective boundary layer) in the tropics and southern midlatitudes and at altitudes above 600 hPa in the northern midlatitudes. It is significantly positive below 850 hPa in all three zones, as might be expected in a mixed layer with rising temperatures over a moist surface. The results are qualitatively consistent with trends in NCEP atmospheric temperatures (which must also be treated with great caution) that show an increase in the stability of the convective boundary layer as the global temperature has risen over the period. The upper-level negative trends in q are inconsistent with climate-model calculations and are largely (but not completely) inconsistent with satellite data. Water vapor feedback in climate models is positive mainly because of their roughly constant relative humidity (i.e., increasing q) in the mid-to-upper troposphere as the planet warms. Negative trends in q as found in the NCEP data would imply that long-term water vapor feedback is negative—that it would reduce rather than amplify the response of the climate system to external forcing such as that from increasing atmospheric CO₂. In this context, it is important to establish what (if any) aspects of the observed trends survive detailed examination of the impact of past changes of radiosonde instrumentation and protocol within the various international networks.     

The turbulent heat flux from arctic leads

The turbulent transfer of heat from Arctic leads in winter is one of the largest terms in the Arctic heat budget. Results from the AIDJEX Lead Experiment (ALEX) suggest that the sensible component of this turbulent heat flux can be predicted from bulk quantities. Both the exponential relation N = 0.14R _x ^0.72 and the linear relation N = 1.6 × 10^–3 R _x+ 1400 fit our data well. In these, N is the Nusselt number formed with the integrated surface heat flux, and R _x is the Reynolds number based on fetch across the lead. Because of the similarity between heat and moisture transfer, these equations also predict the latent heat flux. Over leads in winter, the sensible heat flux is two to four times larger than the latent heat flux.The internal boundary layer (IBL) that develops when cold air encounters the relatively warm lead is most evident in the modified downwind temperature profiles. The height of this boundary layer, , depends on the fetch, x, on the surface roughness of the lead, z ₀ and on both downwind and upwind stability. A tentative, empirical model for boundary layer growth is % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWcaaqaaiabes% 7aKbqaaiaadQhadaWgaaWcbaGaaGimaaqabaaaaOGaeyypa0JaeqOS% di2aaeWaaeaacqGHsisldaWcaaqaaiaadQhadaWgaaWcbaGaaGimaa% qabaaakeaacaWGmbaaaaGaayjkaiaawMcaamaaCaaaleqabaGaaGim% aiaac6cacaaI4aaaaOWaaeWaaeaadaWcaaqaaiaadIhaaeaacaWG6b% WaaSbaaSqaaiaaicdaaeqaaaaaaOGaayjkaiaawMcaamaaCaaaleqa% baGaaGimaiaac6cacaaI0aaaaaaa!472D!\[\frac{\delta }{{z_0 }} = \beta \left( { - \frac{{z_0 }}{L}} \right)^{0.8} \left( {\frac{x}{{z_0 }}} \right)^{0.4} \] where L is the Obukhov length based on the values of the momentum and sensible heat fluxes at the surface of the lead, and is a constant reflecting upwind stability.Velocity profiles over leads are also affected by the surface nonhomogeneity. Besides being warmer than the upwind ice, the surface of the lead is usually somewhat rougher. The velocity profiles therefore tend to decelerate near the surface, accelerate in the mid-region of the IBL because of the intense mixing driven by the upward heat flux, and rejoin the upwind profiles above the boundary layer. The profiles thus have distinctly different shapes for stable and unstable upwind conditions.     

Effect of surface layer thickness on simultaneous transport of heat and water in a bare soil and its implications for land surface schemes

Abstract

A physically‐based multi‐layer numerical model is developed to determine the coupled transport of heat and water in the soil and in the soil‐atmosphere boundary layer. Using inputs of standard weather data and initial soil conditions the model is capable of predicting the surface energy balance components as well as water content and temperature profiles in the soil. It is used to predict these variables for a bare silt loam soil under two tillage treatments, viz. culti‐packed and left loose after disc‐harrowing, and the predicted results are compared with measurements. Very good agreement between the model predictions and measured evaporation and heat fluxes and soil water and temperatures for a ten‐day period shows that the model is capable of simulating the coupled transport of soil heat and soil water and their transfer across the soil surface‐atmosphere interface adequately.

Model predictions were compared with those of CLASS (Canadian Land Surface Scheme). It is shown that CLASS, version 2.6, provides good estimates of evaporation and hence the latent heat flux density, Q_E, under wetter soil conditions, but overestimates Q_E at moderately wet soil conditions and underestimates it under dry soil conditions. Under dry to moderately wet soil conditions the calculation of evaporation from bare soil is very sensitive to the thickness of the top layer particularly as the thickness approaches 10 cm.     

Subfilter-scale Fluxes over a Surface Roughness Transition. Part I: Measured Fluxes and Energy Transfer Rates

A wind-tunnel experiment was designed and carried out to study the effect of a surface roughness transition on subfilter-scale (SFS) physics in a turbulent boundary layer. Specifically, subfilter-scale stresses are evaluated that require parameterizations and are key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. The surface transition considered in this study consists of a sharp change from a rough, wire-mesh covered surface to a smooth surface. The resulting magnitude jump in aerodynamic roughnesses, M = ln(z ₀₁/z ₀₂), where z ₀₁ and z ₀₂ are the upwind and downwind aerodynamic surface roughnesses respectively, is similar to that of past experimental studies in the atmospheric boundary layer. The two-dimensional velocity fields used in this study are measured using particle image velocimetry and are acquired at several positions downwind of the roughness transition as well as over a homogeneous smooth surface. Results show that the SFS stress, resolved strain rate and SFS transfer rate of resolved kinetic energy are dependent on the position within the boundary layer relative to the surface roughness transition. A mismatch is found in the downwind trend of the SFS stress and resolved strain rate with distance from the transition. This difference of behaviour may not be captured by some eddy-viscosity type models that parameterize the SFS stress tensor as proportional to the resolved strain rate tensor. These results can be used as a benchmark to test the ability of existing and new SFS models to capture the spatial variability SFS physics associated with surface roughness heterogeneities.     

Air transport cruise altitude restrictions to minimize contrail formation

Abstract

While quantification of the effects of NO _x and water vapor is still at an early stage there is evidence that contrail formation could make a significant contribution to global warming. This paper builds on previous research that analyzed a policy of restricting air transport cruise altitudes to eliminate contrail formation. Our previous work [Transport. Res. D 7(6) (2002) 451], examined altitude restrictions in European airspace and concluded that this could be a beneficial policy for reducing climate change impacts from aviation. Since most of the flights in European airspace are short-haul flights, this paper evaluates the trade-offs between altitude restrictions, fuel burn and journey times for longer haul flights of up to 6000 nm. Our focus is on the North Atlantic and US airspace and we examine potential contrail fraction to determine optimal cruise altitudes for reducing contrail formation. Changes in fuel burn and travel times associated with flight levels of 18,000 and 31,000 ft for different aircraft types are analyzed. We find that, in most cases, CO₂ emission increases would be unlikely to entirely counteract the benefit of possible reductions in contrail formation. For some aircraft types, the percentage increase in emitted CO₂ was found to be strongly dependent on journey length. In general, journey times appear not to be a major issue except for some aircraft types. Our results suggest that reducing aircraft cruise altitudes could be a beneficial policy for mitigating climate change impacts from the aviation sector. This is clearly dependent on aircraft type and the distances traveled, but more importantly on ambient atmospheric conditions which can vary significantly between regions and due to daily variation. This suggests that real time flight planning to minimize contrail formation should be investigated as a possible climate mitigation policy.     

Molecular dissipation of turbulent fluctuations in the convective mixed layer part II: Height variations of characteristic time scales and experimental test of molecular dissipation models

Extensive turbulence measurements from the Limagne and Beauce experiments were used to compute a characteristic time scale of the turbulence field (Τ = second moment/dissipation rate) for turbulent kinetic energy, temperature and humidity variances, and temperature-humidity covariance. The height variations of these time scales were analysed. The characteristic half-time scale Τ/2 of the turbulent velocity field was found, as expected, to be of the same order of magnitude as the large-eddy time scale Τ _L = Z_i/w*, showing that the turbulence structure is controlled by large eddies in the bulk of the mixed layer. The increase of Τ/2 above z/Z _i ～- 0.7 implies, however, that this time scale is no longer relevant to destruction of turbulent kinetic energy in the statically stable region with negative heat fluxes. An effective time scale Τ_eff, introduced by Zeman (1975), has been computed and its behaviour discussed. The scales for θ′ ², q′², and θ′q′ were found to be much shorter than Τ. Furthermore, a significant difference in behaviour was also revealed between the characteristic time scales of temperature and humidity fields in the stable layer. By using these experimental estimates, we tested some of the models for molecular dissipations, which are currently in use in higher order closure atmospheric boundary-layer models. The parameterized dissipation rates for θ′ ², and q′ ² agree well qualitatively with experimental estimates in the bulk of the mixed layer. In the stable layer, however, the parameterized dissipation rate ε _θ tends to become larger than the experimental ones although the parameterized dissipation rate ε _q still agrees with the experimental ones. For the molecular dissipation of θ′q′, this current model becomes physically inconsistent in the middle part of the mixed layer, because this term may become a production term for temperature-humidity covariance.     

A global three-dimensional model of the tropospheric sulfur cycle

The tropospheric part of the atmospheric sulfur cycle has been simulated in a global three-dimensional model. The model treats the emission, transport, chemistry, and removal processes for three sulfur components; DMS (dimethyl sulfide), SO₂ and SO₄ ^2– (sulfate). These processes are resolved using an Eulerian transport model, the MOGUNTIA model, with a horizontal resolution of 10° longitude by 10° latitude and with 10 layers in the vertical between the surface and 100 hPa. Advection takes place by climatological monthly mean winds. Transport processes occurring on smaller space and time scales are parameterized as eddy diffusion except for transport in deep convective clouds which is treated separately. The simulations are broadly consistent with observations of concentrations in air and precipitation in and over polluted regions in Europe and North America. Oxidation of DMS by OH radicals together with a global emission of 16 Tg DMS-S yr^–1 from the oceans result in DMS concentrations consistent with observations in the marine boundary layer. The average turn-over times were estimated to be 3, 1.2–1.8, and 3.2–6.1 days for DMS, SO₂, and SO₄ ^2– respectively.     

A sensitive fast-response probe to measure turbulence and heat flux from any airplane

The theory, configuration, and accuracy of an inexpensive probe to measure turbulence from a small airplane are presented. The probe employs a nine-hole pressure-sphere design along with inprobe high-frequency pressure, temperature, and acceleration sensors. This sensor suite is specifically designed to extend mass, momentum and energy eddy-flux measurement to the higher frequencies characteristic of marine and nocturnal boundary layers. The probe is part of a mobile flux system, independent of the conveyance, which does not require a separate Inertial Navigation System.The new nine-port pressure sphere turbulence probe allows accurate turbulent velocity measurement with proper probe installation and appropriate computation technique for dynamic pressure. A thermistor in the central pressure port provides simultaneous temperature measurement, at a location symmetrical with respect to the flow, for accurate determination of true airspeed and heat flux. The probemounted temperature sensor gives heat fluxes with variance 5% of the mean in a weakly-turbulent marine boundary layer.Oak Ridge Associated Universities, assigned to NOAA/ATDD.     

Dry-air boundary conditions for correction of eddy flux measurements

For measurements of eddy fluxes in the atmospheric boundary layer of gases (such as CO₂) whose average concentration is very large compared to the fluctuations, corrections for air density fluctuations are required. With the boundary condition of no flux of dry air at the surface, the evaporation correction to eddy fluxes is 2.6 times larger than has been estimated with the boundary condition of no mass flux at all at the surface. The heat flux correction is also increased by a few per cent.