AN EVALUATION OF DISPERSION COEFFICIENTS FOR USE IN AIR QUALITY MODELS

Standard deviations of concentration in horizontal andvertical directions i.e. y andz have been estimated by using fivedifferent schemes based on empirical(due to Pasquill and Briggs)schemes and sophisticated methods(due to Irwin, Draxler, Taylor, Hanna et al.). The fiveschemes are discussed atlength. The purpose of this study is to make use ofmeteorological observations whichare routinely available, to test all the above methods andintercompare the resultswith one another and observations so that the sensitivityof each schemeunder various atmospheric stability conditions could beassessed. It has beenfound that the existing schemes are good enough to providereasonable estimates ofdispersion coefficient (y) during highly unstableconditions (Pasquill stability classes A and B). However, thesame is not true for the case when the stability increasesfrom C to F and turbulencedecreases, specifically during stable atmospheric conditions,when the observedvalues were found to be much higher than all the existingschemes. This suggests thatwhile we continue to use the current methods of estimatingthe dispersion parameters,a rigorous search is required for methods which give predictionswhich are close-to-realityduring such conditions which are represented by lowlevels (in terms of magnitude)of atmospheric turbulence leading to higher levelsof pollution.As one of the sophisticated methods requiresthe use of v and w (standard deviationsof wind velocity fluctuation in y and z directions),these have been estimated andvalidated with observed data (field experiments conductedby EPRI at Kincaid).Statistical evaluation of v and wbased on performance measures indicate a goodperformance of the parameterisations adopted in thisstudy. In the case of w duringunstable conditions a comparison of three differentschemes with observations is made.     

Sodar Detected Top-Down Convection in a Nocturnal Cloud-Topped Boundary Layer: A Case Study

Nocturnal convection, originating in a well-mixed marine cloud-topped boundary layer, advected onshore, was observed using a Doppler sodar on the Tyrrhenian coast in Italy. The horizontal and vertical dimensions of the downdrafts were evaluated. The oscillation frequency triggered by the downdrafts at the inversion layer, derived from the harmonic analysis of the sodar measured vertical velocity (w), is compared with the Brunt-Vaisala frequency, obtained from the rawinsonde temperature profile. A similarity function for the ²_w vertical profile was used to fit the sodar experimental data and to retrieve the depth of the mixing layer and the sensible heat flux at the top of the cloud layer. The results are in agreement with the convection layer depth observed in the sodar echoes facsimile record, and with the energy budget evaluated at the top of the cloud layer using the rawinsonde profiles.     

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.     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.     

An atmospheric dispersion model with continuous vertical variation of dispersion parameters

Summary A dispersion model is proposed to predict the continuous vertical variation of the dispersion parameters _y and _z in case of hot pollutant release to the atmosphere. In such a case, the plume rises far above the ground and is subject to varying levels of turbulence. The framework in this paper can be divided into three approaches: (1) determination of the eddy diffusivitiesK _y (z, _y ) andK _z (z, _z ) as functions of height above ground and plume dimensions, (2) determination of both the plume rise and its vertical velocity using a modified version of Brigg's formula, and (3) numerical solution of actual problems with buoyant plumes at each time step. The model results have been applied to a case of pollutant release from fire destruction of a chemical storehouse roof.With 15 Figures     

Spatial Variability of Turbulent Fluxes in the Roughness Sublayer of an Even-Aged Pine Forest

