Arctic hazes in summer over Greenland and the North American Arctic: II. Nature and concentrations of accumulation-mode and giant particles

Airborne measurements made during August 1985 over Greenland and its environs show that both accumulation-mode (0.1 m D2.0 m) and giant (D2 m) particles were present in relatively high concentrations in arctic haze layers and that the accumulation-mode particles dominated light scattering. Particles with diameters (D) between 1 and 4 m consisted predominately of mixed materials, small and dense inclusions, and probably organic compounds containing sulfur. Many of the particles from 0.1 to 1 m in diameter were also of mixed composition, with sulfuric acid, ammonium sulfate and organics probably the dominant constituents.     

Nitrate in the atmospheric boundary layer of the tropical South Pacific: Implications regarding sources and transport

Weekly bulk aerosol samples collected at Funafuti, Tuvalu (8°30S, 179°12E), American Samoa (14°15S, 170°35W), and Rarotonga (21°15S, 159°45W), from 1983 through most of 1987 have been analyzed for nitrate and other constituents. The mean nitrate concentration is about 0.11 g m^–3 at each of these stations: 0.107±0.011 g m^–3 at Funafuti; 0.116±0.008 at American Samoa; and 0.117±0.010 at Rarotonga. Previous measurements of mineral aerosol and trace metal concentrations at American Samoa are among the lowest ever recorded for the near-surface troposphere and indicate that this region is minimally affected by transport of soil material and pollutants from the continents. Consequently, the nitrate concentration of 0.11 g m^–3 can be regarded as the natural level for the remote marine boundary layer of the tropical South Pacific Ocean. In contrast, over the tropical North Pacific which is significantly impacted by the transport of material from Asia and North America, the mean nitrate concentrations are about three times higher, 0.29 and 0.36 g m^–3 at Midway and Oahu, respectively. The major sources of the nitrate over the tropical South Pacific are still very uncertain. A very significant correlation between the nitrate concentrations at American Samoa and the concentrations of ²¹⁰Pb suggests that transport from continental sources might be important. This continental source could be lightning, which occurs most frequently over the tropical continents. A near-zero correlation with ⁷Be indicates that the stratosphere and upper troposphere are probably not the major sources. A significant biogenic source would be consistent with the higher mean nitrate concentrations, 0.16 to 0.17 g m^–3, found over the equatorial Pacific at Fanning Island (3°55N, 159°20W) and Nauru (0°32S, 166°57E). The lack of correlation between nitrate and nss sulfate at American Samoa does not necessarily preclude an important role for marine biogenic sources.     

On the determination of the height of the Ekman boundary layer

The heighth of the Ekman turbulent boundary layer determined by the momentum flux profile is estimated with the aid of considerations of similarity and an analysis of the dynamic equations. Asymptotic formulae have been obtained showing that, with increasing instability,h increases as ¦¦^1/2 (where is the non-dimensional stratification parameter); with increasing stability, on the other hand,h decreases as ^–1/2. For comparison, a simple estimate of the boundary-layer heighth _u determined by the velocity profile is given. As is shown, in unstable stratification,h _u behaves asymptotically as ¦¦^–1, i.e., in a manner entirely different from that ofh .     

Elemental Composition of Mineral Aerosol Generated from Sudan Sahara Sand

