SEARCH · Search NASA

Results for “cloud optical thickness”

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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Optical scattering and microphysical properties of subvisual cirrus clouds, and climatic implications

The optical and microphysical properties of subvisual cirrus clouds are derived from ground-based polarization lidar, shortwave radiation flux, and solar corona measurements of two approximately 0.75 km deep cirrus located near the tropopause. The first cloud produced no visual manifestations under excellent viewing conditions, and the second appeared to be a persistent aircraft contrail that was generally visible except in the zenith direction. Average lidar linear depolarization ratios and volume backscatter coefficients for the two clouds were 0.19 and 0.35, and 0.6 x 10 to the -3 and 1.4 x 10 to the -3 /km sr, respectively. It is estimated that the zenith-subvisual cirrus contained ice crystals of 25-micron effective diameter at a mean concentration of 25/1 and ice mass content of 0.2 mg/cu m. The threshold cloud optical thickness for visual-versus-invisible cirrus, derived from both broadband shortwave flux and 0.694 micrometer lidar data, is found to be tau sub c approx equal 0.03. Such tau values are comparable to those of 5 to 10 km deep stratospheric aerosol clouds of volcanic origin and polar stratospheric clouds, which are episodic in nature. Hence, we conclude that if these clouds are a fairly common feature of the upper troposphere, as recent SAGE satellite measurements would suggest, then the impact of natural and contrail subvisual cirrus on the planet's radiation balance may be relatively significant.

Sassen, Kenneth↗

NASA NTRS2000

AROTEL Temperature Retrievals Within PSC's Using Raman Scattering

NASA Goddard Space Flight Center's Airborne Raman Ozone Temperature Lidar (AROTEL) made temperature retrievals within Polar Stratospheric Clouds (PSCs) on several flights during the SAGE III Ozone Loss and Validation Experiment (SOLVE) campaign. The location of the PSCs was confirmed using simultaneously acquired data from the NASA Langley Aerosol lidar. Retrievals were made on flight dates 991207, 991210 and 000127 from just above the aircraft to 25 kilometers geometric altitude. Raman temperature retrievals are, to first order, insensitive to Mie interference because the Raman signals are red shifted by 2331 cm(exp -1) from the initial laser wavelength. Backscattering from clouds and aerosols is consequently not observed by the detector; however, extinction does impact the measurement and limits retrievals to optically thin clouds. Comparisons between retrievals employing Rayleigh and Raman scattering show the Raman temperatures to be significantly warmer than those employing Rayleigh scattering. Uncertainties are a function of altitude: at 25 km they were approximately 3 K. Temperature profiles could not be retrieved for optically thick clouds.

Burris, John↗

NASA NTRS1989

Analytical optical scattering in clouds

An analytical optical model for scattering of light due to lightning by clouds of different geometry is being developed. The self-consistent approach and the equivalent medium concept of Twersky was used to treat the case corresponding to outside illumination. Thus, the resulting multiple scattering problem is transformed with the knowledge of the bulk parameters, into scattering by a single obstacle in isolation. Based on the size parameter of a typical water droplet as compared to the incident wave length, the problem for the single scatterer equivalent to the distribution of cloud particles can be solved either by Mie or Rayleigh scattering theory. The super computing code of Wiscombe can be used immediately to produce results that can be compared to the Monte Carlo computer simulation for outside incidence. A fairly reasonable inverse approach using the solution of the outside illumination case was proposed to model analytically the situation for point sources located inside the thick optical cloud. Its mathematical details are still being investigated. When finished, it will provide scientists an enhanced capability to study more realistic clouds. For testing purposes, the direct approach to the inside illumination of clouds by lightning is under consideration. Presently, an analytical solution for the cubic cloud will soon be obtained. For cylindrical or spherical clouds, preliminary results are needed for scattering by bounded obstacles above or below a penetrable surface interface.

Phanord, Dieudonne D.↗

NASA NTRS2014

Global Free Tropospheric NO2 Abundances Derived Using a Cloud Slicing Technique Applied to Satellite Observations from the Aura Ozone Monitoring Instrument (OMI)