The spatial variability of turbulent flow statistics in the roughness sublayer (RSL) of a uniform even-aged 14 m (= h) tall loblolly pine forest was investigated experimentally. Using seven existing walkup towers at this stand, high frequency velocity, temperature, water vapour and carbon dioxide concentrations were measured at 15.5 m above the ground surface from October 6 to 10 in 1997. These seven towers were separated by at least 100m from each other. The objective of this study was to examine whether single tower turbulence statistics measurements represent the flow properties of RSL turbulence above a uniform even-aged managed loblolly pine forest as a best-case scenario for natural forested ecosystems. From the intensive space-time series measurements, it was demonstrated that standard deviations of longitudinal and vertical velocities (_u, _w) and temperature (_T) are more planar homogeneous than their vertical flux of momentum (u_* ²) and sensible heat (H) counterparts. Also, the measured H is more horizontally homogeneous when compared to fluxes of other scalar entities such as CO₂ and water vapour. While the spatial variability in fluxes was significant (>15 %), this unique data set confirmed that single tower measurements represent the canonical structure of single-point RSL turbulence statistics, especially flux-variance relationships. Implications to extending the moving-equilibrium hypothesis for RSL flows are discussed. The spatial variability in all RSL flow variables was not constant in time and varied strongly with spatially averaged friction velocity u_*, especially when u_* was small. It is shown that flow properties derived from two-point temporal statistics such as correlation functions are more sensitive to local variability in leaf area density when compared to single point flow statistics. Specifically, that the local relationship between the reciprocal of the vertical velocity integral time scale (I_w) and the arrival frequency of organized structures (/h) predicted from a mixing-layer theory exhibited dependence on the local leaf area index. The broader implications of these findings to the measurement and modelling of RSL flows are also discussed.     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

The behaviour of the atmospheric aerosol scattering coefficient under varying meteorological conditions

Summary The atmospheric aerosol scattering coefficient _s , measured for more than a year more or less continuously in Vienna, Austria, exhibits unexpected patterns of variation. Apart from the usual ones following changes in relative humidity or traffic characteristics,a distinctive pattern is found before a change in air mass. _s rises by a factor of 1.5 to 2 some hours (usually two or three) before the passage of the front without a corresponding change in emission characteristics or relative humidity and then falls either below or to its previous level. This behaviour of _s occurred at all frontal passages during the sampling period at all times of day and of year except when the wind speeds were very high.An explanation is attempted by examining the mixing heights before a change in airmass since a reduced vertical dispersion due to pre-frontal changes of stability could account for the increase in _s (and thus the aerosol concentration). It has been found that calculated mixing heights are reduced by nearly the same factor as the value of _s is increased before the front. After the front the factors are similar, but then the aerosol concentration depends also on the origin of the air mass.

---

Zusammenfassung Der Streukoeffizient des atmosphärischen Aerosols ( _s ) wurde in Wien mehr als ein Jahr lang mehr oder weniger kontinuierlich gemessen. Dabei zeigten sich unerwartete Änderungen. Abgesehen vom üblichen Tagesgang im Zusammenhang mit Verkehr und relativer Feuchte fand sich vor einem Luftmassenwechsel ein charakteristischer zeitlicher Verlauf. _s steigt einige Stunden (meist zwei oder drei) vor dem Frontdurchgang an, ohne daß sich die relative Feuchte oder die Quellencharakteristik entsprechend ändert, und fällt dann entweder auf oder unter seinen ursprünglichen Wert. Dieses Verhalten trat bei allen Frontdurchgängen zu jeder Tages- und Jahreszeit auf. Die einzigen Ausnahmen waren Fronten mit sehr hohen Windgeschwindigkeiten.In dieser Arbeit wird versucht, das Verhalten des Streukoeffizienten (und damit der Aerosolkonzentration) durch eine Betrachtung der Mischungshöhen zu erklären, da eine Reduktion der vertikalen Ausbreitung durch Stabilitätsänderungen vor der Front den Anstieg von _s bewirken könnte. Eine Berechnung der Mischungshöhen ergab, daß sie vor der Front um fast denselben Faktor abnahmen um den _s anstieg. Nach der Front waren die Änderungsfaktoren einander ähnlich, obwohl die Aerosolkonzentration auch vom Ursprung der Luftmasse abhing.

---

---

With 5 Figures     

The response of a propeller anemometer to turbulent flow with the mean wind vector perpendicular to the axis of rotation

