Studies of the atmospheric stability characteristics during the solar eclipse of February 16, 1980

Observations made with a monostatic sodar and from a 120 m instrumented tower have been used to study the variations in the atmospheric boundary layer at Tarapur (19° 50 N, 72° 41 E) during the solar eclipse of February 16, 1980. Atmospheric instability was reduced below normal values during the eclipse but the atmosphere at no time became stable.     

Measurement of atmospheric total ozone by the filter photometric method

The total ozone content in the atmosphere was determined from the multichannel photometer observations of direct solar radiation made in the urban environment at Pune (18° 32 N, 73° 51E, 559 m ASL) and Sinhagad hill station (18° 22N, 73° 45E, 1305 m ASL) during March 1980-February 1982. The total ozone content of the atmosphere was computed making use of the differential absorption of solar radiation due to ozone at 0.4 and 0.6 m wavelengths in the Chappuis band. The values of the ozone data obtained from the photometer observations at Pune and Sinhagad were compared with the corresponding ozone data obtained from the Dobson spectrophotometer located at Pune. Values of ozone obtained by the photometric method were found to be smaller by 8–18% than the Dobson values when Vigroux's absorption coefficients were used. Similarly, when the absorption coefficients of Inn and Tanaka (1953) were used, the ozone values obtained by the photometric method were smaller by 4–14% than the Dobson values. The ozone values at the hill station obtained from the photometric method were in better agreement (5%) with the Dobson values.     

Maximum precipitation for different short-term intervals

Summary A time series and extreme value analysis of maximum precipitation for distinct time intervals from 10 minutes to 1 day have been deduced for the observatory Zagreb-Gri ( = 45°49 and = 15°59,H _s = 157 m) for the period between 1908 and 1985. The Spearman rank correlation test revealed that the short-term precipitation maxima series for the time intervals of 10 to 30 minutes and 8 to 24 hours indicate no significant trend for the 95 percent probability level. The 40-minute to 4-hour precipitation maxima exhibited an increaese in recent time. Annual maximumt-minute precipitation was estimated using the Gumbel distribution for sets of data gradually prolonged with 10-year steps towards the past. The ratios (q) of maximum precipitation estimates for the shorter periods and of estimates for the reference period (1908–1985) have been calculated. Their positions according to the confidence interval on the reference estimates were determined. Confident extreme value estimates could be obtained using at least the 50–60 year data series.With 3 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.     

The application of a two-dimensional numerical model of the planetary boundary layer to the Nanticoke region on the north shore of Lake Erie

A two-dimensional atmospheric boundary-layer model is applied to the Nanticoke region on the northern shore of Lake Erie (80 ° 03W and 42 ° 50N) to simulate numerically the observed wind and temperature profiles. In general, the profiles predicted by the model agree reasonably well with the observed profiles.     

High-Latitude Springtime Photochemistry. Part I: NOx, PAN and Ozone Relationships

Springtime measurements of NO_x, ozone, PAN,J(NO₂), and other compounds were made near Ny-Ålesund,Svalbard (78°54N, 11°53E), in 1994 and Poker Flat,Alaska (65°08N, 147°29W), in 1995. At Svalbard medianmixing ratios for PAN and NO_x of 237 and 23.7 pptv,respectively, were observed. The median mixing ratios at Poker Flat for PANand NO_x were 79.5 and 85.9 pptv, respectively. These data areused to estimate thermal PAN decomposition using several differentapproaches. At Svalbard PAN decomposition was very small, while at PokerFlat up to 30 pptv/h PAN decomposed. At both sites the NO_x/PANratio increased with temperature between –10 and 20°C implyingthat PAN decomposition is an important NO_x source. In-situozone production was calculated from the measured NO, NO₂,O₃, J(NO₂), and temperature data, using thesteady state assumption Median ozone production was 605 pptv/h at PokerFlat, and one order of magnitude smaller at Svalbard during the daytime.Only at Poker Flat could a direct influence on the diurnal ozone cycle beobserved from in-situ production. These results imply that PAN decompositionis a major source of NO_x in the high latitude troposphere, andthat this contributes to the observed spring maximum in surface ozone.     

Sea breeze activity at an inland station Kharagpur (India) — A case study

