Search NASA⌕ Search

SEARCH · Search NASA

Results for “fall speed”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Observations of Precipitation Size and Fall Speed Characteristics within Coexisting Rain and Wet Snow

Ground-based measurements of particle size and fall speed distributions using a Particle Size and Velocity (PARSIVEL) disdrometer are compa red among samples obtained in mixed precipitation (rain and wet snow) and rain in the Oregon Cascade Mountains and in dry snow in the Rock y Mountains of Colorado. Coexisting rain and snow particles are distinguished using a classification method based on their size and fall sp eed properties. The bimodal distribution of the particles' joint fall speed-size characteristics at air temperatures from 0.5 to 0 C suggests that wet-snow particles quickly make a transition to rain once mel ting has progressed sufficiently. As air temperatures increase to 1.5 C, the reduction in the number of very large aggregates with a diame ter > 10 mm coincides with the appearance of rain particles larger than 6 mm. In this setting. very large raindrops appear to be the result of aggregates melting with minimal breakup rather than formation by c oalescence. In contrast to dry snow and rain, the fall speed for wet snow has a much weaker correlation between increasing size and increasing fall speed. Wet snow has a larger standard deviation of fall spee d (120%-230% relative to dry snow) for a given particle size. The ave rage fall speed for observed wet-snow particles with a diameter great er than or equal to 2.4 mm is 2 m/s with a standard deviation of 0.8 m/s. The large standard deviation is likely related to the coexistence of particles of similar physical size with different percentages of melting. These results suggest that different particle sizes are not required for aggregation since wet-snow particles of the same size can have different fall speeds. Given the large standard deviation of fa ll speeds in wet snow, the collision efficiency for wet snow is likely larger than that of dry snow. For particle sizes between 1 and 10 mm in diameter within mixed precipitation, rain constituted I % of the particles by volume within the isothermal layer at 0 C and 4% of the particles by volume for the region just below the isothermal layer where air temperatures rise from 0" to 0.5"C. As air temperatures increa sed above 0.5 C, the relative proportions of rain versus snow particl es shift dramatically and raindrops become dominant. The value of 0.5 C for the sharp transition in volume fraction from snow to rain is sl ightly lower than the range from 1 .l to 1.7 C often used in hydrolog ical models.

Yuter, Sandra E.↗

Dual measurement terminal fall speeds and multiple Doppler winds

It is shown that radar-derived terminal fall speed measurements can be useful in determining vertical air velocity in the middle troposphere by means of a network of Doppler radars. The theoretical principles of the dual measurement technique are described, and the relationship between measurement accuracies and theoretical estimates of terminal fall speeds is discussed. It is demonstrated that the use of differential reflectivity to estimate terminal fall speeds can reduce the standard error of vertical velocity estimates by 40-50 percent.

Grosh, R. C.↗

Evaluation of the New Version of the Laser-Optical Disdrometer, OTT Parsivel2

A comparative study of raindrop size distribution measurements has been conducted at NASA's Goddard Space Flight Center where the focus was to evaluate the performance of the upgraded laser-optical OTT Particle Size Velocity (Parsivel2; P2) disdrometer. The experimental setup included a collocated pair of tipping-bucket rain gauges, OTT Parsivel (P1) and P2 disdrometers, and Joss-Waldvogel (JW) disdrometers. Excellent agreement between the two collocated rain gauges enabled their use as a relative reference for event rain totals. A comparison of event total showed that the P2 had a 6%absolute bias with respect to the reference gauges, considerably lower than the P1 and JW disdrometers. Good agreement was also evident between the JW and P2 in hourly raindrop spectra for drop diameters between 0.5 and 4 mm. The P2 drop concentrations mostly increased toward small sizes, and the peak concentrations were mostly observed in the first three measurable size bins. The P1, on the other hand, underestimated small drops and overestimated the large drops, particularly in heavy rain rates. From the analysis performed, it appears that the P2 is an improvement over the P1 model for both drop size and rainfall measurements. P2 mean fall velocities follow accepted terminal fall speed relationships at drop sizes less than 1 mm. As a caveat, the P2 had approximately 1ms21 slower mean fall speed with respect to the terminal fall speed near 1 mm, and the difference between the mean measured and terminal fall speeds reduced with increasing drop size. This caveat was recognized as a software bug by the manufacturer and is currently being investigated.

