Flux-profile Relationships for Wind Speed and Temperature in the Stable Atmospheric Boundary Layer

Wind and temperature profiles in the stable boundary layer were analyzed in the context of MoninObukhov similarity. The measurements were made on a 60-m tower in Kansas during October 1999 (CASES-99). Fluxprofile relationships, obtained from these measurements in their integral forms, were established for wind speed and temperature. Use of the integral forms eliminates the uncertainty and accuracy issues resulting from gradient computations. The corresponding stability functions, which were nearly the same for momentum and virtual sensible heat, were found to exhibit different features under weakly stable conditions compared to those under strongly stable conditions. The gradient stability functions were found to be linear, namely _m = 1+ 5.8 and _h = 1 + 5.4 up to a limit of the MoninObukhov stability parameter = 0.8; this is consistent with earlier findings. However, for stronger stabilities beyond a transition range, both functions were observed gradually to approach a constant, with a value of approximately 7. To link these two distinct regimes, a general but pliable functional form with only two parameters is proposed for the stability functions, covering the entire stability range from neutral to very stable conditions.     

A non-linear climate oscillator controlled by biogeochemical cycling in the ocean: an alternative model of Quaternary ice age cycles

A new, biogeochemical model of ice age cycles is developed and applied which explains major features of climate variations in the late Quaternary —rapid ice age terminations, large glacial-interglacial amplitudes and 100-kyr cycles — in a way consistent with the paleorecord. Existing models which invoke non-linear, ice-sheet-earth-crust dynamics to explain ice age cycles are not consistent with simultaneous terminations in both hemispheres and other phase relationships implied by the paleorecord. The present model relates climate change to oscillations of oceanic primary (new) production controlled by the availability of inorganic nitrogen. Large oscillations follow shelf erosion events triggered by small sea-level drops. These drops are due to glacial buildup associated with a minimum in Northern Hemisphere insolation. Rapid global warming at terminations is initiated by open ocean denitrification events leading to new production crashes and rapid modification of atmospheric trace gas concentrations (CO₂, DMS, N₂O). Other feedbacks of the land-ice-atmosphere-ocean system control the rest of the climate cycle. 100-kyr cycles derive from orbital pacemaking of the strong, low-frequency model response. Results suggest that the climate regime transition near 800 kyr B.P. may be related to changes in the continental shelf slope, that existing chronologies based on orbital tuning may need to be revised and that temporary increases in atmospheric N₂O concentrations at terminations, due to the denitrification events, may have caused significant greenhouse warming. A spike of elevated N₂O concentration at terminations may be recorded in polar ice.     

Orbital forcing of Arctic climate: mechanisms of climate response and implications for continental glaciation

Progress in understanding how terrestrial ice volume is linked to Earths orbital configuration has been impeded by the cost of simulating climate system processes relevant to glaciation over orbital time scales (10³–10⁵ years). A compromise is usually made to represent the climate system by models that are averaged over one or more spatial dimensions or by three-dimensional models that are limited to simulating particular snapshots in time. We take advantage of the short equilibration time (10 years) of a climate model consisting of a three-dimensional atmosphere coupled to a simple slab ocean to derive the equilibrium climate response to accelerated variations in Earths orbital configuration over the past 165,000 years. Prominent decreases in ice melt and increases in snowfall are simulated during three time intervals near 26, 73, and 117 thousand years ago (ka) when aphelion was in late spring and obliquity was low. There were also significant decreases in ice melt and increases in snowfall near 97 and 142 ka when eccentricity was relatively large, aphelion was in late spring, and obliquity was high or near its long term mean. These glaciation-friendly time intervals correspond to prominent and secondary phases of terrestrial ice growth seen within the marine ¹⁸O record. Both dynamical and thermal effects contribute to the increases in snowfall during these periods, through increases in storm activity and the fraction of precipitation falling as snow. The majority of the mid- to high latitude response to orbital forcing is organized by the properties of sea ice, through its influence on radiative feedbacks that nearly double the size of the orbital forcing as well as its influence on the seasonal evolution of the latitudinal temperature gradient.     

Three million years of monsoon variability over the northern Sahara