We derive free-tropospheric NO2 volume mixing ratios (VMRs) and stratospheric column amounts of NO2 by applying a cloud slicing technique to data from the Ozone Monitoring Instrument (OMI) on the Aura satellite. In the cloud-slicing approach, the slope of the above-cloud NO2 column versus the cloud scene pressure is proportional to the NO2 VMR. In this work, we use a sample of nearby OMI pixel data from a single orbit for the linear fit. The OMI data include cloud scene pressures from the rotational-Raman algorithm and above-cloud NO2 vertical column density (VCD) (defined as the NO2 column from the cloud scene pressure to the top-of-the-atmosphere) from a differential optical absorption spectroscopy (DOAS) algorithm. Estimates of stratospheric column NO2 are obtained by extrapolating the linear fits to the tropopause. We compare OMI-derived NO2 VMRs with in situ aircraft profiles measured during the NASA Intercontinental Chemical Transport Experiment Phase B (INTEX-B) campaign in 2006. The agreement is generally within the estimated uncertainties when appropriate data screening is applied. We then derive a global seasonal climatology of free-tropospheric NO2 VMR in cloudy conditions. Enhanced NO2 in the free troposphere commonly appears near polluted urban locations where NO2 produced in the boundary layer may be transported vertically out of the boundary layer and then horizontally away from the source. Signatures of lightning NO2 are also shown throughout low and middle latitude regions in summer months. A profile analysis of our cloud slicing data indicates signatures of uplifted and transported anthropogenic NO2 in the middle troposphere as well as lightning-generated NO2 in the upper troposphere. Comparison of the climatology with simulations from the Global Modeling Initiative (GMI) for cloudy conditions (cloud optical thicknesses > 10) shows similarities in the spatial patterns of continental pollution outflow. However, there are also some differences in the seasonal variation of free-tropospheric NO2 VMRs near highly populated regions and in areas affected by lightning-generated NOx. Stratospheric column NO2 obtained from cloud slicing agrees well with other independently-generated estimates, providing further confidence in the free-tropospheric results.

NO2↗

NASA NTRSDate not supplied

Machine Learning Application for Improving Cloud Detection and Phase Determination Over Sunglint Regions for Geostationary Satellites

Cloud detection and phase determination over sunglint regions has been a challenge, especially for geostationary (GEO) satellites. Sunglint is observed when the sunlight specular reflection is at the same viewing angle of the satellite sensor. This intense reflection in the visible channels (VIS) is often comparable to that from optically thick clouds. It also contaminates the shortwave infrared channels (SWIR). Consequently, VIS and SWIR channels become less useful - or not useful- when they are saturated, hampering the detection of cloudy and clear-sky pixels. Sunglint contamination happens frequently and exists nearly in every daytime GEO full disk satellite images. However, sunglint intensity and region are difficult to model due to variable viewing geometry and ocean surface conditions. Moreover, existing physical models do not meet the accuracy required for operational GEO satellite cloud detection. We developed a machine learning algorithm to improve cloud detection in sunglint conditions for the NASA Langley’s Satellite ClOud and radiation Property retrieval System (SatCORPS). This poster presents our recent progress in the algorithm development, validation and applications. The algorithm is validated using collocated SatCORPS GOES-East and GOES-West cloud products. We demonstrate that the machine learning cloud detection in sunglint regions is superior to the traditional approach by improving temporal consistency between sunglint and non-sunglint conditions.

Machine Learning, Cloud detection, Sunglint, SatCO↗

NASA NTRS2002

The Radiative Impact of Smoke Aerosols on Clouds in Southeast Asia

The impact of smoke aerosols generated from biomass burning activities in Southeast Asia on the total (direct and indirect) reflected solar radiation from clouds was investigated using satellite data. Narrowband measurements from UV to near-infrared wavelengths (from SeaWiFS and TOMS) were combined with broadband radiation measurements (from CERES). Using this information, we quantified how smoke aerosols change the cloud forcing spectrally and as a whole in the Southeast Asia region. In this region our results show that smoke is present over large areas of cloud-covered regions, and that the frequency of such occurrences is high in the boreal spring. Depending on the thickness of the smoke aerosol, the reflected solar radiation from clouds could he reduced by as much as 100 Watt/sq m, on average over the March 2000 data. We also found that the reduction in the reflectance of the clouds at 670 nm is large enough to lead to significant errors in cloud optical thickness retrievals from satellites such as AVHRR and MODIS.

Hsu, N. Christina↗

NASA NTRSDate not supplied

International Satellite Cloud Climatology Project (ISCCP) Stage D2 Monthly Cloud Product - Revised Algorithm in Hierarchical Data Format (ISCCP_D2)