Surfaces fluxes, turbulent kinetic energy and Flux Richardson number are calculated for three typical sea breeze days characterizing three types of sea breeze onset at an inland station Kharagpur (22°21 N, 87°19 E), 80 km inland, one of the sites for MONTBLEX (MONsoon Trough Boundary Layer EXperiment), in India. The sea breeze onset is associated with a decrease in momentum and heat fluxes and an increase in moisture flux. The turbulent kinetic energy shows quite a significant value even in the late afternoon. The surface layer becomes nearly stable even before sunset, due to the passage of the sea breeze.     

Variations in atmospheric CO2 at the Meteorological Research Institute, Tsukuba, Japan

Since April 1986, measurements of the CO₂ concentration in the surface air have been conducted at the Meteorological Research Institure (MRI, 36°04 N, 140°07 E, 25 m above sea level) in Tsukuba, located 50 km northeast of Tokyo, Japan. The CO₂ data measured over times between 11:00 Japan Standard Time (JST) and 16:00 JST (C _N ) were considered to be representative of the air (within a few ppmv) in the planetary boundary layer. To evaluate the representative CO₂ level on a spatial scale larger than that of the C _N record, the CO₂ data with hour-to-hour variation less than 1 ppmv were selected (C _P ). Comparison of these data with those of Ryori (39°02 N, 141°50 E), a continental station operated by the Japan Meteorological Agency, indicates that the C _P record provides a representative CO₂ level in the air on spatial scales of at least a few hundred kilometers.The C _N record allows an investigation of the internanual changes in photosynthesis/respiration against changes in climatological parameters. Within a small temperature anomaly (ca.±1 °C) respiration is sensitive to the temperature change, while photosynthesis is less sensitive. When the temperature anomaly is large, however, photosynthesis and respiration tend to be competitive.     

Trajectory Curvature As A Selection Criterion For valid Lagrangian Stochastic Dispersion Models

Among well-mixed multi-dimensional Lagrangian stochastic (LS) dispersion models, we observe that those in poorest agreement with observations produce spiralling trajectories, with an associated reduction in dispersion. We therefore investigate statistics of increments d ' to the orientation '= arctan(W'/U') of the Lagrangian velocity-fluctuation vector – as a possible means to distinguish the better LS models within the well-mixed class. Zero-spin models, having d' = 0, are found to provide best agreement with observations. It is not clear however, whether imposition of the zero-spin property selects (in conjunction with the well-mixed condition) a unique model.     

Calculation of velocity skewness in real and artificial plant canopies

Wind velocities within a plant canopy are much more strongly skewed than those of the air flow above. We have examined the governing Eulerian equations for the velocity products u _i, u _j u_k using data from a wind tunnel study with an artificial canopy consisting of an array of 5 cm lengths of monofilament fishing line, and from measurements in corn (Zea mays L).Simple parameterizations for pressure-velocity correlations, and for the quadruple velocity products allowed reasonably accurate calculations of the third moments using measured profiles of the mean velocity, variance and covariance fields. Comparisons of individual terms in the rate equations for ovu _i, u _j u _krevealed that diffusion (from above) and mean shear were most important in creating large skewness in the canopy. A drag term also contributed but was of lesser importance. These terms were balanced by return-to-isotropy and a turbulence interaction term. A quasi-Gaussian approximation considerably underestimated the magnitude of the fourth moments within the canopy.     

Cold surges along the front range of the Rocky mountains: Development of a classification scheme

Summary A classification scheme for cold surges occurring along the Front Range of the Rocky Mountains has been developed using NMC surface and upper air data. A separation is seen in the data between cases with parent anticyclones that moved from Canada south into the United States (Northern type) and those which moved from the Eastern Pacific Ocean east into the central US (Western type). The Northern surges are further broken down into three types, A, B and B, when the degree of southward motion of the anticyclone and the leading edge of the cold air is considered. An anticyclone movement south of 40°N distinguishes type A, while anticyclone tracks remaining north of 40°N characterize types B and B. The leading edge of the surge (the surge line) associated with the type A anticyclone travels south of 25N, while the type B anticyclones produce two distinct surge line movements, one remaining north of 25°N and the other moving south of 25°N. It is not possible to break down the Western surges into distinct types. As a rule, the anticyclones associated with Western surges remain between 50°N and 30°N, with the surge line moving south to 18°N on average.Though the classification scheme is developed in a relatively subjective manner, the statistical robustness of the scheme is verified through application of the nonparametric Mann-Whitney test. Results of the test clearly provide support for the division of the surges that was obtained.With 16 Figures     

On the enhanced smoothing over topography in some mesometeorological models