We present a 3 million year record of aeolian dust supply into the eastern Mediterranean Sea, based on hematite contents derived from magnetic properties of sediments from Ocean Drilling Program Site 967. Our record has an average temporal resolution of 400 years. Geochemical data validate this record of hematite content as a proxy for the supply of aeolian dust from the Sahara. We deduce that the aeolian hematite in eastern Mediterranean sediments derives from the eastern Algerian, Libyan, and western Egyptian lowlands located north of the central Saharan watershed (21°N). In corroboration of earlier work, we relate dust flux minima to penetration of the African summer monsoon front to the north of the central Saharan watershed. This would have enhanced soil humidity and vegetation cover in the source regions, in agreement with results from green Sahara climate models. Our results indicate that this northward monsoon penetration recurred during insolation maxima throughout the last 3 million years. As would be expected, this orbital precession-scale mechanism is modulated on both short (100-kyr) and long (400-kyr) eccentricity time scales. We also observe a strong expression of the 41-kyr (obliquity) cycle, which we discuss in terms of high- and low-latitude mechanisms that involve Southern Hemisphere meridional temperature contrasts and shifts in the latitudes of the tropics, respectively. We also observe a marked increase in sub-Milankovitch variability around the mid-Pleistocene transition (0.95 Ma), which suggests a link between millennial-scale climate variability, including monsoon dynamics, and the size of northern hemisphere ice sheets.     

Re-analysis of the nuclear winter phenomenon

Summary An analysis of the report of the (U.S.) National Academy of Sciences (NAS) on atmospheric effects of a nuclear exchange leads to conclusions that differ from those of the NAS and of the earlier TTAPS and AMBIO studies. Any cooling of the earth's surface is likely to beshort-lived because of rapid removal of the smoke clouds originating from nuclear burst-initiated fires, andminor because of appreciable green-house effects due to several distinct physical causes. (One of these, neglected in prior analyses, is the infrared absorption from cirrus clouds produced directly by the nuclear bursts.) Taken together, these effects may even induce slight surface warming (nuclear summer) instead of cooling (nuclear winter). The consequences to atmospheric ozone are similarly ambiguous; depending on the detailed nuclear scenario, the net ozone content may increase-rather than decrease as argued by TTAPS. Experiments could settle some uncertainties.With 2 Figures     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

Bioaccumulation of Cobalt in Parmelia sulcata

Six locations across mainland Portugal were selected for exposing Parmelia sulcata, for a one-year period (8 months for one site), with simultaneous measurement of total (dry + wet) deposition (one-month periods). The exposed lichens and the total (dry + wet) deposition were analysed for cobalt contents by INAA (instrumental neutron activation analysis) and ICP-MS (inductively coupled plasma mass spectroscopy), respectively. The designated wet deposition was evaluated through the collected water volume; the designated dry deposition was assessed after the (dried) residual mass of the wet deposition. An excellent agreement between Co contents in exposed lichens and the cumulative (1) Co contents in the dry deposition, (2) dry deposition, and (3) wet deposition has been found for the locations with alternate drought and precipitation months, high dry deposition, and high Co contents in the latter. Continuous rainfall was found to hinder the Co accumulation in the lichen due to its release from the lichen and/or lower Co contents in the dry deposition. At three locations, P. sulcata Co contents, after subtraction of the background (before exposure), equalled or exceeded the Co contents in the cumulative dry deposition at the end of the exposure time. The optimal exposure period for this species likely depends on the exposure conditions.     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

Thermal comfort of man in different urban environments

Summary On July 29, 1985, a hot summer day, biometeorological measurements were performed simultaneously in three different urban structures within the city of Munich and in the trunk space of a nearby tall spruce forest. Based on the results of these experiments the following thermophysiologically relevant biometeorological indices were calculated: Predicted mean vote, skin wettedness and physiologically equivalent temperature. These three indices are derived from different models for the human energy balance. They allow the assessment of the thermal components of the microclimates at the selected sites with regard to application in urban planning. The results quantitatively show the great heat stress in the urban structure street canyon, exposed to south, whereas in the trunk space of the tall spruce forest there is nearly an optimal climate even on hot summer days. Between these extremes the results for street canyon, exposed to north show a little higher heat load than for backyard with trees.

---

