Bulk characteristics of heat transfer in the unstable,baroclinic atmospheric boundary layer

Field data for the unstable, baroclinic, atmospheric boundary layer over land and over the sea are considered in the context of a general similarity theory of vertical heat transfer. The dependence of δθ/θ_* upon logarithmic functions of h _c z _T and stability (through the similarity function C) is clearly demonstrated in the data. The combined data support the conventional formulation for the heat transfer coefficient δθ/θ_* when,

1. the surface scaling length is z _T (« z ₀), the height at which the surface temperature over land is obtained by extrapolation of the temperature profile
2. the height scale is taken as the depth of convective mixing h _c
3. the temperature profile equivalent of the von Karman constant is taken as 0.41
4. areal average, rather than single point, values of δθ are employed in strongly baroclinic conditions. No significant effect of baroclinity or the height scale ratio as proposed in the general theory is found. Variations in C about a linear regression relation against stability are most probably due to uncertainties in the areal surface temperature and to experimental errors in general temperature measurements.

     

Experiments on scalar dispersion within a model plant canopy part II: An elevated plane source

This paper describes a wind-tunnel experiment on the dispersion of trace heat from an effectively planar source within a model plant canopy, the source height being h _s = 0.80 h _c , where h _c is the canopy height. A sensor assembly consisting of three coplanar hot wires and one cold wire was used to make simultaneous measurements of the temperature and the streamwise and vertical velocity components. It was found that:

1. The thermal layer consisted of two parts with different length scales, an inner sublayer (scaling with h _s and h _c ) which quickly reached streamwise equilibrium downstream of the leading edge of the source, and an outer sublayer which was self-preserving with a length scale proportional to the depth of the thermal layer.
2. Below 2h _c , the vertical eddy diffusivity for heat from the plane source (K _HP ) was substantially less than the far-field limit of the corresponding diffusivity for heat from a lateral line source at the same height as the plane source. This shows that dispersion from plane or other distributed sources in canopies is influenced, near the canopy, by turbulence ‘memory’ and must be considered as a superposition of both near-field and far-field processes. Hence, one-dimensional models for scalar transport from distributed sources in canopies are wrong in principle, irrespective of the order of closure.
3. In the budgets for temperature variance, and for the vertical and streamwise components of the turbulent heat flux, turbulent transport was a major loss between h _s and h _c and a principal gain mechanism below h _s , as also observed in the budgets for turbulent energy and shear stress.
4. Quadrant analysis of the vertical heat flux showed that sweeps and ejections contributed about equal amounts to the heat flux between h _s and h _c , though among the more intense events, sweeps were dominant. Below h _s , almost all the heat was transported by sweeps.

     

Experiments on scalar dispersion within a model plant canopy,part III: An elevated line source

An experiment is reported in which heat was released as a passive tracer from an elevated lateral line source within a model plant canopy, with h _s = 0.85 h _c (h _s and h _c being the source and canopy heights, respectively). A sensor assembly consisting of three coplanar hot wires and one cold wire was used to measure profiles of mean temperature % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikamaana% aabaGaeqiUdehaaiaacMcaaaa!390C!\[(\overline \theta )\], temperature variance (Σ_θ ²), vertical and streamwise turbulent heat fluxes, and third moments of wind and temperature fluctuations. Conclusions were:

1. Despite the very heterogeneous flow within the canopy, the observed dispersive heat flux (due to spatial correlation between time-averaged temperature and vertical velocity) was small. However, there is evidence from the plume centroid (which was lower than h _s at the source) of systematic recirculating motions within the canopy.
2. The ratio % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS% baaSqaaiabeI7aXjaab2gacaqGHbGaaeiEaaqabaGccaGGVaWaa0aa% aeaacqaH4oqCaaWaaSbaaSqaaiaab2gacaqGHbGaaeiEaaqabaaaaa!41DF!\[\sigma _{\theta {\text{max}}} /\overline \theta _{{\text{max}}} \] (of maximum values on vertical profiles) decreased from 1 near the source to an asymptotic value of 0.4 far downstream, in good agreement with previous experimental and theoretical work for concentration fluctuations in the surface layer well above the canopy.
3. The eddy diffusivity for heat from the line source (K _HL ) increased, downstream of the source, to a nearly constant ‘far-field’ vertical profile. Within the canopy, the far-field K _HL was an order of magnitude larger than K _HP , the equivalent diffusivity for a plane source; well above the canopy, the two were equal. The time scale defined by (far-field K _HL )/(vertical velocity variance) was independent of height within the canopy.
4. Budgets for temperature variance, vertical heat flux and streamwise heat flux are remarkably similar to the equivalent budgets for an elevated line source in the surface layer well above the canopy, except in the lower part of the canopy in the far field, where vertical transport is much more important than in the surface layer.
5. A random flight simulation of the mean height and depth of the temperature plume was generally in good agreement with experiment. However, details of the temperature and streamwise turbulent heat flux profiles were not correct, suggesting that the model formulation needs to be improved.

     

Wind profiles at the Boulder Tower

Analysis of wind profiles at the Boulder Tower (BAO) leads to these conclusions:

1. The variation of roughness with wind direction found earlier is confirmed. Roughness lengths measured on the tower are larger than those measured close to the surface.
2. The profiles and measurements of Reynolds stress are consistent with a von-Karman constant of 0.35.
3. The form φ_m=(1?15z/L)^-1/3 fits best in the range -0.6 < z/L < 0. In the range 0 < z/L < 0.5, θ _m ～ 1 + 4.7z/L provides a good fit to the observations. For z/L < 0.1, φ _m also depends on h, the thickness of the PBL. For z/L < -0.6, Φ _m approaches the constant 0.5, in contrast to all previous suggestions. For larger stabilities, the upper level is usually not in the surface layer, and wind ratios become independent of z/L.
4. With snow cover, the effective roughness diminishes to about 1 cm, even for directions for which the roughness length without snow is large.
5. Estimation of winds at 100 or 150 m from information near the surface is best for similarity theory provided that the ratio of height to Monin-Obukhov L is less than 0.1. For larger z/L, simple power laws seem more appropriate.

     

Atmospheric turbulence within and above a Douglas-fir stand. Part I: Statistical properties of the velocity field

This is the first of two papers reporting the results of a study of the turbulence regimes and exchange processes within and above an extensive Douglas-fir stand. The experiment was conducted on Vancouver Island during a two-week rainless period in July and August 1990. The experimental site was located on a 5^o slope. The stand, which was planted in 1962, and thinned and pruned uniformly in 1988, had a (projected) leaf area index of 5.4 and a heighth=16.7 m. Two eddy correlation units were operated in the daytime to measure the fluctuations in the three velocity components, air temperature and water vapour density, with one mounted permanently at a height of 23.0m (z/h=1.38) and the other at various heights in the stand with two to three 8-hour periods of measurement at each level. Humidity and radiation regimes both above and beneath the overstory and profiles of wind speed and air temperature were also measured. The most important findings are:

1. A marked secondary maximum in the wind speed profile occurred in the middle of the trunk space (aroundz/h=0.12). The turbulence intensities for the longitudinal and lateral velocity components increased with decreasing height, but the intensity for the vertical velocity component had a maximum atz/h=0.60 (middle of the canopy layer). Magnitudes of the higher order moments (skewness and kurtosis) for the three velocity components were higher in the canopy layer than in the trunk space and above the stand.
2. There was a 20% reduction in Reynolds stress fromz/h=1.00 to 1.38. Negative Reynolds stress or upward momentum flux perisistently occurred atz/h=0.12 and 0.42 (base of the canopy), and was correlated with negative wind speed gradients at the two heights. The longitudinal pressure gradient due to the land-sea/upslope-downslope circulations was believed to be the main factor responsible for the negative Reynolds stress.
3. Momentum transfer was highly intermittent. Sweep and ejection events dominated the transfer atz/h=0.60, 1.00 and 1.38, with sweeps playing the more important role of the two atz/h=0.60 and 1.00 and the less important role atz/h=1.38. But interaction events were of greater magnitude than sweep and ejection events atz/h=0.12 and 0.42.

     