The system transfer function ¦H(v)¦² at frequencyv (units of Hz) for a vertical velocity propeller anemometer in a statistically stationary and horizontally homogeneous turbulent flow is determined from: (1) experimental estimates of propeller velocity spectra; and (2) estimates of Eulerian vertical velocity spectra based on the hypothesis that degradation of the input vertical velocity Fourier components occurs in the inertial subrange. The experimental estimates of ¦H(v)¦² were adequately summarized with the mathematical expression for the system transfer function of a first-order system with parameterT which has units of time and is analogous to the time constant of a horizontal velocity propeller anemometer. Dimensional analysis techniques and the Monin-Obukhov similarity hypothesis were used to construct a model for the system parameterT which yielded the result that _w /D ₁ ( _w /)^1/3, where _w , andD ₁ denote the standard deviation of the input vertical velocity fluctuations, the horizontal mean wind speed, and the diameter of the propeller, respectively. The system parameterT is interpreted in terms of the time required for the propeller velocity statistics to become asymptotically independent of time upon being released from rest in a statistically stationary turbulent flow.Currently on leave of absence from the Indian Institute of Technology, New Delhi, India.     

Short-Term Variability of Atmospheric DMS and Its Oxidation Products at Amsterdam Island during Summer Time

A one-month experiment was performed at Amsterdam Island in January 1998, to investigate the factors controlling the short-term variations of atmospheric dimethylsulfide (DMS) and its oxidation products in the mid-latitudes remote marine atmosphere. High mixing ratios of DMS, sulfur dioxide (SO₂) and dimethylsulfoxide (DMSO) have been observed during this experiment, with mean concentrations of 395 parts per trillion by volume (pptv) (standard deviation, = 285, n = 500), 114 pptv ( = 125, n = 12) and 3 pptv ( = 1.2, n = 167), respectively. Wind speed and direction were identified as the major factors controlling atmospheric DMS levels. Changes in air temperature/air masses origin were found to strongly influence the dimethylsulfoxide (DMSO)/DMS and SO₂/DMS molar ratios, in line with recent laboratory data. Methanesulfonic acid (MSA) and non-sea-salt sulfate (nss-SO₄ ^2–) mean concentrations in aerosols during this experiment were 12.2± 6.5 pptv (1, n=47) and 59 ± 33 pptv (1, n=47), respectively. Evidence of vertical entrainment was reported following frontal passages, with injection of moisture-poor, ozone-rich air. High MSA/ nss-SO₄ ^2– molar ratios (mean 0.44) were calculated during these events. Finally following frontal passages, few spots in condensation nuclei (CN) concentration were also observed.     

Variability of surface roughness and turbulence intensities at a coastal site in India

Surface-layer features with different prevailing wind directions for two distinct seasons (Southwest Monsoon and Northeast Monsoon) on the west coast of India are studied using data obtained from tower-based sensors at a site located about 500 m from the coast. Only daytime runs have been used for the present analysis. The surface boundary-layer fluxes have been estimated using the eddy correlation method. The surface roughnessz ₀ obtained using the stability-corrected wind profiles (Paulson, 1970) has been found to be low for the Southwest monsson season. For the other season,z ₀ is relatively high. The drag coefficientC _D varies with height in the NE monsoon season but not in the season with lowz ₀. This aspect is reflected in the wind profiles for the two seasons and is discussed in detail. The scaling behaviour of friction velocityu _* and the turbulence intensity of longitudinal, lateral and vertical winds _u, _v and _w, respectively) are further examined to study their dependence on fetch. Our study shows that for the non-dimensional case, _u/u_* and _v/u_* do not show any surface roughness dependence in either season. On the other hand, for _w/u_* for the season with lowz ₀, the values are seen to agree well with that of Panofskyet al. (1977) for homogeneous terrain whereas for the other season with highz ₀, the results seem to conform more to the values observed by Smedman and Högström (1983) for coastal terrain. The results are discussed in the light of observations by other investigators.     

An isentropic vertical coordinate model: Design and application to atmospheric frontogenesis studies