Zusammenfassung An einem heißen Sommertag, dem 29. Juli 1985, wurden in drei Stadtstrukturen in München und im Stammraum eines nahegelegenen Fichtenhochwaldes zeitgleich biometeorologische Messungen durchgeführt. Mit den Meßergebnissen wurden folgende thermophysiologisch relevante biometeorologische Indizes berechnet: Predicted mean vote, Hautbenetzungsgrad und physiologisch äquivalente Temperatur. Diese drei Indizes beruhen auf verschiedenen Modellen zur menschlichen Energiebilanz. Mit den drei Indizes wurden die thermischen Komponenten der Mikroklimate an den ausgewählten Meßplätzen im Hinblick auf Stadtplanungsaufgaben bewertet. Die Ergebnisse zeigen quantitativ die relativ große Hitzebelastung bei der Stadtstruktur Straßenschlucht, nach Süd exponiert, während im Stammraum des Fichtenhochwaldes selbst an heißen Sommertagen nahezu optimale Bedingungen herrschen. Zwischen diesen Extremen liegen die Ergebnisse für die anderen Meßplätze, wobei für Straßenschlucht, nach Nord exponiert die Wärmebelastung etwas höher als für Innenhof mit Bäumen ist.

---

---

With 6 Figures     

Sensitivity of freezing drizzle formation in stably stratified clouds to ice processes

Summary This study examines the impact of ice formation and growth processes on freezing drizzle formation in stably stratified clouds. In particular we investigate the reason why freezing drizzle is rarely observed in clouds with top temperatures less than –15°C. We also investigate the sensitivity of freezing drizzle formation to the Hallett Mossop secondary ice process (Hallet and Mossop, 1974). The evaluation is performed by simulating cloud formation over a two-dimensional idealized mountain using a detailed microphysical scheme. The height and width of the two-dimensional mountain were designed to produce an updraft pattern with extent and magnitude similar to documented freezing drizzle cases. The simulations show that: (i) drizzle formation is very sensitive to the ice crystal concentration, with a significant reduction in the area over which drizzle forms and the maximum drizzle water content as the cloud top temperature decreases below –10°C, and (ii) secondary ice crystal formation has a significant effect on drizzle formation at cloud top temperatures below –10°C.The above two factors are likely the main cause for the lack of freezing drizzle at cloud top temperatures less than –15°C. We also found that neglecting the depletion of ice forming nuclei resulted in considerable overestimation of the ice crystal concentration and suppression of drizzle, even for the –10°C case.     

Numerical case study of the Altai-Sayan lee cyclogenesis over East Asia

Summary Climatological studies show that the Altai-Sayan lee side is one of the major cyclogenesis areas in the Northern Hemisphere. In case of the Altai-Sayan lee cyclogenesis, the surface cyclone is generated when a primary cyclone is swept north of the mountains. In the mid-troposphere, a trough develops and finally turns into a cutoff low within 48 h. The main synoptic features are similar to those of Alpine cyclogenesis. Numerical simulations are performed to assess the effect of different representation of orography on the Altai-Sayan cyclogenesis. Two experiments are performed, a step-mountain (ETA) and an envelope orography (SGM) experiment. The ETA experiment produced the cyclogenesis in a way similar to that in the analysis both at the surface and at mid-troposphere. The SGM experiment failed in the simulation of the upper cutoff low. The difference in predicted pressure between the ETA and the SGM experiment shows a dipolar structure suggesting that the blocking effect of the mountains is essential in the development of the Altai-Sayan lee cyclogenesis.With 7 Figures     

Recurrence statistics of concentration fluctuations in plumes within a near-neutral atmospheric surface layer

This study examines the statistical properties of the concentration derivative, , for a dispersing plume in a near-neutrally stratified atmospheric surface layer. Towards this goal, the probability density function (pdf) of , and the conditional pdf of given a fixed concentration level, , have been measured. These pdfs are found to be modeled well by a generalizedq-Gaussian (gqG) distribution with intermittency exponent,q, equal to 0.3 and 3/4, respectively. These results highlight the strong intermittency effect (patchiness) of the small-scale concentration eddy structures in the plume. The distribution of time intervals between successive high peaks in the squared derivative process, x², is found to be well approximated by a power-law distribution, implying that occurrences of these high peaks are much more clustered than would be predicted by a Poisson or shot-noise process. The results are used to improve models for the joint pdf of and , and for the expected number of upcrossings per unit time interval of a fixed concentration level that have been proposed by Kristensenet al. (1989). The predictions of the improved models are in accord with observations, and suggest that the intercorrelation between and must be explicitly incorporated if good estimates of the upcrossing intensity are to be obtained.     