Multiple windbreaks: An aeolean ensemble

Near-neutral measurements of the turbulent wind field within and above a sequence of 15 parallel windbreaks on a flat pastoral site are presented. The windbreak fences each had a porosity of 60% and were equally-spaced at 6 times their height (h = 2 m). The following conclusions seem justified for wind directions within 10 ° of the normal to the array:

1. Above the windbreaks (2h), mean windspeeds first decreased and then increased asymptotically to a value in equilibrium with the new surface roughness. At 0.5h, windspeeds exhibited a slow increase down the entire array.
2. Reflecting differences in approach flows, the drag on the initial fence was almost twice that on barriers farther downstream. This reduction in momentum extraction per windbreak was associated with an elevation in the zero-plane displacement to a level equal to 0.8h.
3. At positions well-removed from the initial fences, mean windspeeds were reduced throughout the entire region below shelter height. In this region, the flow became increasingly dominated by downward moving air with velocities much greater than the local average. The zone of reduced turbulence was small, extending only 2h downstream of a barrier at a height of 0.25h. This corresponded with the region excluded from smoke trails released at the top of windbreaks.
4. An approximate TKE budget mid-way between windbreaks 7 and 8 suggests that shear and wake production peak near z = h and that production is balanced by dissipation and vertical transport components. Advective and inertial interaction terms are negligible at this midway position but are likely to be major sources of TKE closer to the windbreak. Local equilibrium is attained above z = 1.5h implying the existence of a constant-stress layer.

The measurements show the practical difficulty of simultaneously reducing both mean windspeeds and turbulence levels with repeated windbreaks at conventional spacings for horticultural applications.     

A comparison of ABL heights inferred routinely from lidar and radiosondes at noontime

The height of the atmospheric boundary layer (ABL) obtained with lidar and radiosondes is compared for a data set of 43 noon (12.00 GMT) cases in 1984. The data were selected to represent the synoptic circulation types appropriately. Lidar vertical profiles at 1064 nm were used to obtain three estimates for the ABL height (h _lid), based on the first gradient in the back-scatter profile, namely, at the beginning, middle and top of the gradient. The boundary-layer height obtained with the radiosondes (h _s) was determined with the dry-parcel-intersection method in unstable conditions. As a first guess for near-neutral and stable conditions, the height of the first significant level in the potential temperature profile was taken. Overall, the boundary-layer thickness estimates agree surprisingly well (regression lineh _lidb=h_s:cc.=0.93 and the standard error=121 m). However, in 10% of the cases, the lidar estimate was significantly lower (difference>400 m) than the routinely inferredh _s. These outliers are discussed separately. For stable conditions, an estimate of ABL height (h _N) is also made based on the friction velocity and the Brunt-Väisälä frequency. The agreement betweenh _Nandh _lidbis good. Discrepancies between the two methods are caused by:

1. rapid growth of the boundary layer arround the measurement time;
2. the presence of a deep entrainment layer leading to a large zone in which quantities are not well mixed;
3. a large systematic error of 100–200 m in the estimate of boundary-layer height obtained from the radiosonde due to the way that profiles are recorded, as a series of significant points.

     

Experiments on scalar dispersion within a model plant canopy part I: The turbulence structure