disdrometers↗

On the Measurement Errors of the Joss-Waldvogel Disdrometer

The Joss-Waldvogel (JW) disdrometer is considered to be a reference instrument for drop size distribution measurements. It has been widely used in many field campaigns as part of validation efforts of radar rainfall estimation. It has also been incorporated in radar rain gauge rainfall observation networks at several ground validation sites for NASA s Tropical Rainfall Measuring Mission (TRMM). It is anticipated that the Joss-Waldvogel disdrometer will be one of the key instruments for ground validation for the upcoming Global Precipitation Measurement (GPM) mission. The JW is an impact type disdrometer and has several shortcomings. One such shortcoming is that it underestimates the number of small drops in heavy rain due to the disdrometer dead time. The detection of smaller drops is also suppressed in the presence of background noise. Further, drops larger than 5.0 to 5.5 mm diameter cannot be distinguished by the disdrometer. The JW assumes that all raindrops fall at their terminal fall speed. Ignoring the influence of vertical air motion on raindrop fall speed results in errors in determining the raindrop size. Also, the bulk descriptors of rainfall that requires the fall speed of the drops will be overestimated or underestimated due to errors in measured size and assumed fall velocity. Long-term observations from a two-dimensional video disdrometer are employed to simulate the JW disdrometer and assess how it s shortcomings affect radar rainfall estimation. Data collected from collocated JW disdrometers were also incorporated in this study.

Tokay, Ali↗

Comparison of Raindrop Size Distribution Measurements by Collocated Disdrometers

An impact-type Joss-Waldvogel disdrometer (JWD), a two-dimensional video disdrometer (2DVD), and a laser optical OTT Particle Size and Velocity (PARSIVEL) disdrometer (PD) were used to measure the raindrop size distribution (DSD) over a 6-month period in Huntsville, Alabama. Comparisons indicate event rain totals for all three disdrometers that were in reasonable agreement with a reference rain gauge. In a relative sense, hourly composite DSDs revealed that the JWD was more sensitive to small drops (,1 mm), while the PD appeared to severely underestimate small drops less than 0.76mm in diameter. The JWD and 2DVD measured comparable number concentrations of midsize drops (1-3mm) and large drops (3-5 mm), while the PD tended to measure relatively higher drop concentrations at sizes larger than 2.44mm in diameter. This concentration disparity tended to occur when hourly rain rates and drop counts exceeded 2.5mm/h and 400/min, respectively. Based on interactions with the PD manufacturer, the partially inhomogeneous laser beam is considered the cause of the PD drop count overestimation. PD drop fall speeds followed the expected terminal fall speed relationship quite well, while the 2DVD occasionally measured slower drops for diameters larger than 2.4mm, coinciding with events where wind speeds were greater than 4m/s. The underestimation of small drops by the PD had a pronounced effect on the intercept and shape of parameters of gamma-fitted DSDs, while the overestimation of midsize and larger drops resulted in higher mean values for PD integral rain parameters

disdrometers↗

GCSS Idealized Cirrus Model Comparison Project