Time-dependent greenhouse warming computations with a coupled ocean-atmosphere model

Climate changes during the next 100 years caused by anthropogenic emissions of greenhouse gases have been simulated for the Intergovernmental Panel on Climate Change Scenarios A (business as usual) and D (accelerated policies) using a coupled ocean-atmosphere general circulation model. In the global average, the near-surface temperature rises by 2.6 K in Scenario A and by 0.6 K in Scenario D. The global patterns of climate change for both IPCC scenarios and for a third step-function 2 x CO₂ experiment were found to be very similar. The warming delay over the oceans is larger than found in simulations with atmospheric general circulation models coupled to mixed-layer models, leading to a more pronounced land-sea contrast and a weaker warming (and in some regions even an initial cooling) in the Southern Ocean. During the first forty years, the global warming and sea level rise due to the thermal expansion of the ocean are significantly slower than estimated previously from box-diffusion-upwelling models, but the major part of this delay can be attributed to the previous warming history prior to the start of present coupled ocean-atmosphere model integration (cold start).     

General circulation model experiments with surface albedo changes

In 1975, Charney proposed a biogeophysical feedback mechanism to partly explain the droughts that occur in desert border areas. He showed that a perturbation of albedo (due to a natural or anthropogenic decrease of vegetation) can be unstable and lead to a variation of precipitation in the region where albedo is changed.Several numerical experiments have been achieved with general circulation models to study the sensitivity of climate to surface albedo. We compare the GLAS and LMD model results for the Sahel. For all models, rainfall decreases when albedo increases and net radiative heating of soil decreases. We show the variations of circulation simulated by the LMD model that we obtain when albedo is increased. These changes are compared to the weakening of Easterly Jet at 200 mb observed during dry years.     

Surface-Layer Fluxes in Stable Conditions

Micrometeorological tower data from the Microfronts experiment are analyzed. Scale-dependencies of the flux and flux sampling error are combined to automatically determine Reynolds turbulence cut-off time scales for computing fluxes from time series. The computed downward heat flux at the 3 m height averaged over nine nights with 7.3 hours each night is 20% greater than the downward heat flux computed at the 10 m height. In contrast, there is only a 1.2% difference between 3 m and 10 m heat fluxes averaged over daytime periods, and there is less than a 2% difference between 3 m and 10 m momentum fluxes whether averaged over nighttime or daytime periods.Stability functions, M(z/L) and H(z/L) are extended to z/L up to 10, where z is the observational height and L is the Obukhov length. For 0.01 < z/L < 1 the estimated functions generally agree with Businger-Dyer formulations, though the H estimates include more scatter compared to the M estimates. For 1 < z/L < 10, the flux intermittency increases, the flux Richardson number exceeds 0.2, and the number of flux samples decreases. Nonetheless the estimates of the stability function M based on 3-m fluxes are closer to the formula proposed by Beljaars and Holtslag in 1991 while the M functions based on 10-m fluxes appears to be closer to the formula proposed by Businger et al. in 1971. The stability function H levels off at z/L = 0.5.     

The Effect of Atmospheric Stability on the Surface-Layer Characteristics in a Low-Wind Area of Tropical West Africa

This study details the observed effects ofatmospheric stability on characteristics of thesurface layer in a low wind speed (U = 1.5 m s-1)regime of tropical West Africa. Theaerodynamic roughness length, z0, anddisplacement height, d, obtained from profilewind-speed data at our bush land site (height 2 m)have values of 0.24 ± 0.10 m and 1.54 ± 0.04 mrespectively. In the unstable range (-2.5 < Ri < -0.1; Riis gradient Richardson number), thestandard deviation in wind speed fluctuations, u, increased from 0.57 ± 0.19 m s-1 toa maximum of 0.7 ± 0.2 m s-1 in near-neutralconditions, and in the stable range, the parameterdecreased rapidly to 0.41 ± 0.15 m s-1 at Ri 0.2.In the same stability range, the horizontal winddispersion, , decreased withincreasing stability from 19 ± 8 deg. to 13 ± 5 deg.The surface-layer integral quantity, u/u*, when plottedas a function of stability, is in agreement with theempirical results. The ratio ofsensible heat flux (estimated) to the net radiationranged between 0.1 and 0.2 at nighttime,increasing to about 0.5 during the daytime, and showeda strong dependency on season.     