This is the first of a series of three papers describing experiments on the dispersion of trace heat from elevated line and plane sources within a model plant canopy in a wind tunnel. Here we consider the wind field and turbulence structure. The model canopy consisted of bluff elements 60 mm high and 10 mm wide in a diamond array with frontal area index 0.23; streamwise and vertical velocity components were measured with a special three-hot-wire anemometer designed for optimum performance in flows of high turbulence intensity. We found that:

1. The momentum flux due to spatial correlations between time-averaged streamwise and vertical velocity components (the dispersive flux) was negligible, at heights near and above the top of the canopy.
2. In the turbulent energy budget, turbulent transport was a major loss (of about one-third of local production) near the top of the canopy, and was the principal gain mechanism lower down. Wake production was greater than shear production throughout the canopy. Pressure transport just above the canopy, inferred by difference, appeared to be a gain in approximate balance with the turbulent transport loss.
3. In the shear stress budget, wake production was negligible. The role of turbulent transport was equivalent to that in the turbulent energy budget, though smaller.
4. Velocity spectra above and within the canopy showed the dominance of large eddies occupying much of the boundary layer and moving downstream with a height-independent convection velocity. Within the canopy, much of the vertical but relatively little of the streamwise variance occurred at frequencies characteristic of wake turbulence.
5. Quadrant analysis of the shear stress showed only a slight excess of sweeps over ejections near the top of the canopy, in contrast with previous studies. This is a result of improved measurement techniques; it suggests some reappraisal of inferences previously drawn from quadrant analysis.

     

The scavenging of nitrate by clouds and precipitation

A model with spectral microphysics was developed to describe the scavenging of nitrate aerosol particles and HNO₃ gas. This model was incorporated into the dynamic framework of an entraining air parcel model with which we computed the uptake of nitrate by cloud drops whose size distribution changes with time because of condensation, collision-coalescence and break-up. Significant differences were found between the scavenging behavior of nitrate and our former results on the scavenging behavior of sulfate. These reflect the following chemical and microphysical differences between the two systems:

1. nitrate particles occur in a larger size range than sulfate particles.
2. HNO₃ has a much greater solubility than SO₂ and is taken up irreversibly inside the drops in contrast to SO₂.
3. nitric acid in the cloud water is formed directly on uptake of HNO₃ gas whereas on uptake of SO₂ sulfuric acid is formed only after the reaction with oxidizing agents such as e.g., H₂O₂ or O₃.
4. nitrate resulting from uptake of HNO₃ is confined mainly to small drops, whereas sulfate resulting from uptake of SO₂ is most concentrated in the largest, oldest drops, which have had the greatest time for reaction.

Sensitivity studies showed that the nitrate concentration of small drops is significantly affected by the mass accommodation coefficient.     

Atmospheric boundary layer research at Cabauw

At Cabauw, The Netherlands, a 213 m high mast specifically built for meteorological research has been operational since 1973. Its site, construction, instrumentation and observation programs are reviewed. Regarding analysis of the boundary layer at Cabauw, the following subjects are discussed:

- terrain roughness;

- Monin-Obukhov theory in practice;

- the structure of stable boundary layers;

- observed evolution of fog layers;

- inversion rise and early morning entrainment;

- use of the geostrophic wind as a predictor for wind profiles;

- height variation of wind climate statistics;

- air pollution applications: long range transport and short range dispersion;

- dependence of sound wave propagation on boundary-layer structure;

- testing of weather and climate models.

     

Variability of energy fluxes in relation to the net-radiation of urban and suburban areas: a case study