Since 1983 an international group of institutions has collected and analyzed satellite radiance measurements from up to five geostationary and two polar orbiting satellites to infer the global distribution of cloud properties and their diurnal, seasonal and interannual variations. The primary focus of the first phase of the project (1983-1995) was the elucidation of the role of clouds in the radiation budget (top of the atmosphere and surface). In the second phase of the project (1995 onwards) the analysis also concerns improving understanding of clouds in the global hydrological cycle. [Location=TROPOSPHERE] [Temporal_Coverage: Start_Date=1983-07-01; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=280 Km; Longitude_Resolution=280 Km; Temporal_Resolution=Monthly].

CLOUD TOP PRESSURE↗

NASA NTRS1974

Atmosphere of Venus: Implications of Venera 8 sunlight measurements

Venera 8 measurements of solar illumination within the atmosphere of Venus are quantitatively analyzed by using a multilayer model atmosphere. The analysis shows that there are at least three different scattering layers in the atmosphere of Venus and the total cloud optical thickness is about 10 or greater. However, because of the nature of the observations, it is not possible to determine the vertical distribution of absorbed solar energy, which would reveal the drive for the atmospheric dynamics and the strength of the greenhouse effect. Future spacecraft observations should be designed to (1) measure both upward and downward solar fluxes, (2) include measurements of the highest cloud layers, and (3) employ narrow-band and broad-band sensors.-

Lacis, A. A.↗

NASA NTRS1990

Cloud optical parameters as derived from the multispectral cloud radiometer

Here, researchers analyze the data taken in the marine stratocumulus intensive field observation component of the First ISCCP Regional Experiment (FIRE), conducted off the California coast in July, 1987 using the Multispectral Cloud Radiometer (MCR). Researchers installed an efficient retrieval scheme for deriving the cloud optical thickness and droplet mode radius, and have described the behavior of the retrieval error. Generally the scheme can retrieve the optical thickness and mode radius adequately unless they are too small; the use of the forward scattering region is more sound than the use of the backscattering portion. Applying the method to real data, researchers derived the two-dimensional distribution of optical thickness and mode radius for a portion of one of the FIRE marine stratocumulus missions. In this case study, the droplet size showed a more uniform distribution than optical thickness with some correlation between large droplet size and small optical thickness. Although researchers can find microphysical reasons for these tendencies, they suspect that the remotely sensed droplet size may be overestimated. As a future problem, they will compare results with in situ data of the droplet size distribution. Also it will be very important to check several reasons why the droplet radius might be overestimated, e.g., soot contamination, effect of escaping photons from the lateral sides of broken clouds, and so on.

Nakajima, Eruyuki↗

NASA NTRSDate not supplied

International Satellite Cloud Climatology Project (ISCCP) Stage D2 Monthly Cloud Product - Revised Algorithm in Native (NAT) Data Format (ISCCP_D2_NAT)

CLOUD OPTICAL THICKNESS↗

NASA NTRS2000

Design and Performance of McRas in SCMs and GEOS I/II GCMs

The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.

Sud, Yogesh C.↗

NASA NTRS2012

The Operational MODIS Cloud Optical and Microphysical Property Product: Overview of the Collection 6 Algorithm and Preliminary Results

Operational Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals of cloud optical and microphysical properties (part of the archived products MOD06 and MYD06, for MODIS Terra and Aqua, respectively) are currently being reprocessed along with other MODIS Atmosphere Team products. The latest "Collection 6" processing stream, which is expected to begin production by summer 2012, includes updates to the previous cloud retrieval algorithm along with new capabilities. The 1 km retrievals, based on well-known solar reflectance techniques, include cloud optical thickness, effective particle radius, and water path, as well as thermodynamic phase derived from a combination of solar and infrared tests. Being both global and of high spatial resolution requires an algorithm that is computationally efficient and can perform over all surface types. Collection 6 additions and enhancements include: (i) absolute effective particle radius retrievals derived separately from the 1.6 and 3.7 !-lm bands (instead of differences relative to the standard 2.1 !-lm retrieval), (ii) comprehensive look-up tables for cloud reflectance and emissivity (no asymptotic theory) with a wind-speed interpolated Cox-Munk BRDF for ocean surfaces, (iii) retrievals for both liquid water and ice phases for each pixel, and a subsequent determination of the phase based, in part, on effective radius retrieval outcomes for the two phases, (iv) new ice cloud radiative models using roughened particles with a specified habit, (v) updated spatially-complete global spectral surface albedo maps derived from MODIS Collection 5, (vi) enhanced pixel-level uncertainty calculations incorporating additional radiative error sources including the MODIS L1 B uncertainty index for assessing band and scene-dependent radiometric uncertainties, (v) and use of a new 1 km cloud top pressure/temperature algorithm (also part of MOD06) for atmospheric corrections and low cloud non-unity emissivity temperature adjustments.