Eighteen soil samples from central Sudan were fractionated by dry sieving ina size fraction from <45 m to >300 m while aerosols generatedfrom these soils were fractionated in the particle size range from 0.25 mto >16 m. The elemental concentrations of soil samples were determinedby energy-dispersive X-ray fluorescence, while the elemental concentrationsof generated aerosols were analysed by particle-induced X-ray emission. Theelements Al, K and Rb show a slight positive fractionation with decreasingparticle size throughout the particle size range studied. The concentrationsof Ca, Mn, Fe, Sr and Y are maximum in the small soil size fraction (<45m) and decrease for the coarse soil size fractions, while in the mineralaerosol particle sizes (0.25– > 16 m) the concentrations remainmore or less constant. The size distributions for Cr, Ti and Zr show a maximumin the particle size range 45–100 m and the concentrations of theseelements decrease sharply in the aerosol fraction down to 16 m to remainconstant in the smaller aerosol fractions.Enrichment factors for the elements were calculated relative to five referencematerials: average crustal rock, average soil, the investigated Sahara bulksoil, the finest fraction of this soil and the aerosol generated from thissoil, and using four reference elements: Al, Si, Ti and Fe. The enrichmentfactors were found to vary significantly depending on the choice of thereference material or the reference element. The enrichment factors for theSudan mineral aerosol were almost identical to those for Khartoum atmosphericaerosol but different from those for Namib mineral aerosol and Israelatmospheric aerosol following dust storms. Multivariate display methods(cluster analysis, principal component analysis and linear discriminantanalysis) were applied to the element ratios in the mineral aerosol from theSahara and Namib and this showed that these mineral aerosol can bedifferentiated into different groups. An attempt was also made to relate themineral aerosol to its parent soil through the use of these multivariatetechniques and the elemental ratios in both the mineral aerosols and the bulksoils (Namib and Sahara). It was also possible using the elemental ratios andthe multivariate display methods to associate the crustal component to themineral aerosol generated from the Sahara.     

On the similarity in atmospheric fluctuations of carbon dioxide,water vapor and temperature over vegetated fields

This paper describes the similarity between atmospheric fluctuations of carbon dioxide, water vapor and temperature using data which cover a wide range of instability (0.02 < < 10). The is the Monin-Obukhov stability parameter including the humidity effect.The spectral analysis shows that the coherency between fluctuations of carbon dioxide and water vapor or temperature is very close to unity, and the phase difference is basically out of phase for whole frequency ranges analyzed. The stability dependence of the normalized standard deviation of carbon dioxide is very similar to those of water vapor and temperature. The normalized standard deviation is about 2.5 under near neutral conditions, and it decreases with increasing instability following the -1/3; power law as (-)^-1/3. The skewness factors of carbon dioxide, water vapor and temperature show a systematic departure with increasing instabilities for 0.02 < s- < 1, and level off at high instabilities for 1 < -\s < 10. The stability dependence of the flatness factors is not so clear as that noted in the standrard deviation and skewness factors. Dissipation rates of carbon dioxide, water vapor and temperature variance are well related to the spectral peak wavelength. This seems to be real since the local production and local dissipation rates are the main terms, almost balancing one another in the variance budget equations for scalar entities.     

Evaporation from a blanket bog in a foggy coastal environment

Frequent fog severely restricts evaporation from blanket bogs in Newfoundland because it more than halves the radiant energy input, and it eliminates the vapor pressure deficit, resulting in evaporation at the equilibrium rate (average = 0.99 during fog). During these periods, there is no surface resistance to evaporation because the bog has been wetted by fog drip, and although the latent heat flux dominates over sensible heat (average = 0.8), both are small. In contrast, the surface dries during clear periods, increasing the surface resistance to evaporation so that sensible heat becomes more important ( = 1.05). When the mosses are dry, evaporation is below the equilibrium rate ( = 0.87), although the higher available energy ensures that actual evaporation is higher. During clear periods, daily evaporation averaged 2.5 mm, compared to 1.1 and 0.7 mm for fog and rain, respectively. The suppressed evaporation at this site is important in maintaining appropriate hydrological conditions for blanket bog development.     

Convective Profile Constants Revisited