The GCSS Working Group on Cirrus Cloud Systems (WG2) is conducting a systematic comparison and evaluation of cirrus cloud models. This fundamental activity seeks to support the improvement of models used for climate simulation and numerical weather prediction through assessment and improvement of the "process" models underlying parametric treatments of cirrus cloud processes in large-scale models. The WG2 Idealized Cirrus Model Comparison Project is an initial comparison of cirrus cloud simulations by a variety of cloud models for a series of idealized situations with relatively simple initial conditions and forcing. The models (16) represent the state-of-the-art and include 3-dimensional large eddy simulation (LES) models, two-dimensional cloud resolving models (CRMs), and single column model (SCM) versions of GCMs. The model microphysical components are similarly varied, ranging from single-moment bulk (relative humidity) schemes to fully size-resolved (bin) treatments where ice crystal growth is explicitly calculated. Radiative processes are included in the physics package of each model. The baseline simulations include "warm" and "cold" cirrus cases where cloud top initially occurs at about -47C and -66C, respectively. All simulations are for nighttime conditions (no solar radiation) where the cloud is generated in an ice supersaturated layer, about 1 km in depth, with an ice pseudoadiabatic thermal stratification (neutral). Continuing cloud formation is forced via an imposed diabatic cooling representing a 3 cm/s uplift over a 4-hour time span followed by a 2-hour dissipation stage with no cooling. Variations of these baseline cases include no-radiation and stable-thermal-stratification cases. Preliminary results indicated the great importance of ice crystal fallout in determining even the gross cloud characteristics, such as average vertically-integrated ice water path (IWP). Significant inter-model differences were found. Ice water fall speed is directly related to the shape of the particle size distribution and the habits of the ice crystal population, whether assumed or explicitly calculated. In order to isolate the fall speed effect from that of the associated ice crystal population, simulations were also performed where ice water fall speed was set to the same constant value everywhere in each model. Values of 20 and 60 cm/s were assumed. Current results of the project will be described and implications will be drawn. In particular, this exercise is found to strongly focus the definition of issues resulting in observed inter-model differences and to suggest possible strategies for observational validation of the models. The next step in this project is to perform similar comparisons for well observed case studies with sufficient high quality data to adequately define model initiation and forcing specifications and to support quantitative validation of the results.

Starr, David OC.↗

Comparing Aircraft Observations of Snowfall to Forecasts Using Single or Two Moment Bulk Water Microphysics Schemes

High resolution weather forecast models with explicit prediction of hydrometeor type, size distribution, and fall speed may be useful in the development of precipitation retrievals, by providing representative characteristics of frozen hydrometeors. Several single or double-moment microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, allowing for the prediction of up to three ice species. Each scheme incorporates different assumptions regarding the characteristics of their ice classes, particularly in terms of size distribution, density, and fall speed. In addition to the prediction of hydrometeor content, these schemes must accurately represent the vertical profile of water vapor to account for possible attenuation, along with the size distribution, density, and shape characteristics of ice crystals that are relevant to microwave scattering. An evaluation of a particular scheme requires the availability of field campaign measurements. The Canadian CloudSat/CALIPSO Validation Project (C3VP) obtained measurements of ice crystal shapes, size distributions, fall speeds, and precipitation during several intensive observation periods. In this study, C3VP observations obtained during the 22 January 2007 synoptic-scale snowfall event are compared against WRF model output, based upon forecasts using four single-moment and two double-moment schemes available as of version 3.1. Schemes are compared against aircraft observations by examining differences in size distribution, density, and content. In addition to direct measurements from aircraft probes, simulated precipitation can also be converted to equivalent, remotely sensed characteristics through the use of the NASA Goddard Satellite Data Simulator Unit. Outputs from high resolution forecasts are compared against radar and satellite observations emphasizing differences in assumed crystal shape and size distribution characteristics.

Molthan, Andrew L.↗

Use of 2d-video Disdrometer to Derive Mean Density-size and Ze-SR Relations: Four Snow Cases from the Light Precipitation Validation Experiment