This paper addresses the relation between the net-radiation (Q ^*) and the ground heat flux (Q _G), the energy stored in the soil ( $\Updelta Q_{\rm S}$ ), and the residual of the energy partition (R = Q ^* ? Q _H  ? Q _E ) of urban and suburban areas of Oklahoma City, USA. These three forms of energy were observed or estimated from observations taken during Joint Urban 2003 Campaign. The database includes net-radiation, soil temperature, ground heat flux, and turbulent fluxes. In most cases the estimates of the energy stored in the soil were obtained by assuming roughly a certain type of soil and an effective soil depth. From the residuals it seems to be possible to distinguish the urban boundary layer from the suburban boundary layer when plotted as a function of net-radiation. Hysteresis coefficients were computed for fits of net-radiation against R, $\Updelta Q_{\rm S}$ and Q _G. In particular, the hysteresis patterns show that Q ^* vs. R represents clearer urban areas or suburban areas under the influence of an urban “plume”. On the other hand, hysteresis curves obtained from $\Updelta Q_{\rm S}$ or Q _G account for better the ground composition. A possible consequence is that the land use of urban areas could be roughly inferred from curve shapes such as Q ^* vs. R, or Q ^* versus another input variable representing the storage term. The objective is to show the variability of the subsurface-related energy fluxes across an urban area using these three different quantities and also to show that $\Updelta Q_{\rm S}, \,Q_{\rm G}$ , or R (and their corresponding hysteresis curves) are likely to be quantitatively different, which have not been clearly stated in the literature.     

Bare-ground surface heat and water exchanges under dry conditions: Observations and parameterization

A simplified land-surface parameterization is tested against bare-soil data collected during the EFEDA experiment conducted in Spain in June 1991. A complete data set, made up of soil properties as well as hydrological and atmospheric measurements, is described and discussed. The 11-day data set is characterized by very dry conditions and high surface temperatures during the day. Large values of sensible and soil heat fluxes and small values of surface evaporation (≈1 mm/day) were observed. This data set was modelled, leading to the following conclusions:

1. In the model, the parameterization provides values of the soil thermal properties and subsequently of the predicted soil heat fluxes which are overestimated when compared with the observations.
2. Following the literature, a value of the ratio between the roughness lengths for momentumZ _oand heatZ _ohof close to 10 for fairly homogeneous areas of bare soil and vegetation is used. This value leads to a fair prediction of the surface temperature. If the roughness lengths were taken to be equal, as is often assumed in atmospheric modelling, a poorer prediction results.
3. Finally, the vapor phase transfer mode is found dominant close to the surface and a modified parameterization including this effect is proposed. It allows a fair prediction of both surface evaporation and near-surface water content.

     

Surface influence upon vertical profiles in the nocturnal boundary layer

Near-surface wind profiles in the nocturnal boundary layer, depth h, above relatively flat, tree-covered terrain are described in the context of the analysis of Garratt (1980) for the unstable atmospheric boundary layer. The observations at two sites imply a surface-based transition layer, of depth z _*, within which the observed non-dimensional profiles Φ _M ⁰ are a modified form of the inertial sub-layer relation \(\Phi _M \left( {{z \mathord{\left/ {\vphantom {z L}} \right. \kern-0em} L}} \right) = \left( {{{1 + 5_Z } \mathord{\left/ {\vphantom {{1 + 5_Z } L}} \right. \kern-0em} L}} \right)\) according to $$\Phi _M^{\text{0}} \simeq \left( {{{1 + 5z} \mathord{\left/ {\vphantom {{1 + 5z} L}} \right. \kern-\nulldelimiterspace} L}} \right)\exp \left[ { - 0.7\left( {{{1 - z} \mathord{\left/ {\vphantom {{1 - z} z}} \right. \kern-\nulldelimiterspace} z}_ * } \right)} \right]$$ , where z is height above the zero-plane displacement and L is the Monin-Obukhov length. At both sites the depth z _* is significantly smaller than the appropriate neutral value (z _*N ) found from the previous analysis, as might be expected in the presence of a buoyant sink for turbulent kinetic energy.     

Assessing the Consequences of Climate Change for Food and forest Resources: A View from the IPCC

Important findings on the consequences of climate change for agriculture and forestry from the recently completed Third Assessment Report (TAR) of the Intergovernmental Panel on Climate Change (IPCC) are reviewed, with emphasis on new knowledge that emerged since the Second Assessment Report (SAR). The State-Pressure-Response-Adaptation model is used to organize the review. The major findings are:

• Constant or declining food prices are expected for at least the next 25 yr, although food security problems will persist in many developing countries as those countries deal with population increases, political crisis, poor resource endowments, and steady environmental degradation. Most economic model projections suggest that low relative food prices will extend beyond the next 25 yr, although our confidence in these projections erodes farther out into the 21st century.
• Although deforestation rates may have decreased since the early 1990s, degradation with a loss of forest productivity and biomass has occurred at large spatial scales as a result of fragmentation, non-sustainable practices and infrastructure development.
• According to United Nations estimates, approximately 23% of all forest and agricultural lands were classified as degraded over the period since World War II.
• At a worldwide scale, global change pressures (climate change, land-use practices and changes in atmospheric chemistry) are increasingly affecting the supply of goods and services from forests.
• The most realistic experiments to date – free air experiments in an irrigated environment – indicate that C₃ agricultural crops in particular respond favorably to gradually increasing atmospheric CO₂ concentrations (e.g., wheat yield increases by an average of 28%), although extrapolation of experimental results to real world production where several factors (e.g., nutrients, temperature, precipitation, and others) are likely to be limiting at one time or another remains problematic. Moreover, little is known of crop response to elevated CO₂ in the tropics, as most of the research has been conducted in the mid-latitudes.
• Research suggests that for some crops, for example rice, CO₂ benefits may decline quickly as temperatures warm beyond optimum photosynthetic levels. However, crop plant growth may benefit relatively more from CO₂ enrichment in drought conditions than in wet conditions.
• The unambiguous separation of the relative influences of elevated ambient CO₂ levels, climate change responses, and direct human influences (such as present and historical land-use change) on trees at the global and regional scales is still problematic. In some regions such as the temperate and boreal forests, climate change impacts, direct human interventions (including nitrogen-bearing pollution), and the legacy of past human activities (land-use change) appear to be more significant than CO₂ fertilization effects. This subject is, however an area of continuing scientific debate, although there does appear to be consensus that any CO₂ fertilization effect will saturate (disappear) in the coming century.
• Modeling studies suggest that any warming above current temperatures will diminish crop yields in the tropics while up to 2–3 ^°C of warming in the mid-latitudes may be tolerated by crops, especially if accompanied by increasing precipitation. The preponderance of developing countries lies in or near the tropics; this finding does not bode well for food production in those countries.
• Where direct human pressures do not mask them, there is increasing evidence of the impacts of climate change on forests associated with changes in natural disturbance regimes, growing season length, and local climatic extremes.
• Recent advances in modeling of vegetation response suggest that transient effects associated with dynamically responding ecosystems to climate change will increasingly dominate over the next century and that during these changes the global forest resource is likely to be adversely affected.
• The ability of livestock producers to adapt their herds to the physiological stress of climate change appears encouraging due to a variety of techniques for dealing with climate stress, but this issue is not well constrained, in part because of the general lack of experimentation and simulations of livestock adaptation to climate change.
• Crop and livestock farmers who have sufficient access to capital and technologies should be able to adapt their farming systems to climate change. Substantial changes in their mix of crops and livestock production may be necessary, however, as considerable costs could be involved in this process because investments in learning and gaining experience with different crops or irrigation.
• Impacts of climate change on agriculture after adaptation are estimated to result in small percentage changes in overall global income. Nations with large resource endowments (i.e., developed countries) will fare better in adapting to climate change than those with poor resource endowments (i.e., developing countries and countries in transition, especially in the tropics and subtropics) which will fare worse. This, in turn, could worsen income disparities between developed and developing countries.
• Although local forest ecosystems will be highly affected, with potentially significant local economic impacts, it is believed that, at regional and global scales, the global supply of timber and non-wood goods and services will adapt through changes in the global market place. However, there will be regional shifts in market share associated with changes in forest productivity with climate change: in contrast to the findings of the SAR, recent studies suggest that the changes will favor producers in developing countries, possibly at the expense of temperate and boreal suppliers.
• Global agricultural vulnerability is assessed by the anticipated effects of climate change on food prices. Based on the accumulated evidence of modeling studies, a global temperature rise of greater than 2.5 ^°C is likely to reverse the trend of falling real food prices. This would greatly stress food security in many developing countries.

     