This paper examines the interpolation betweenBusinger–Dyer (Kansas-type) formulae,_u = (1 -1 6 )^-1/4 and_t = (1 - 16 )^-1/2, and free convection forms. Based on matching constraints, the constants, a_u and a_t, in the convective flux-gradient relations, _u = (1 - a_u )^-1/3 and _t = (1 - a_t )^-1/3, are determined. It isshown that a_u and a_t cannot be completely independent if convective forms are blended with theKansas formulae. In other words, these relationships already carryinformation about a_u and a_t. This follows because the Kansas relations cover a wide stability range (up to = - 2), which includes a lower part of the convective sublayer (about 0.1 < - < 2). Thus, there is a subrange where both Kansas and convective formulae are valid. Matching Kansas formulae and free convection relations within thesubrange 0.1 < - < 2 and independently smoothing ofthe blending function are used to determine a_u and a_t. The values a_u = 10 for velocity and a_t = 34for scalars (temperature and humidity) give a good fit. This new approacheliminates the need for additional independent model constants and yields a`smooth' blending between Kansas and free-convection profileforms in the COARE bulk algorithm.     

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.     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

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     

Size-Resolved Characterisation of Soluble Ions in the Particles in the Tropospheric Plume of Masaya Volcano, Nicaragua: Origins and Plume Processing

We present the first application of a multi-stage impactor to study volcanic particle emissions to the troposphere from Masaya volcano, Nicaragua. Concentrations of soluble SO₄ ^2–,Cl^–, F^–, NO₃ ^–, K⁺, Na⁺,NH₄ ⁺, Ca²⁺ and Mg²⁺ were determined in 11 size bins from 0.07 m to >25.5 m. The near-source size distributions showed major modes at 0.5m (SO₄ ^2–, H⁺,NH₄ ⁺); 0.2 m and 5.0 m (Cl^–) and 2.0–5.0 m(F^–). K⁺ and Na⁺ mirrored the SO₄ ^2– size-resolvedconcentrations closely, suggesting that these were transported primarily asK₂SO₄ and Na₂SO₄ in acidic solution, while Mg²⁺ andCa²⁺ presented modes in both <1 m and >1 m particles. Changes in relative humidity were studied by comparing daytime (transparent plume) and night-time (condensed plume) results. Enhanced particle growth rates were observed in the night-time plume as well as preferential scavenging of soluble gases, such as HCl, by condensed water. Neutralisation of the acidic aerosol by background ammonia was observed at the crater rim and to a greater extent approximately 15 km downwind of the active crater. We report measurements of re-suspended near-source volcanic dust, which may form a component of the plume downwind. Elevated levels ofSO₄ ^2–, Cl^–, F^–,H⁺, Na⁺, K⁺ and Mg²⁺ were observed around the 10 m particle diameter in this dust. The volcanic SO₄ ^2– flux leaving the craterwas 0.07 kg s^–1.     

Final results of the CONDORS convective diffusion experiment

TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF ₆ values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, _y and _z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations _vand _w 0.6w^* (the convective scale velocity), as measured atz>100m. Surface release _y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger _v and reduced wind speed measured near the surface. Decreased _y at small dimensionless average times is also noted. Finally, convectively scaled dy, C _y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC _y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.     

Free convection frequency spectra of atmospheric turbulence over the sea

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

Concentration of Peroxides and Formaldehyde in Air and Rain and Gas-Rain Partitioning

Rain and air of Florence have been collected in a continuous way andanalysed by flow analysis spectrofluorimetric methods for formaldehydeand hydrogen peroxide. Diurnal and seasonal variations were observed;the mean/maximum concentrations of all data (as gm^–3) are 3.3/23.4 for HCHO and 0.4/4.93 forH₂O₂. The effect of external sources and ofphotochemical reactions produces periods of positive and negativecorrelations for these compounds. The mean/maximum rain concentration ofall data are 98/443 g l^–1 for HCHO and 84/685 g l^–1 for H₂O₂. Concentrationratios rain/air and discrepancies to Henry's Law equilibrium arediscussed.     

An Alternative Function For The Wind And Temperature Gradients In Unstable Surface Layers

The function ()=(1+|z/L|^2/3)^1/2,where z is the height, L the Obukhov length, and a constant,is proposed for the nondimensional wind speed and temperaturegradients (flux-profile relationships) in anunstable surface layer. This function agrees quite well withboth wind speed and temperature data,has the theoretically correct behaviour in convective conditions,and leads to simple results when integrated to produce the mean profiles.     

Zur Struktur der Wolken

Zusammenfassung Auf dem Feldberg (1493 m NN) im südlichen Schwarzwald wurden von Januar 1953 bis Mai 1953 und von September 1954 bis Februar 1955 Messungen von Tropfengröße, Wassergehalt und Sichtweite in Wolken durchgeführt. Die Tropfen der einzelnen Spektren haben einen mittleren Durchmesser zwischen 6 und 30 ; die Durchmesser der größten gemessenen Tropfen liegen zwischen 50 und 60 . Der Wassergehalt hat Werte zwischen 0,003 und 1,54 g/m³, wobei über 90% aller Messungen einen Wassergehalt unter 0,50 g/m³ zeigen. Ein einfacher Zusammenhang zwischen Sichtweite, Tropfengröße und Wassergehalt nach derTrabertschen Formel besteht nicht.Eine enge beziehung wurde zwischen der Form der Spektren und dem Niederschlag gefunden. In Zeiten ohne Niederschlag und bei geringfügigem Niederschlag überwiegen schmale ( < 12,0 ), in Zeiten mit länger andauerndem Niederschlag breite Spektren ( > 12,0 ). Eine Ausnahme machen alle Messungen in Hochnebel, die breite Spektren mit zwischen 13 und 17 zeigen, ohne daß Niederschlag fällt.Das gleichzeitige Vorhandensein der drei Phasen Wasserdampf, flüssiges Wasser und Eisteilchen über längere Zeiträume hinweg wurde festgestellt; somit kann die Eiskristall-Niederschlagstheorie vonBergeron und vonFindeisen für den Mechanismus der Niederschlagsbildung nicht im allgemein angenommenen Ausmaß zutreffen. Es erscheint notwendig, die Frage der Wolkenphysik nicht nur thermodynamisch, sondern auch in vermehrtem maß kolloidphysikalisch im Sinne vonSchmauss zu untersuchen.

---

Summary Measurements of drop size, water content and visibility in clouds have been carried out from January till May 1953 and from September 1954 to February 1955 on the Feldberg (1493 m. a. s. 1.) in the southern Black Forest. The drops of the different spectra have a mean diameter between 6 and 30 ; the diameters of the largest measured drops lie between 50 and 60 . The water content shows values between 0.003 and 1.54g/m³ whereby more than 90 per cent of all measurements lead to a water content below 0.50 g/m³. A simple relation between visibility, drop size and water content according to the formula ofTrabert does not exist.A closed connection between the form of the spectra and the precipitations has been found. In periods without or with only slight precipitation narrow spectra ( < 12.0 ) prevail, in periods with persistent precipitation broad ones ( > 12.0 ). All measurements in elevanted fog make an exception, showing broad spectra with between 13 and 17 without any precipitation.The simultaneous existance of the 3 phases: water vapour, liquid water and ice particles over alonger periods has been stated; the ice crystall precipitation theory ofBergeron andFindeisen for the mechanism of precipitation can therefore not hold in the generally accepted extent. The examination of the problems of cloud physics appears to be necessary, not only in thermodynamical respect, but also in the sense of coloid physics according toSchmauss.

---

Résumé On a effectué de janvier à mai 1953 et de septembre 1954 à février 1955 au Feldberg (1493 m.) en Forêt-Noire des mesures de grosseur de gouttelettes, de teneur en eau et de visibilité dans les nuages. Les gouttelettes ont und diamètre moyen compris entre 6 et 30 ; les plus grosses mesurent 50 à 60 . La teneur en eau varie de 0,003 à 1,54 g/m³; dans plus de 90% des mesures, elle reste inférieure a 0,5g/m³. Il n'y a pas de relation simple, qui suivrait la formula deTrabert, entre visibilité, grosseur des gouttelettes et teneur en eau.Il existe par contre une relation étroite entre la forme du spectre de fréquence et les précipitations. Lorsque celles-ci sont nulles ou très faibles ce sont les spectres étroit ( < 12,0 ) qui dominent; par précipitations durables, ce sont les spectres larges ( > ). Dans le stratus, sans pluie, les spectres larges ( entre 13 et 17 ) prédominent.On a constaté l'existence simultanée des trois phases solide, liquide et gazeuse pendant d'assez longs intervalles de temps, de sorte que le théorie deBergeron et deFindeisen concernant la genène des précipitations n'a pas la généralité qu'on lui prête. Il apparaît nécessaire d'aborder les problèmes de la physique des nuages non seulement du point de vue thermodynamique, mais aussi de celui de l'état colloïdal au sens deSchmauss.

---

---

Mit 9 Textabbildungen.

---

Hinsichtlich der Wiedergabe von Tabellen und Tropfenspektren stark verkürzte Fassung einer bei der Technischen Hochschule in Karlsruhe eingereichten und genehmigten Dissertation; die hier weggelssene Dokumentation kann zur Einsichtnahme bezogen werden bei den Bibliotheken der Technischen Hochschule Karlsruhe und des Deutschen Wetterdienstes, Offenbach (Main).     

Infrared observations and numerical modelling of grassland fires in the Northern Territory, Australia

Summary This paper reports on a small-scale pilot experiment held early in the dry season near Darwin, Australia, in which fine-scale observations of several prescribed fires were made using infrared digital video. Infrared imaging is used routinely to locate fires as infrared radiation suffers little attenuation as it propagates through the smoke that normally obscures visible imagery. However, until now, little use has been made of digital video imagery in analyzing the convective-scale structure of prescribed (or wild) fires. The advantage of digital video imagery is that the individual frames can be objectively analyzed to determine the convective motion in the plane viewed by the camera. The infrared imagery shows mostly rising plumes, much like convective clouds. The flow is highly convective, and the vertical transport of heat is confined to relatively narrow thermals. The updrafts range from a few ms^–1 to around 15ms^–1. A numerical model is used to simulate one of the prescribed fires at very high-resolution. For the most part, the model predictions compare well to the observations. The model produces plumes that are around 7m high, and spaced around 5m apart, which is similar to that observed. The model correctly predicts the mean rate of spread of the fire to be 1.3ms^–1. Perhaps the most serious limitations to using infrared observations of the type presented here are the difficulties in interpreting precisely the relationship between the observed infrared temperature field and the air temperature calculated by the model, and the exact connection between the infrared camera derived flow field and that calculated by the model.     

Parameterizing turbulent diffusion through the joint probability density

The convective mass flux parameterization often used in meteorological modeling expresses the vertical flux of a transported scalar as proportional to the product of the difference in mean values of the scalar in updrafts and downdrafts and their characteristic velocity. The proportionality factor is a constant to be specified. We show that this proportionality factor also appears in the relaxed eddy accumulation technique of Businger and Oncley. That associates the surface-layer flux of a scalar with the product of the standard deviation of vertical velocity and the mean concentration difference between updrafts and downdrafts.We show that this constant (b) is determined uniquely by the joint probability density (jpd) of vertical velocity and the scalar. Using large-eddy simulation, we generate this jpd for a conservative scalar diffusing through a convective boundary layer. It has quite different forms in top-down and bottom-up diffusion geometries. The bottom-up jpd is fairly well represented by a jointly Gaussian form and implies b ~ 0.6, in good agreement with the surface-layer value reported by Businger and Oncley. The top-down jpd is strikingly non-Gaussian and gives b ~ 0.47. Updrafts carry the bulk of the scalar flux - 70% in the bottom-up case, 60% in the top-down case.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

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.