The application of the 2D-video disdrometer to measure fall speed and snow size distribution and to derive liquid equivalent snow rate, mean density-size and reflectivity-snow rate power law is described. Inversion of the methodology proposed by Böhm provides the pathway to use measured fall speed, area ratio and '3D' size measurement to estimate the mass of each particle. Four snow cases from the Light Precipitation Validation Experiment are analyzed with supporting data from other instruments such as Precipitation Occurrence Sensor System (POSS), Snow Video Imager (SVI), a network of seven snow gauges and three scanning C9 band radars. The radar-based snow accumulations using the 2DVD-derived Ze-SR relation are in good agreement with a network of seven snow gauges and outperform the accumulations derived from a climatological Ze-SR relation used by the Finnish Meteorological Institute (FMI). The normalized bias between radar-derived and gauge accumulation is reduced from 96% when using the fixed FMI relation to 28% when using the Ze-SR relations based on 2DVD data. The normalized standard error is also reduced significantly from 66% to 31%. For two of the days with widely different coefficients of the Ze-SR power law, the reflectivity structure showed significant differences in spatial variability. Liquid water path estimates from radiometric data also showed significant differences between the two cases. Examination of SVI particle images at the measurement site corroborated these differences in terms of unrimed versus rimed snow particles. The findings reported herein support the application of Böhm's methodology for deriving the mean density-size and Ze-SR power laws using data from 2D-video disdrometer.

radar↗

Intercomparison of Large-Eddy Simulations of Arctic Mixed-Phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds.

clouds↗

Derivation of Physical and Optical Properties of Midlatitude Cirrus Ice Crystals for a Size-Resolved Cloud Microphysics Model

Single-crystal images collected in mid-latitude cirrus are analyzed to provide internally consistent ice physical and optical properties for a size-resolved cloud microphysics model, including single-particle mass, projected area, fall speed, capacitance, single-scattering albedo, and asymmetry parameter. Using measurements gathered during two flights through a widespread synoptic cirrus shield, bullet rosettes are found to be the dominant identifiable habit among ice crystals with maximum dimension (Dmax) greater than 100μm. Properties are therefore first derived for bullet rosettes based on measurements of arm lengths and widths, then for aggregates of bullet rosettes and for unclassified (irregular) crystals. Derived bullet rosette masses are substantially greater than reported in existing literature, whereas measured projected areas are similar or lesser, resulting in factors of 1.5-2 greater fall speeds, and, in the limit of large Dmax, near-infrared single-scattering albedo and asymmetry parameter (g) greater by approx. 0.2 and 0.05, respectively. A model that includes commonly imaged side plane growth on bullet rosettes exhibits relatively little difference in microphysical and optical properties aside from approx. 0:05 increase in mid-visible g primarily attributable to plate aspect ratio. In parcel simulations, ice size distribution, and g are sensitive to assumed ice properties.

asymmetry↗

Evaluating the Performance of Single and Double Moment Microphysics Schemes During a Synoptic-Scale Snowfall Event

Increases in computing resources have allowed for the utilization of high-resolution weather forecast models capable of resolving cloud microphysical and precipitation processes among varying numbers of hydrometeor categories. Several microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, ranging from single-moment predictions of precipitation content to double-moment predictions that include a prediction of particle number concentrations. Each scheme incorporates several assumptions related to the size distribution, shape, and fall speed relationships of ice crystals in order to simulate cold-cloud processes and resulting precipitation. Field campaign data offer a means of evaluating the assumptions present within each scheme. The Canadian CloudSat/CALIPSO Validation Project (C3VP) represented collaboration among the CloudSat, CALIPSO, and NASA Global Precipitation Measurement mission communities, to observe cold season precipitation processes relevant to forecast model evaluation and the eventual development of satellite retrievals of cloud properties and precipitation rates. During the C3VP campaign, widespread snowfall occurred on 22 January 2007, sampled by aircraft and surface instrumentation that provided particle size distributions, ice water content, and fall speed estimations along with traditional surface measurements of temperature and precipitation. In this study, four single-moment and two double-moment microphysics schemes were utilized to generate hypothetical WRF forecasts of the event, with C3VP data used in evaluation of their varying assumptions. Schemes that incorporate flexibility in size distribution parameters and density assumptions are shown to be preferable to fixed constants, and that a double-moment representation of the snow category may be beneficial when representing the effects of aggregation. These results may guide forecast centers in optimal configurations of their forecast models for winter weather and identify best practices present within these various schemes.

Molthan, Andrew L.↗

Multiyear Simulations of the Martian Water Cycle with the Ames General Circulation Model