Limitations of an eddy-correlation technique for the determination of the carbon dioxide and sensible heat fluxes

A modified infrared CO₂ gas analyzer, a small thermocouple assembly, a heated-thermocouple anemometer for horizontal wind, and a propeller-type vertical wind sensor were used to measure the eddy fluxes of heat and CO₂ above a corn crop. Experimental results of these fluxes are discussed. The main sources of errors of the eddy fluxes using these instruments were estimated:

1. Sensors with a time constant of 0.5 s appear to be fast enough to detect most of the vertical CO₂ transfer as long as the sensors are located at least one meter above the crop surface.
2. The deviation from steady-state conditions for 10-min periods was found to have a significant effect on the eddy flux estimates.
3. Temperature fluctuations of the air sample passing through the CO₂ infrared gas analyzer were found to be non-negligible but could be easily corrected.
4. A 1° misalignment of the vertical anemometer affected these eddy fluxes by less than 10% under all circumstances studied.

     

Hourly estimation of soil heat flux density at the soil surface with three models and two field methods

Heat flux density at the soil surface (G ₀) was evaluated hourly on a vegetal cover 0.08 m high, with a leaf area index of 1.07 m² m^?2, during daylight hours, using Choudhury et al. (Agric For Meteorol 39:283–297, 1987) ( $ G_0^{\text{rn}} $ ), Santanello and Friedl (J Appl Meteorol 42:851–862, 2003) ( $ G_0^{\text{s}} $ ), and force-restore ( $ G_0^{\text{fr}} $ ) models and the plate calorimetry methodology ( $ G_0^{\text{pco}} $ ), where the gradient calorimetry methodology (G _0R ) served as a reference for determining G ₀. It was found that the peak of G _0R was at 1 p.m., with values that ranged between 60 and 100 W m^?2 and that the G ₀/Rn relation varied during the day with values close to zero in the early hours of the morning and close to 0.25 in the last hours of daylight. The $ G_0^{\text{s}} $ model presented the best performance, followed by the $ G_0^{\text{rn}} $ and $ G_0^{\text{fr}} $ models. The plate calorimetry methodology showed a similar behavior to that of the gradient calorimetry referential methodology.     

Study of the mean monthly air temperatures during the successive sunspot cycles

In the first part of this paper (I–III) the mean monthly air temperatures during successive sunspot cycles are examined for six stations of Central and Northwestern Europe. Based upon the observed data the variations between two consecutive 11-year periods of solar activity are studied for the following elements:

1. a)
   The differenceT _h?T_c between the mean temperatures of the four warmer (May–August) and the four colder (November–February) months of the year.     
   
   

Aircraft measurements of oxides of nitrogen along the eastern rim of the Asian continent: Winter observations

Aircraft observations of oxides of nitrogen (NO _y ), measured with a ferrous sulfate converter, over the sea surrounding the Japanese islands (30–43° N, 131–141° E) were carried out in the winter of 1983 and 1984 at altitudes mostly between 3 and 8 km. NO _y defined here is the sum of NO, NO₂, and other unstable oxides of nitrogen that are converted to NO by ferrous sulfate. The main observations were:

1. Over the Pacific Ocean between the latitudes of 30–35° N, the observed NO _y mixing ratio between 3 and 8 km was a fairly constant 200 pptv. The NO mixing ratio increased with altitude from 15 pptv at 3 km to 35 pptv at 7 km.
2. Over the Sea of Japan, tropospheric NO _y mesured between 1 and 6 km started increasing with latitude North of 35° N and reached about 1000 pptv at 40° N.
3. NO _y was measured in an air mass transported from the stratosphere near a tropopause fold region. When the ozone mixing ratio was between 80 and 140 ppbv, the NO _y mixing ratio was about 200 pptv.

     

Increases of wet deposition at remote sites in Japan from 1991 to 2009

Temporal trends in wet deposition of major ions were explored at nationwide remote sites in Japan from April 1991 to March 2009 by using the seasonal Kendall slope estimator and the nonparametric seasonal Kendall test. For the trend analysis, datasets from eight remote sites (Rishiri, Echizenmisaki, Oki, Ogasawara, Shionomisaki, Goto, Yakushima, and Amami) were selected from the Japanese Acid Deposition Survey (JADS) conducted by the Ministry of the Environment. Deposition of H⁺ has been increasing at remote sites in Japan on a national scale. Significant (p????0.05) increases in H⁺ deposition were detected with changes of +3?C+9?%?year^?1 at seven sites, while insignificant increases were observed at one site. Depositions of non-sea salt (nss)-SO ₄ ^2? and NO ₃ ^? significantly increased at four and six sites, respectively, with changes of +1?C+3?%?year^?1. Significant increases in precipitation at four sites would have contributed to the increase in depositions of H⁺, nss-SO ₄ ^2? , and NO ₃ ^? . The emission trends of SO₂ and NO_x did not corresponded to the deposition trends of nss-SO ₄ ^2? and NO ₃ ^? . The different trends indicated that temporal variation of precipitation amount trend dominated the deposition trends.     

Measurements of internal waves and turbulence in two-dimensional stratified shear flows

A turbulent stratified shear flow is generated in a towing tank by towing a grid or a circular cylinder through a tank of stratified salt water. The internal waves and turbulence generated in these flows are visualized with shadowgraphs and measured with quartz-coated hot-film probes (up to four probes for velocity fluctuations) and single-electrode conductivity probes (up to four probes for salinity fluctuations) which are towed at the same speed as the obstacle. The velocity and salinity signals are recorded on magnetic tapes. A portion of these signals is processed directly-on-line with a digital computer. From these shadowgraphs and probe measurements, we observe that

1. Far downstream of the obstacle where the turbulence has already subsided, the stratified fluid always has a layered structure. This layered structure persists for a long time, and is a result of the convection of turbulently mixed layers by the mean flow. These results indicate that in the regions of a stably stratified atmosphere and ocean where the turbulence has subsided, one could often find layered structure.
2. There are spectral peaks and valleys in the measured velocity and salinity autospectra when the stratifications are sufficiently strong. Under certain conditions, these spectral peaks tend to lift up the spectral curves to show substantialf ^?5/3 subranges, although the turbulence Reynolds numbers are too low for the flows to have recognizable inertial subranges. This anomalousf ^?5/3 subrange demonstrates the pitfalls of using spectral measurements in thef ^?5/3 subrange to predict the turbulent energy dissipation rate through the Kolmogorov hypothesis.
3. A diagnostic method is developed for distinguishing internal waves from turbulence, utilizing their phase characteristics. The phase characteristics can be conveniently examined from the cospectra and quadrature spectra measurements of: (a), two vertically separated velocity probes; (b), two vertically separated density probes; and (c), a velocity probe and a density probe. This method is demonstrated to be useful in the laboratory and can be applied directly to atmospheric and oceanic measurements to distinguish internal waves from turbulence.
4. From the coherency measurements, it is found that the entire turbulent stratified wake is actually whipping up and down at a frequency corresponding to the Brunt-Väisälä frequency. This indicates that similar stratified shear flows in the atmosphere and in the ocean, such as the jet streams in the atmosphere and the Cromwell current in the ocean, may oscillate vertically, which in turn can induce horizontal oscillation and meandering.