Platnick, Steven↗

NASA NTRS2014

The Atmospheric Infrared Sounder Version 6 Cloud Products

The version 6 cloud products of the Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Sounding Unit (AMSU) instrument suite are described. The cloud top temperature, pressure, and height and effective cloud fraction are now reported at the AIRS field-of-view (FOV) resolution. Significant improvements in cloud height assignment over version 5 are shown with FOV-scale comparisons to cloud vertical structure observed by the CloudSat 94 GHz radar and the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP). Cloud thermodynamic phase (ice, liquid, and unknown phase), ice cloud effective diameter D(sub e), and ice cloud optical thickness (τ) are derived using an optimal estimation methodology for AIRS FOVs, and global distributions for 2007 are presented. The largest values of tau are found in the storm tracks and near convection in the tropics, while D(sub e) is largest on the equatorial side of the midlatitude storm tracks in both hemispheres, and lowest in tropical thin cirrus and the winter polar atmosphere. Over the Maritime Continent the diurnal variability of tau is significantly larger than for the total cloud fraction, ice cloud frequency, and D(sub e), and is anchored to the island archipelago morphology. Important differences are described between northern and southern hemispheric midlatitude cyclones using storm center composites. The infrared-based cloud retrievals of AIRS provide unique, decadal-scale and global observations of clouds over portions of the diurnal and annual cycles, and capture variability within the mesoscale and synoptic scales at all latitudes.

infrared instruments↗

NASA NTRS1992

Lidar multiple scattering: Dependence on atmospheric parameters

Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).

Wiegner, Matthias↗

NASA NTRSDate not supplied

Development of a Consistent GEOsat Cloud Property Dataset for the CERES Climate Data Record

Cloud properties are critical for understanding the Earth’s radiation budget and cloud feedbacks. At NASA Langley Research Center, the Satellite ClOud and Radiative Property retrieval System (SatCORPS) provides real-time and historical analyses of clouds derived from Geostationary satellite (GEOsat) data for weather and climate applications. For the Clouds and the Earth’s Radiant Energy System (CERES) program, the global constellation of GEOsats has been analyzed since 2000 to help characterize and account for the diurnal cycle of clouds and their radiative impacts in the CERES climate data record. Obtaining consistent cloud properties over the GEOsat data record during the CERES era is a major objective but a significant challenge considering the diversity of imaging capabilities deployed during that time. The GEOsat data analysis approach for the current CERES Edition-4 (Ed4) data products was focused on accuracy and consistency with MODIS by employing as much spectral information as possible from each satellite. However, the inconsistent use of spectral information across GEOsats led to marked discontinuities in the spatial and temporal record of cloud properties that had to be accounted for post facto in downstream CERES processing. This paper reports progress in developing a new GEOsat analysis system for the next CERES edition (Ed5) that has potential to improve cross-platform consistency and continuity. In this approach, the spectral channel complement is limited to just 3-channels during daytime, ~0.65 µm (VIS), ~3.9 µm (NIR), and ~10.8 µm (IR), common to nearly all of the satellites in the record. At night, a 2-channel approach is taken with the NIR and IR, and ~6.7 µm bands that includes a machine learning approach for optically thick cloud properties. A tradeoff is the potential for reduced accuracy particularly using data from the more advanced satellites that have more spectral channels (e.g. SEVIRI, AHI and ABI) that are known to help improve thin cirrus detection, cloud-aerosol discrimination and estimates in other difficult conditions that challenge cloud remote sensing. The new continuity approach is applied to one month of global GEOSat data for each year of the CERES record since 2000 and compared with the Ed4 GEO and MODIS cloud property time series in order to evaluate the level of improved consistency in the GEOsat record and to assess the accuracy impacts. Cloud fraction will also be assessed with CALIPSO data. Outstanding issues and challenges will be discussed. The results are expected to guide future work needed to develop a more robust GEOsat cloud data record for CERES.

CERES CDR↗

NASA NTRS1992