An equation is derived for the components of the horizontal (turbulent) frictional force in the -coordinate system with special attention to mesometeorological flow models. The starting point is the horizontal equation of motion in its flux-form in the -system in which we replace (following Reynolds' procedure) the velocity components u,v and % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbai % aaaaa!37B8! \[ \dot \sigma \] aswell as other relevant quantities by terms of the form u = + u,..., = ± + % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbai % Gbauaaaaa!37C3! \[ \dot \sigma ' \] , etc. ( = time average of u; u = fluctuating part of u.) Next, the equation is averaged with respect to time and terms which we believe are small in mesometeorological flows, are neglected. On expressing by an appropriate expression that involves w, the result shows the appearance of two new terms which, have not been considered previously in the published literature. While the expression earlier used in the literature involved the -derivative of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG3bGbauaaaaaaaa!380B!\[\overline {u'w'} \] alone, the new terms add the -derivatives of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaadaahaaWcbeqaaiaaikdaaaaaaaaa!37EC!\[\overline {u'^2 } \] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG2bGbauaaaaaaaa!380A!\[\overline {u'v'} \] for the x-component of the force, and the -derivatives of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG2bGbauaadaahaaWcbeqaaiaaikdaaaaaaaaa!37ED!\[\overline {v'^2 } \]} and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG2bGbauaaaaaaaa!380A!\[\overline {u'v'} \] for the y-component, where and are the slopes of the -surfaces in the x- and y-directions, respectively. Further, a few numerical simulations of the sea-breeze over topography are carried out with and without the correction terms. It is shown that when corrections terms are not included the effective smoothing is stronger above the sloping regions and may amount to as high as 50 percent of the convergence with slopes of ~.04. The ìnclusìon of the new terms does not lead to any special computational difficulties and for that reason there is no compelling reason to neglect them, all the more so because, as is shown, the addition of the new terms results in a consistent apportioning of the degree of horizontal diffusion.On leave from CIMMS, Norman, OK.Now visiting Dept. of Met., Helsinki, Finland.     

Time series analyses of concentration and wind fluctuations

Analyses of concentration fluctuation (C) spectra from boundary-layer smoke plume experiments at six separate locations show that the spectra from these experiments generally exhibit an inertial subrange at high frequencies with a slope of -5/3 and indicate peak energy at a time period of about 50 to 100 s. These periods of peak energy are a factor of two to five less than those for the peak of the wind speed fluctuation (u or v) spectra. A general spectral formula fits normalized spectra from the U.S. and Australia, where the frequency, n, is made dimensionless by multiplying by the plume dispersion parameter, _y , and dividing by the wind speed, u. Peak energy occurs at a dimensionless frequency of n _y/u equal to about 0.15. The Kolmogorov constant in the inertial subrange is estimated from a set of averaged spectra. Cross-spectra indicate little relation between concentration and wind fluctuations. However, most of the correlation that exists is due to periods larger than about 10 or 20 s.     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.     

The Value Of The β Coefficient In The Relaxed Eddy Accumulation Method In Terms Of Fourth-Order Moments

The commonly measured value of in the relaxed eddy accumulationmethod of about 0.56is shown to arise from the non-Gaussiannature of turbulence. Fourth-orderGram–Charlier functions forthe two-dimensional probability distributionsof variation in the horizontal component of wind velocityand concentrations of water vapour, carbondioxide and methane with respect to thevertical component of wind velocity are used to examinethe value of .An analytical solution for ispresented in terms of fourth-order moments.Under mean conditions, this solution givesa value for of0.557. Variation of is shown to be controlledprimarily by the ratio of the mean ofc'w3 (where c'is relevant to the entity of interest andw' is vertical component of windvelocity) to the correlationcoefficient between the entity concentrationand vertical component of wind velocity.     

Seasonal variation of atmospheric dimethylsulfide at Amsterdam Island in the southern Indian Ocean

Daily measurements of atmospheric concentrations of dimethylsulfide (DMS) were carried out for two years in a marine site at remote area: the Amsterdam Island (37°50S–77°31E) located in the southern Indian Ocean. DMS concentrations were also measured in seawater. A seasonal variation is observed for both DMS in the atmosphere and in the sea-surface. The monthly averages of DMS concentrations in the surface coastal seawater and in the atmosphere ranged, respectively, from 0.3 to 2.0 nmol l^-1 and from 1.4 to 11.3 nmol m^-3 (34 to 274 pptv), with the highest values in summer. The monthly variation of sea-to-air flux of DMS from the southern Indian Ocean ranges from 0.7 to 4.4 mol m^-2 d^-1. A factor of 2.3 is observed between summer and winter with mean DMS fluxes of 3.0 and 1.3 mol m^-2 d^-1, respectively.     