Mars atmosphere is carbon dioxide dominated with non-negligible amounts of water vapor and suspended dust particles. The atmospheric dust plays an important role in the heating and cooling of the planet through absorption and emission of radiation. Small dust particles can potentially be carried to great altitudes and affect the temperatures there. Water vapor condensing onto the dust grains can affect the radiative properties of both, as well as their vertical extent. The condensation of water onto a dust grain will change the grain s fall speed and diminish the possibility of dust obtaining high altitudes. In this capacity, water becomes a controlling agent with regard to the vertical distribution of dust. Similarly, the atmosphere s water vapor holding capacity is affected by the amount of dust in the atmosphere. Dust is an excellent green house catalyst; it raises the temperature of the atmosphere, and thus, its water vapor holding capacity. There is, therefore, a potentially significant interplay between the Martian dust and water cycles. Previous research done using global, 3-D computer modeling to better understand the Martian atmosphere treat the dust and the water cycles as two separate and independent processes. The existing Ames numerical model will be employed to simulate the relationship between the Martian dust and water cycles by actually coupling the two cycles. Water will condense onto the dust, allowing the particle's radiative characteristics, fall speeds, and as a result, their vertical distribution to change. Data obtained from the Viking, Mars Pathfinder, and especially the Mars Global Surveyor missions will be used to determine the accuracy of the model results.

R M Haberle↗

Data Comparison: Satellite and Falling Sphere Temperatures

Small meteorological rocketsondes providing temperature data have beam used for comparison with, and validation of measurements from satellite-borne instruments. A significant number of rocket-borne falling spheres were launched in conjunction with the Upper Atmosphere Research Satellite (UARS) for validation of the Halogen Occultation Experiment (HALOE), High Resolution Doppler Interferometer (HRDI), and the Microwave Limb Sounder (MLS) instruments. Upper stratosphere and mesosphere temperatures measured with these instruments on UARS are compared with inflatable spheres launched from Wallops Island (1992-1999), Brazil (1994), Hawaii (1992), Norway (1992), and Sweden (1993 and 1996). Time and space differences varied between the satellite measurement and the rocketsonde launch, for example HALOE overpasses occurred within 5 days and in some cases there were spatial differences of up to 30 degrees longitude. Validation measurements of the HRDI instrument occurred at Wallops Island when it passed within 20 minutes and 330 kilometers of the launch site. Because of discontinuity in the falling sphere drag coefficients when fall speed neared MACH 1 falling sphere temperatures near 70 kilometers attitude are biased toward lower temperatures. Availability of improved software and a new atmospheric model have helped to reduce this bias. The validated remote instrument measurements permit a new perspective of atmospheric structure to be formed, not always possible with the limited number of falling sphere measurements. Features of the remote measurement temperature profiles and their possible use to extend the climatological data base at the rocketsonde sites will be discussed.

Schmidlin, Francis J.↗

The role of entrainment by falling raindrops in microbursts

The numerical model of Krueger et al. (1986) for dry microburst simulations is used to study the role of entrainment by falling raindrops. Two series of numerical simulations were conducted: a control series, and a series with the raindrop fall speed set to zero so that the rain moved with the air instead of falling through it. The results show that entrainment due to falling raindrops helps microbursts with large raindrop mixing ratios to form in stable stratifications. It is found that entrainment appears to contribute to the small spatial and temporal scales that characterize microburst outflows.

Krueger, Steven K.↗

Spaceborne Radar Measurements of Rainfall Vertical Velocity

This paper studies the performance of a spaceborne precipitation radar in measuring vertical Doppler velocity of rainfall. As far as a downward pointing precipitation radar is concerned, one of the major problems affecting Doppler measurement at the nadir direction arises from the Non-Uniform Beam-Filling effect (NUBF). That is, when significant variation in rain rate is present within the radar IFOV (Instrument Field of View) in the along track direction. the Doppler shift caused by the radial component of the horizontal speed of the satellite is weighted differently among the portions of IFOV. The effects of this non-uniform weighting may dominate any other contribution. Under this condition, shape, average value and width of the Doppler spectrum may not be directly correlated with the vertical velocity of the precipitating particles. However, by using an inversion technique which over-samples the radar measurements in the along track direction, we show that the shift due to NUBF can be evaluated, and that the NUBF induced errors on average fall speed can be reduced.