Does UV instrumentation effectively measure ozone abundance?

Measurements of O3 on Mars provide significant information about the chemistry and composition of the atmosphere, including long-term changes. The most extensive and accurate data were inferred from the Mariner 9 UV spectrometer experiment. Mars O3 shows strong seasonal and latitudinal variation, with column abundances ranging from 0.2 microns at equatorial latitudes to 60 microns over the northern winter polar latitudes (1 micron-atm is a column abundance of 2.689 x 10(exp 15) molecules cm(exp-2)). The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms in one of its two spectral channels every 3 seconds with a spectral resolution of 15 Angstroms and an effective field-of-view of approximately 300 km(exp 2). Measurements were made for almost half a Martian year, with winter and spring in the Northern Hemisphere and summer and fall in the Southern Hemisphere. The detectability limit of the spectrometer was approximately 3 microns of ozone. The UV spectrometer on Mariner 9 was incapable of penetrating the dust during dust storms; the single-scattering albedo and phase function of airborne dust and cloud ice are not known to the degree required to extract the small UV signal reflected up from near the surface. The reflectance spectroscopy technique would also have difficulty detecting the total column abundance of O3 in cases where large dust abundances exist together with the polar hood, especially at high latitudes where large solar zenith angles magnify those optical depths; yet these cases would contain the maximum O3, based on theoretical results. It is quite possible that the maximum O3 column abundance observed by Mariner 9 of 60 microns is common. In fact, larger quantities may exist in some of the colder areas with optically thick clouds and dust. As the Viking period often had more atmospheric dust loading than did that of Mariner 9, the reflectance spectroscopy technique may even have been incapable of detecting the entire O3 column abundance during much of the Mars year that Viking observed, particularly at high latitudes.

Lindner, Bernhard Lee↗

NASA NTRS1998

Estimates of Spectral UV Irradiance from Earth Probe TOMS: Comparisons with Ground-Based Measurements

Estimates of the spectral UV flux incident on the Earth's surface are calculated based on total column ozone and cloud optical thickness determined from the Earth Probe Total Ozone Mapping Spectrometer (TOMS) satellite data. Spectral fluxes are calculated between 300 and 325 nanometers using a weighting function similar to the measured function of the Brewer Spectrophotometer. Comparisons of the TOMS-derived values with ground-based Brewer Spectrophotometer measurements at 3 locations from the Canadian Brewer network (Toronto, Saskatoon and Saturna Island) and 3 locations from the United States EPA network (Boston, Gaithersburg, and Boulder) show reasonably good agreement over a variety of clear and cloudy conditions. Some systematic differences are apparent, particularly when snow/ice is present and the TOMS instrument cannot distinguish between clouds and snow covered ground, thus leading to an underestimation of UV flux. The presence of absorbing aerosols near the ground or sub-pixel clouds can also lead to significant errors. There are also possible errors in the Brewer data due to radiometric calibration uncertainties and uncorrected cosine response.

Labow, G. J.↗

NASA NTRS2013

Global Weather States and Their Properties from Passive and Active Satellite Cloud Retrievals

In this study, the authors apply a clustering algorithm to International Satellite Cloud Climatology Project (ISCCP) cloud optical thickness-cloud top pressure histograms in order to derive weather states (WSs) for the global domain. The cloud property distribution within each WS is examined and the geographical variability of each WS is mapped. Once the global WSs are derived, a combination of CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) vertical cloud structure retrievals is used to derive the vertical distribution of the cloud field within each WS. Finally, the dynamic environment and the radiative signature of the WSs are derived and their variability is examined. The cluster analysis produces a comprehensive description of global atmospheric conditions through the derivation of 11 WSs, each representing a distinct cloud structure characterized by the horizontal distribution of cloud optical depth and cloud top pressure. Matching those distinct WSs with cloud vertical profiles derived from CloudSat and CALIPSO retrievals shows that the ISCCP WSs exhibit unique distributions of vertical layering that correspond well to the horizontal structure of cloud properties. Matching the derived WSs with vertical velocity measurements shows a normal progression in dynamic regime when moving from the most convective to the least convective WS. Time trend analysis of the WSs shows a sharp increase of the fair-weather WS in the 1990s and a flattening of that increase in the 2000s. The fact that the fair-weather WS is the one with the lowest cloud radiative cooling capability implies that this behavior has contributed excess radiative warming to the global radiative budget during the 1990s.

histograms↗