Trace gas exchange at the air/water interface: Measurements of mass accommodation coefficients

A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO₃ (0.01), N₂O₅ (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.     

A new method of cluster analysis for numerical classification of climate

Summary In this paper, a new method of cluster analysis incorporating a statistical test, was proposed in order to decide whether to accept or reject the assumption of a numerical classification of climate according to a ² test. On the basis of normality of errors, the ² test has been derived from the variance of the Euclidean distance or from that of the reduced distance.This new method of cluster analysis consists of uniform clustering and central clustering. The former is adapted to a more detailed numerical classification of climate, while the latter is appropriate for the large-scale numerical classification of climate. Examples of calculation and application of this new method of cluster analysis are presented in order to illustrate calculations undertaken in the processes of climate classification.     

Charakteristische Zirkulationstypen in mittleren Breiten der nördlichen Hemisphäre

Zusammenfassung In der vorliegenden Arbeit werden neue Klassifikationsprinzipien für Großwetterlagen entwickelt. Bisher wurde bei Wetterlagenklassifikationen das Druckfeld zugrunde gelegt, wobei quasistationären Druckzentren eine nicht berechtigte Vorrangstellung eingeräumt wurde. In der hier versuchten Klassifikation wird vom Strömungsfeld ausgegangen, das in elementare Formen zerlegt wird. Eine zu diesem Zweck durchgeführte statistische Untersuchung ergab, daß alle im Strömungsfeld auftretenden Zirkulationstypen auf drei Grundformen zurückgeführt werden können. Diese Grundformen sind:Driften, Wellen undWirbel.Die Untersuchung ergab im einzelnen, daß in mittleren Breiten der nördlichen Hemisphäre bei 49% aller untersuchten Fälle Driften, bei 23% Wellen und bei 28% Wirbel auftraten.In der hier durchgeführten Klassifikation wird das Druckfeld durch das Strömungsfeld und der Begriff Großwetterlage durch den umfassenderen Begriff des Zirkulationstyps ersetzt. Damit wird der unberechtigte Vorrang der Druckformen bei der Wetterlagenklassifikation aufgegeben. Die Klassifizierung der Zirkulationstypen ergibt sich schließlich durch Kombination der drei Zirkulationselemente: Drift, Welle und Wirbel.

---

Summary New principles of classification for large-scale weather situations are outlined in this paper. Hitherto the pressure-field has been taken as a basis for such classifications by conceding a precedence of an unjustified position to semi-permanent centres of pressure. The new classification starts from the field of large-scale motions, which is dissected in elementary models. A statistical test yielded the possibility to reduce all types of atmospheric circulations in the following three elementary models:drifts, waves andeddies.In detail it was found out, that drifts occur in 49%, waves in 23% and eddies in 28% of all cases investigated.In the new classification the term pressure-field is substituted by field of motion and the expression large-scale weather situation by the more comprehensive conception type of circulation. By that the unjustified priority of pressure-centers in classifying weather situations is abolished. At last the classification of the types of circulation follows from a combination of the three elementary models: drift, wave and eddy.

---

résumé La présente étude développe de nouveaux principes de classification des situations météorologiques. Alors que jusqu'ici on s'est fondé sur le champ de pression ce qui conduisait à attribuer aux centres d'action quasi stationnaires un rôle trop important, l'auteur part ici du champ de mouvement décomposé en formes élémentaires. Un examen statistique lui a montré que tous les types de circulation peuvent se ramener à trois formes fondamentales:courants, ondulations ettourbillons.Aux latitudes moyennes de l'hémisphère Nord les courants représentent le 49%, les ondulations le 23% et les tourbillons le 28%.Au champ de pression se substitue donc le champ de courant, et les situations météorologiques se groupent en types de circulation ce qui supprime le rôle prépondérant des formes isobariques. Le classement final des types de circulation résulte de la combinaison des trois types mentionnés:courants, ondulations ettourbillons.

---

---

Mit 15 Textabbildungen     

Application of a “Big-Tree” model to regional climate modeling: a sensitivity study

Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m²) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m², respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m² respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.