Summary The isentropic vertical coordinate model developed at UCLA is briefly reviewed. The review includes an outline of the approach used to overcome technical difficulties in handling model layers with small mass.The model's performance is demonstrated by simulating the evolution of a middle-latitude baroclinic disturbance. During the evolution of the disturbance, sharp frontal zones are generated in the upper and middle troposphere with realistic tropopause folding. The extent to which different dynamical processes contribute to frontogenesis is analyzed.While the model successfully simulates frontogenesis in the upper and middle troposphere, it has a difficulty in simulating surface fronts. The difficulty arises due to the lack of degrees of freedom in surface temperatures since an isentropic vertical coordinate model requires a large number ofvertical layers to obtain a highhorizontal resolution at the lower boundary. This suggests the potential of a hybrid vertical coordinate, which approaches at upper levels and at lower levels.With 12 Figures     

Diurnal variation of horizontal wind direction fluctuations in complex terrain at Geysers,Cal.

Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km² Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s^-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.     

Boundary-layer parameterizations for applied dispersion modeling over urban areas

Many applied dispersion models require the knowledge of boundary-layer parameters such as sensible heat flux,Q _H , friction velocity,u _*, and turbulent energy components, _w and _v . Formulas are suggested for calculating these parameters over a wide variety of types of ground surfaces, based on simple observations of wind speed near the ground and fractional cloud cover, and specification of constants such as roughness length, albedo, and soil moisture availability. Observations ofu _*,Q _H , _w , and _v during field experiments in St. Louis and Indianapolis are used to test the formulas for urban sites. Relative errors of about ±20% in the predictions are seen to occur whenu _*,Q _H , _w , and _v are large. However, when these quantities are small (e.g.,u _* < 0.2 m/s), the errors in the predictions are as large as the mean value of the quantity itself.In addition, it is concluded from studies of available field data and theories that the magnitude of _w is not well-known at elevations above about 100m during the late afternoon and night. Some simple parameterizations for _w . are suggested that are consistent with the observed steady decrease in ground-level concentration in the afternoon and the sudden increase in concentration that can occur a few hours after sunset due to wind shears associated with a low-level jet, for continuous plumes emitted from moderate to tall stacks.     

Laboratory Simulation Of Vertical Plume Dispersion Within A Convective Boundary Layer

Mesoscale measurements of the vertical dispersion coefficient ₂ by using a composite turbulence water tank were validated through a comparison with CONDORS (Convective Diffusion Observed with Remote Sensors) field data, and were analysed with respect to the intensity of the thermal flux, mechanical turbulence, and plume release height.It seems possible to correct the plume _z values for different release heights below 0.5z_i (z_i is the mixing height) by applying an equation expressing the height dependency of turbulence intensity. The downwind distance where the plume's mass centre height approaches its final level was also analysed with respect to the above three parameters, and an empirical equation to estimate the downwinddistance derived.     

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     

High resolution spectral features of a series of aromatic hydrocarbons and BrO: Potential interferences in atmospheric OH-measurements

Naphthalene (C₁₀H₈), several other hydrocarbons, mostly derivates of naphthalene, and bromine oxide (BrO) were analyzed for narrow band (0.01 nm) absorption lines in the wavelength range between 307.7 and 308.3 nm to study their potential impact on OH radical measurements by differential absorption spectroscopy.Only naphthalene showed narrow band absorption lines in this wavelength region. From nine naphthalene lines the differential absorption cross-section was determined.The strongest naphthalene line at 308.002 nm is close to the Q ₁(2) OH line, but about a factor of 200 weaker (=(65.2±15.3)×10^-20 cm²/molec). The corresponding detection limit for naphthalene is about 15 ppt. We re-evaluated some spectra of our OH measurement campaign in July 1987 with respect to naphthalene and obtained an upper limit of 30 ppt for its concentration.BrO was recorded in the larger wavelength interval between 307.7 and 308.7 nm. Structured absorptions were only observed at wavelengths above 308.2 nm and no significant structures were found in the vicinity of the Q ₁(2) and Q ₁(3) OH lines.     

Markov chain simulations of vertical dispersion in the neutral surface layer for surface and elevated releases