Wind characteristics in the lowest 340 m of the atmospheric boundary layer at Pune,India

Wind characteristics in the lowest 340 m (agl) of the atmospheric boundary layer at Pune (18°32N, 73°51E, 559 m ASL) have been investigated using the pilot balloon wind observations obtained during the summer monsoon seasons of 1976, 1979 and 1980. Variations in the zonal and meridional components of wind at the surface, 40, 150 and 340 m (agl) have been described. Frequency distribution of the zonal component showed an unimodal character at the surface, which gradually approached a multimodal character at 340 m. The distribution pattern of the wind components was, by and large, normal. Spectral analysis of the wind components showed that the spectral energy was predominantly shared by 5–7 and 2–3 day periods.     

Glaciers and Climate in the Andes between the Equator and 30° S: What is Recorded under Extreme Environmental Conditions?

Sublimation and melt disturb the environmental information obtained from ice core records in the Andes. In two case studies we demonstrate to what extent these post-depositional processes may remove major parts of the accumulated snow cover. Dark ash layers from the Tungurahua eruption changed the albedo of surface snow on Chimborazo glacier (6268 m, 1°30 S,78°36 W, Ecuador) between two ice core drilling campaigns and forced substantial melt. Re-distribution and washout of the chemical constituents shifted the concentration profiles obtained in December 1999 as compared to an equivalent core drilled in December 2000. The stable isotope records showed that approximately the water equivalent (weq) of an annual layer had melted, and that the percolating melt water penetrated within the firn layer to a depth of at least 16.5 m without refreezing. In the second example, from a site on the dry axis between the tropical and extra-tropical precipitation belts, significant loss of accumulated snow layers occurred by sublimation. A surface experiment at Cerro Tapado glacier (5536 m, 30°08 S,69°55 W, Chile) revealed that losses of 2 mm weq (5 mm snow) per day occurred during the dry period following the 1997/98 El Niño. This loss generally included the entire surface layer enriched in stable isotopes, and thus caused minimal disturbance of the isotopic signature (and hence climatic information) of the net accumulation, yet chemical constituents again experienced considerable changes in concentration. From annual layer counting and direct dating it is obvious that the major part of the accumulated ice on both glaciers is younger than 100 years; however, isotopic and chemical variations at least in the basal ice from Cerro Tapado clearly point to climate conditions different from the recent centuries. This evidence is supported by mass balance considerations derived from a glacier-climate model. The possibility of a third type of disturbance aside from sublimation and melting – in this case a significant hiatus in the environmental chronology – also deserves consideration for other icecore records from this region. Potential disruptions or discontinuities need to be carefully evaluated given the profound changes in climatic and glaciological conditions since the Last Glacial Maximum throughout Holocene times.     

The Henry's Law Constants of the Haloacetic Acids

Henry's law constants K_H (mol kg^-1 atm^-1) have been measured between 278.15 K and 308.15 K for the following organic acids: CH₂FCOOH (ln(K_H[298.15 K]) = 11.3 ± 0.2), CH₂ClCOOH (11.59 ± 0.14), CH₂BrCOOH (11.94 ± 0.21), CHF₂COOH (10.32 ± 0.10), CHCl₂COOH (11.69 ± 0.11), CHBr₂COOH (12.33 ± 0.29), CBr₃COOH (12.61 ± 0.21), and CClF₂COOH (10.11 ± 0.12). The variation of K_H with temperature was determined for all acids except CH₂FCOOH and CBr₃COOH, with _r H° for the dissolution reaction ranging from –85.2 ± 2.6 to –57.1 ± 2.5 kJ mol^-1, meaning that their solubility is generally more sensitive to temperature than is the case for the simple carboxylic acids. The Henry's law constants show consistent trends with halogen substitution and, together with their high solubility compared to the parent (acetic) acid (ln(K_H[298.15 K]) = 8.61), present a severe test of current predictive models based upon molecular structure. The solubility of haloacetic acids and strong dissociation at normal pH mean that they will partition almost entirely into cloud and fog in the atmosphere (0.05–1.0 g H₂O m^-3), but can reside in both phases for the liquid water contents typical of aerosols (10^-5-10^-4 g H₂O m^-3).     

Observed variations of σθ and σΦ in the nocturnal drainage flow in a deep valley

The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms^-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and