Im, Eastwood↗

Estimating the path-average rainwater content and updraft speed along a microwave link

There is a scarcity of methods for accurately estimating the mass of rainwater rather than its flux. A recently proposed technique uses the difference between the observed rates of attenuation A with increasing distance at 38 and 25 GHz, A(38-25), to estimate the rainwater content W. Unfortunately, this approach is still somewhat sensitive to the form of the drop-size distribution. An alternative proposed here uses the ratio A38/A25 to estimate the mass-weighted average raindrop size Dm. Rainwater content is then estimated from measurements of polarization propagation differential phase shift (Phi-DP) divided by (1-R), where R is the mass-weighted mean axis ratio of the raindrops computed from Dm. This paper investigates these two water-content estimators using results from a numerical simulation of observations along a microwave link. From these calculations, it appears that the combination (R, Phi-DP) produces more accurate estimates of W than does A38-25. In addition, by combining microwave estimates of W and the rate of rainfall in still air with the mass-weighted mean terminal fall speed derived using A38/A25, it is possible to detect the potential influence of vertical air motion on the raingage-microwave rainfall comparisons.

Jameson, Arthur R.↗

Observations of Subvisible Cirrus Clouds and Gravity Waves at the Tropical Tropopause

Thin, subvisible cirrus (SVC) clouds at the tropical tropopause have been observed by a number of methods in a variety of observational programs, including in situ sampling and aircraft and space-based lidar. Modeling studies suggest that these clouds play an important role in dehydrating tropospheric air as it enters the stratosphere. This is because particles large enough to have significant fall speeds can form under the conditions of slow cooling that are implied by the large horizontal extent of the SVC sheets. The IR radiation that these clouds absorb, and the upward vertical motion this implies, also makes them candidates for a tropical troposphere-to-stratosphere mass transfer mechanism. They may also play a role in the earth's radiation budget. These sheets were observed on five flights during the Tropical Ozone Transport Experiment (TOTE) by the NASA Langley DIAL lidar aboard NASA's DC-8 research aircraft operating during December 1995 and February 1996 south of Hawaii. The relationship of the SVC's observed during TOTE to convection was not a simple one. One class of SVC's are within 1000 km of the persistent strong convection near 15S (the SPCZ). Trajectory analyses indicated that the SVC air masses have in fact passed through the SPCZ within a few days of observation. These clouds are very close to the tropopause, with maximum potential temperatures not much higher than 370K, consistent with in situ water and total water measurements near the tropopause made during the Stratosphere Troposphere Exchange Project in January 1987 at Darwin, Australia. A second class of SVC's are not immediately downstream of convection. These clouds tend to be higher, reaching potential temperatures of 390K or more. Trajectory analyses indicate that the air in these SVC's originates either in the equatorial western Pacific or along the subtropical jet. In any case, the warm temperatures the SVC air masses encounter just prior to the observation time along the back trajectory imply that the clouds cannot be residual particles from cirrus blowoff, but must form locally as the air move upward and equatorward south of Hawaii. Since all the parcels have encountered colder temperatures than those at the time of observation early in their history, subsynoptic scale temperatures colder than the analysis temperatures appear to be required to explain the formation of ice particles. In fact, the sloping shapes of the SVC's do suggest that they are gravity or inertia-gravity waves. In situ meteorological measurements made by the ER-2 within a day of the DC-8 remote lidar observations show a gravity wave structure near the equator with an estimated period of about 30 hours. This is sufficiently long to allow large particles to form and fall out (thus allowing dehydration). Other ER-2 flights south of Hawaii at other times of year show gravity and inertia-gravity waves with a poleward wavenumber component and significant (5 degrees peak to peak) temperature perturbation.

Pfister, Leonhard↗