Vertical dispersion in the neutral surface layer is investigated using a Markov Chain simulation procedure. The conceptual basis of the procedure is discussed and computation procedures outlined. Wind and turbulence parameterizations appropriate to the neutral surface layer are considered with emphasis on the Lagrangian time scale. Computations for a surface release are compared with field data. Good agreement is found for the variation of surface concentration and cloud height to distances 500 m downwind of the source. The functional form of the vertical concentration profile is examined and an exponential with exponent 1.6 is found to give the best fit with simulations.For elevated releases, it is demonstrated that an initial dip of the mass mean height from the simulation can be normalized for various release heights using a non-dimensionalized downwind coordinate incorporating advective wind speed and wind shear. The vertical distribution standard deviation ( _z ), as employed in Gaussian models, shows a fair degree of independence with source height but close examination reveals an optimum source height for maximum _z at a given downwind distance,x. This source height increases with downwind distance. Also the simulations indicate that vertical wind shear is more important than vertical variation of Lagrangian time scale close to the source, with a reverse effect farther downwind.     

Three-dimensional modeling of synthetic cold fronts interacting with northern Alpine foehn

Summary The effect of the Alpine orography on prototype cold fronts approaching from the west is investigated by three-dimensional numerical model simulations. The numerical experiments cover a range of parameter constellations which govern the prefrontal environment of the front. Especially, the appearance and intensity of prefrontal northern Alpine foehn varies from case to case.The behaviour of a cold front north of the Alps depends much on the prefrontal condition it encounters. It is found that prefrontal foehn can either accelerate or retard the approaching front.An important feature is the pressure depression along the northern Alpine rim that results from the southerly foehn flow. In cases where this depression compensates the eastward directed pressure gradient associated with the largescale flow, the front tends to accelerate and the foehn breaks down as soon as the front passes. In contrast, the foehn prevents the front from a rapid eastward propagation if it is connected with a strong southerly wind component.No-foehn experiments are performed for comparison, where either the mountains are removed, or the static stability is set to neutral. Also shown are effects of different crossfrontal temperature contrasts.List of Symbols c _F propagation speed of a front - x, y horizontal grid spacing (cartesian system) - , horizontal grid spacing (geographic system) - t time step - z vertical grid spacing (cartesian system) - cross-frontal potential temperature difference - _i potential temperature step at an inversion - E turbulent kinetic energy - f Coriolis parameter - FGP frontogenesis parameter (see section 2.2) - g gravity acceleration (g=9.81 m s^–2) - vertical gradient of potential temperature - h terrain elevation (above MSL) - h _i height of an inversion (h _i =1000 m MSL) - H height of model lid (H=9000 m MSL) - K _M exchange coefficient of momentum - K _H exchange coefficient of heat and moisture - longitude - N Brunt-Väisäla-frequency - p pressure - Exner function (=T/) - latitude - q _v specific humidity - R _d gas constant of dry air (R _d =287.06 J kg^–1 K^–1) - density of dry air - t time - T temperature - potential temperature - TFP thermal front parameter (see section 2.2) - u, v, w cartesian wind components - u _g ,v _g geostrophic wind components - horizontal wind vector - x, y, z cartesian coordinates Abbreviations GND (above) ground level - MSL (above) mean sea level - UTC universal time coordinated With 20 Figures     

Statistical properties of 49 years of rainfall rate events

Summary The paper presents a detailed investigation of magnitudes and properties pertaining to the large population of rainfall events recorded during 49 years at the Fabra Observatory in Barcelona. The study includes a statistical analysis of event durationT and rainfall quantityQ together with the statistical rainfall rate parameters: , ²(R) and maximum Rm within an event. The decorrelation time is also analysed. It is found thatQ, T, and can be well modelled by a lognormal distribution, but , ² and Rm are only so for a limited range of precipitation ratesR. The regression analysis between pairs of logarithms of the magnitudes investigated is generally good and a regression coefficient is often better than 0.9. Comparison with published work is also carried out. An attempt is made to discriminate between heavy and non-heavy rainfall rates, and the 50 mm/h threshold is used for the study. The twelve-monthly running average indicates that the rainfall amount has a small increasing trend over the fifty year period. However, this trend is reversed when considering heavy rains. Finally, the return period in years to exceeding a thresholdR within an event is also investigated and the distribution of the population of annual extremes is found to be Gumbel II.With 6 Figures