Engineering topics

Cahalan, R. F.

Publications and source records attributed to Cahalan, R. F..

At least 19 records

NASA NTRS2008

Effects of Doubled CO2 on Tropical Sea-Surface Temperature (SSTs) for Onset of Deep Convection and Maximum SST-GCM Simulations Based Inferences

A primary concern of CO2-induced warming is the associated rise of tropical (10S-10N) seasurface temperatures (SSTs). GISS Model-E was used to produce two sets of simulations-one with the present-day and one with doubled CO2 in the atmosphere. The intrinsic usefulness of model guidance in the tropics was confirmed when the model simulated realistic convective coupling between SSTs and atmospheric soundings and that the simulated-data correlations between SSTs and 300 hPa moiststatic energies were found to be similar to the observed. Model predicted SST limits: (i) one for the onset of deep convection and (ii) one for maximum SST, increased in the doubled C02 case. Changes in cloud heights, cloud frequencies, and cloud mass-fractions showed that convective-cloud changes increased the SSTs, while warmer mixed-layer of the doubled CO2 contained approximately 10% more water vapor; clearly that would be conducive to more intense storms and hurricanes.

Sud, Y. C.↗

NASA NTRS2003

Performance of Goddard Earth Observing System GCM Column Radiation Models under Heterogeneous Cloud Conditions

We test the performance of the shortwave (SW) and longwave (LW) Column Radiation Models (CORAMs) of Chou and collaborators with heterogeneous cloud fields from a global single-day dataset produced by NCAR's Community Atmospheric Model with a 2-D CRM installed in each gridbox. The original SW version of the CORAM performs quite well compared to reference Independent Column Approximation (ICA) calculations for boundary fluxes, largely due to the success of a combined overlap and cloud scaling parameterization scheme. The absolute magnitude of errors relative to ICA are even smaller for the LW CORAM which applies similar overlap. The vertical distribution of heating and cooling within the atmosphere is also simulated quite well with daily-averaged zonal errors always below 0.3 K/d for SW heating rates and 0.6 K/d for LW cooling rates. The SW CORAM's performance improves by introducing a scheme that accounts for cloud inhomogeneity. These results suggest that previous studies demonstrating the inaccuracy of plane-parallel models may have unfairly focused on worst scenario cases, and that current radiative transfer algorithms of General Circulation Models (GCMs) may be more capable than previously thought in estimating realistic spatial and temporal averages of radiative fluxes, as long as they are provided with correct mean cloud profiles. However, even if the errors of the particular CORAMs are small, they seem to be systematic, and the impact of the biases can be fully assessed only with GCM climate simulations.

Oreopoulos, L.↗

NASA NTRS2003

Future Satellite Observations of Solar Irradiance

Required solar irradiance measurements for climate studies include those now being made by the Total Irradiance Monitor (TIM) and the Spectral Irradiance Monitor (SIM) onboard the SORCE satellite, part of the Earth Observing System fleet of NASA satellites. Equivalent or better measures of Total Solar Irradiance (TSI) and Spectral Solar Irradiance (SSI, 200 to 2000 nm) are planned for the post-2010 satellites of the National Polar-orbiting Operational Environmental Satellite System ("OESS). The design life of SORCE is 5 years, so a "Solar Irradiance Gap Filler" EOS mission is being planned for launch in the 2007 time frame, to include the same TSI and SSI measurements. Besides avoiding any gap, overlap of the data sources is also necessary for determination of possible multi-decadal trends in solar irradiance. We discuss these requirements and the impacts of data gaps, and data overlaps, that may occur in the monitoring of the critical solar radiative forcing.

Cahalan, R. F.↗

NASA NTRS2002

The International Intercomparison of 3-Dimensional Radiation Codes

I3RC (International Intercomparison of 3-dimensional Radiation Codes) has as its primary goal to compare a wide variety of three-dimensional (3D) radiative transfer methods applied to Earth's atmosphere, with a few selected cloud fields as input, and a few selected radiative quantities as output. Phases 1 and 2 are now complete, and participants represented institutions in Canada, France, Germany, Russia, the United Kingdom, and the USA, who met for two workshops in Tucson, Arizona USA, and compared results from 5 cloud fields of varying complexity, beginning with simplified atmosphere and surface, and proceeding to more realistic cases. Phase 3 is now underway, focusing on improvement and sharing of 3D radiation code, aided by working groups on "Approximations" and "Open Source". The "Approximations" group has so far focused on diffusive approximate methods in an attempt to gain advantages in execution time, and also to advance the understanding of 3D radiation processes. The "Open Source" subgroup is developing a Monte Carlo radiative transfer toolkit that makes state-of-the-art techniques available to a wide range of users. Activities of both subgroups are further explained at the I3RC website http://i3rc.gsfc.nasa.gov. Participants in 13RC are forming a 3D Working Group under the auspices of the International Radiation Commission, and will meet for this and related activities at a workshop in Tucson in November 2002.

Cahalan, R. F.↗

NASA NTRS2002

Consistency of ARESE II Cloud Absorption Estimates and Sampling Issues

Data from three cloudy days (March 3, 21, 29, 2000) of the ARM Enhanced Shortwave Experiment II (ARESE II) were analyzed. Grand averages of broadband absorptance among three sets of instruments were compared. Fractional solar absorptances were approx. 0.21-0.22 with the exception of March 3 when two sets of instruments gave values smaller by approx. 0.03-0.04. The robustness of these values was investigated by looking into possible sampling problems with the aid of 500 nm spectral fluxes. Grand averages of 500 nm apparent absorptance cover a wide range of values for these three days, namely from a large positive (approx. 0.011) average for March 3, to a small negative (approximately -0.03) for March 21, to near zero (approx. 0.01) for March 29. We present evidence suggesting that a large part of the discrepancies among the three days is due to the different nature of clouds and their non-uniform sampling. Hence, corrections to the grand average broadband absorptance values may be necessary. However, application of the known correction techniques may be precarious due to the sparsity of collocated flux measurements above and below the clouds. Our analysis leads to the conclusion that only March 29 fulfills all requirements for reliable estimates of cloud absorption, that is, the presence of thick, overcast, homogeneous clouds.

Oreopoulos, L.↗

NASA NTRS2002

Cloud Thickness from Offbeam Returns (THOR) Validation Campaign on NASA's P3B Over the ARM/SGP

Physical thickness of a cloud layer, sometimes multiple cloud layers, is a crucial controller of solar heating of the Earth- atmosphere system, which drives the convective processes that produce storm systems. Yet clouds of average optical thickness are opaque to conventional lidar, so their thickness is well estimated only by combining a lidar above and another below cloud, or a radar and lidar on the same side, dual facilities not widely available. Here we report initial observations of a new airborne multiple field of view lidar, capable of determining physical thickness of cloud layers from time signatures of off-beam returns from a I kHz micropulse lidar at 540 rim. For a single layer, the time delay of light returning from the outer diffuse halo of light surrounding the beam entry point, relative to the time delay at beam center, determines the cloud physical thickness. The delay combined with the pulse stretch gives the optical thickness. This halo method requires cloud optical thickness exceeding 2, and improves with cloud thickness, thus complimenting conventional lidar, which cannot penetrate thick clouds. Results are presented from March 25, 2002, when THOR flew a butterfly pattern over the ARM site at 8.3 km, above a thin ice cloud at 5 km, and a thick boundary-layer stratus deck with top at 1.3 km, as shown by THOR channel 1, and a base at about 0.3 km as shown by the ground-based MPL. Additional information is included in the original extended abstract.

Cahalan, R. F.↗

NASA NTRS2000

Influence of the Surface and Cloud Nonuniformities in the Solar Energy Fluxes in the Arctic

Solar energy fluxes reaching the surface and absorbed by it are basic components of the energy balance of the Arctic. They depend mainly on the solar zenith angle, a state of the atmosphere, especially the cloudiness, and the surface albedo. However, they can also be modified by variabilities in the surface albedo and cloud optical thickness. The surface of the Arctic can be highly nonuniform. The surface of the Arctic Ocean, which covers the huge part of the Arctic can be view as a mosaic of sea water, sea ice, snow and, in the melting period, melting ponds. In our paper, results are presented of Monte Carlo simulations of the expected influence of nonuniform cloud structure and nonuniform surface albedo on radiative fluxes at the Arctic surface. In particular, the plane parallel biases in the surface absorptance and atmospheric transmittance are studied. The bias is defined as the difference between the real absorptance or transmittance (i.e. nonuniform conditions) averaged over a given area, and the uniform or plane parallel case with the same mean cloud optical thickness and the same mean surface albedo. The dependence of the biases is analysed with respect to the following: domain averaged values of the cloud optical thickness and surface albedo, scales of their spatial variabilities, correlation between cloud optical thickness and cloud albedo variabilities, cloud height, and the solar zenith angle. Ranges of means and standard deviations of the input parameters typical of Arctic conditions are obtained from the SHEBA experiment.

Rozwadowska, A.↗

NASA NTRS1999

Clouds in ECMWF's 30 KM Resolution Global Atmospheric Forecast Model (TL639)

Global models of the general circulation of the atmosphere resolve a wide range of length scales, and in particular cloud structures extend from planetary scales to the smallest scales resolvable, now down to 30 km in state-of-the-art models. Even the highest resolution models do not resolve small-scale cloud phenomena seen, for example, in Landsat and other high-resolution satellite images of clouds. Unresolved small-scale disturbances often grow into larger ones through non-linear processes that transfer energy upscale. Understanding upscale cascades is of crucial importance in predicting current weather, and in parameterizing cloud-radiative processes that control long term climate. Several movie animations provide examples of the temporal and spatial variation of cloud fields produced in 4-day runs of the forecast model at the European Centre for Medium-Range Weather Forecasts (ECMWF) in Reading, England, at particular times and locations of simultaneous measurement field campaigns. model resolution is approximately 30 km horizontally (triangular truncation TL639) with 31 vertical levels from surface to stratosphere. Timestep of the model is about 10 minutes, but animation frames are 3 hours apart, at timesteps when the radiation is computed. The animations were prepared from an archive of several 4-day runs at the highest available model resolution, and archived at ECMWF. Cloud, wind and temperature fields in an approximately 1000 km X 1000 km box were retrieved from the archive, then approximately 60 Mb Vis5d files were prepared with the help of Graeme Kelly of ECMWF, and were compressed into MPEG files each less than 3 Mb. We discuss the interaction of clouds and radiation in the model, and compare the variability of cloud liquid as a function of scale to that seen in cloud observations made in intensive field campaigns. Comparison of high-resolution global runs to cloud-resolving models, and to lower resolution climate models is leading to better understanding of the upscale cascade and suggesting new cloud-radiation parameterizations for climate models.

Cahalan, R. F.↗

NASA NTRS1999

Cloud Power Spectra-Dependence on Solar Zenith Angle and Wavelength, Implications for Cloud Optical Property Retrievals

Scale breaks (spatial scales at which power-law exponent changes occur) observed in Landsat radiances have proven to be useful indicators of radiative interactions, and have aided the development of improved techniques in the remote sensing of clouds. This work extends previous theoretical studies to absorbing wavelengths by using both Landsat Thematic Mapper (TM) observations and Monte Carlo (MC) simulations to infer the systematic dependencies of power spectral shape on cloud characteristics, illumination conditions, and wavelength. We show that MC simulations operating on a simple fractal model of horizontally inhomogeneous clouds produce power spectra that qualitatively resemble observed spectra. We also show that the decrease in the spectra power-law exponent seen at intermediate scales (referred to as "roughening") as the Sun becomes more oblique is more pronounced at absorbing wavelengths. An automated procedure designed to detect the small scale break location is unable to find systematic differences between TM Band 4 and Band 7, despite the fact that MC simulations point to systematic differences in horizontal fluxes. The effect of these qualitative characteristics of the spatial spectra on the retrieval of cloud optical properties is examined by comparing power spectra of nadir radiances with power spectra of optical properties retrieved using either traditional Independent Pixel Approximation approaches or modifications based on normalized radiance indices and the inverse Non-local Independent Pixel Approximation. Assuming that the actual cloud properties follow perfect scaling behavior at all scales, we show the improvement of the proposed retrieval modifications.

Oreopoulos, Lazaros↗

NASA NTRS1999

Diffuse Reflection of Laser Light From Clouds

Laser light reflected from an aqueous suspension of particles or "cloud" with known thickness and particle size distribution defines the "cloud radiative Green's function", G. G is sensitive to cloud thickness, allowing retrieval of that important quantity. We describe a laboratory simulation of G, useful in design of an offbeam Lidar instrument for remote sensing of cloud thickness. Clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The size distribution extends from 0.5 microns to 25 microns, roughly lognormal, similar to real clouds. Density of suspended spheres is adjusted so photon mean-free-path is about 10 cm, 1000 times smaller than in real clouds. The light source is a Nd:YAG laser at 530 nm. Detectors are flux and photon-counting PMTs, with a glass probe for precise positioning. A Labview 5 VI controls position and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider,and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns. Step size is selectable. Far from the beam, the rate of exponential increase in the beam direction scales as expected from diffusion theory, linearly with cloud thickness, and inversely as the square root of the reduced optical thickness, independent of particle size. Nearer the beam the signal increases faster than exponential and depends on particle size. Results verify 3D Monte Carlo simulations that demonstrate detectability of remotely sensed offbeam returns, without filters at night, with narrow bandpass filter in day.

Cahalan, R. F.↗

NASA NTRS1999

Path Radiance Technique for Retrieving Aerosol Optical Thickness over Land

The key issue in retrieving aerosol optical thickness over land from short-wave satellite radiances is to identify and separate the signal due to scattering by a largely transparent aerosol layer from the noise due to reflection by the background surface, where the signal is relatively uniform compared to the highly inhomogeneous surface contribution. Sensitivity studies in aerosol optical thickness retrievals reveal that the apparent reflectance at the top of the atmosphere is very susceptible to the surface reflectance, especially when aerosol optical thickness is small. Uncertainties associated with surface reflectance estimation can greatly amplify the error of the aerosol optical thickness retrieval. To reduce these uncertainties, we have developed a "path radiance" method to retrieve aerosol optical thickness over land by extending the traditional technique that uses the "dark object" approach to extract the aerosol signal. This method uses the signature of the correlation of visible and mid-IR reflectance at the surface, and couples the correlation with the atmospheric effect. We have applied this method to a Landsat TM (Thematic Mapper) image acquired over the Oklahoma Southern Great Plains (SGP) site of DoE's ARM (Department of Energy's Atmospheric Radiation Measurement) program on September 27, 1997, a very clear day (aerosol optical thickness of 0.07 at 0.5 pm) during the first Landsat IOP (Intensive Observation Period). The retrieved mean aerosol optical thickness for TM band 1 at 0.49 pm and band 3 at 0.66 pm agree very well with the ground-based sun-photometer measurements at the ARM site. The ability to retrieve small aerosol optical thickness makes this path radiance technique promising. More importantly, the path radiance is relatively insensitive to surface inhomogeneity. The retrieved mean path radiances in reflectance units have very small standard deviations for both TM blue and red bands. This small variability of path radiance further supports the current aerosol retrieval method.

Tsay, S.-C.↗

NASA NTRS1999

Path Radiance Technique for Retrieving Aerosol Optical Thickness Over Land

The key issue in retrieving aerosol optical thickness over land from shortwave satellite radiances is to identify and separate the signal due to scattering by a largely transparent aerosol layer from the noise due to reflection by the background surface, where the signal is relatively uniform compared to the highly inhomogeneous surface contribution. Sensitivity studies in aerosol optical thickness retrievals reveal that the apparent reflectance at the top of the atmosphere is very susceptible to the surface reflectance, especially when aerosol optical thickness is small. Uncertainties associated with surface reflectance estimation can greatly amplify the error of the aerosol optical thickness retrieval. To reduce these uncertainties, we have developed a "path radiance" method to retrieve aerosol optical thickness over land by extending the traditional technique that uses the "dark object" approach to extract the aerosol signal. This method uses the signature of the correlation of visible and mid-IR reflectance at the surface, and couples the correlation with the atmospheric effect. We have applied this method to a Landsat TM (Thematic Mapper) image acquired over the Oklahoma Southern Great Plains (SGP) site of DoE's ARM (Department of Energy's Atmospheric Radiation Measurement) program on September 27,1997, a very clear day (aerosol optical thickness of 0.07 at 0.5 microns) during the first Landsat Intensive Observation Period (IOP). The retrieved mean aerosol optical thickness for TM band 1 at 0.49 micron and band 3 at 0.66 micron agree very well with the ground-based sun-photometer measurements at the ARM site. The ability to retrieve small aerosol optical thickness makes this path radiance technique promising. More importantly, the path radiance is relatively insensitive to surface inhomogeneity. The retrieved mean path radiances in reflectance units have very small standard deviations for both TM blue and red bands. This small variability of path radiance further supports the current aerosol retrieval method.

Tsay, S.-C.↗

NASA NTRS1998

Cloud Scaling Properties and Cloud Parameterization

Cloud liquid and cloud traction variability is studied as a function of horizontal scale in the ECMWF forecast model during several 10-day runs at the highest available model resolution, recently refined from approximately 60 km (T213) down to approximately 20 km (T639). At higher resolutions, model plane-parallel albedo biases are reduced, so that models may be tuned to have larger, more realistic, cloud liquid water amounts, However, the distribution of cloud liquid assumed -within- each gridbox, for radiative and thermodynamic computations, depends on ad hoc assumptions that are not necessarily consistent with observed scaling properties, or with scaling properties produced by the model at larger scales. To study the larger-scale cloud properties, ten locations on the Earth are chosen to coincide with locations having considerable surface data available for validation, and representing a variety of climatic regimes, scaling exponents are determined from a range or scales down to model resolution, and are re-computed every three hours, separately for low, medium and high clouds, as well as column-integrated cloudiness. Cloud variability fluctuates in time, due to diurnal, synoptic and other' processes, but scaling exponents are found to be relatively stable. various approaches are considered for applying computed cloud scaling to subgrid cloud distributions used for radiation, beyond simple random or maximal overlap now in common use. Considerably more work is needed to compare model cloud scaling with observations. This will be aided by increased availability of high-resolution surface, aircraft and satellite data, and by the increasing resolution of global models,

Cahalan, R. F.↗

NASA NTRS1983

Interannual variability and climatic noise in satellite-observed outgoing longwave radiation

Upwelling-IR observations of the North Pacific by polar orbiters NOAA 3, 4, 5, and 6 and TIROS-N from 1974 to 1981 are analyzed statistically in terms of interannual variability (IAV) in monthly averages and climatic noise due to short-term weather fluctuations. It is found that although the daily variance in the observations is the same in summer and winter months, and although IAV in winter is smaller than that in summer, the climatic noise in winter is so much smaller that a greater fraction of winter anomalies are statistically significant. The smaller winter climatic noise level is shown to be due to shorter autocorrelation times. It is demonstrated that increasing averaging area does not reduce the climatic noise level, suggesting that continuing collection of high-resolution satellite IR data on a global basis is necessary if better models of short-term variability are to be constructed.

Short, D. A.↗

NASA NTRS1983

Climatological studies of cloudiness

The influence of subresolution clouds on estimates of synoptic-scale averages and the effect of synoptic-scale clouds on estimates of planetary scale interannual variability are discussed. It is shown that climatic noise is strongly dependent on the season and geographic region, so that the largest anomaly is not always the most significant. The role of advection in determining the time-scales and degrees of freedom of synoptic scale cloudiness is empirically demonstrated and described in terms of a simple model of stochastic 'cloud patches'.

Cahalan, R. F.↗

NASA NTRS1982

Sampling errors in the estimation of empirical orthogonal functions

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is 'close' to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

North, G. R.↗

NASA NTRS1982

The latitude dependence of the variance of zonally averaged quantities

Geometric characteristics of the spherical earth are shown to be responsible for the increase of variance with latitude of zonally averaged meteorological statistics. An analytic model is constructed to display the effect of a spherical geometry on zonal averages, employing a sphere labeled with radial unit vectors in a real, stochastic field expanded in complex spherical harmonics. The variance of a zonally averaged field is found to be expressible in terms of the spectrum of the vector field of the spherical harmonics. A maximum variance is then located at the poles, and the ratio of the variance to the zonally averaged grid-point variance, weighted by the cosine of the latitude, yields the zonal correlation typical of the latitude. An example is provided for the 500 mb level in the Northern Hemisphere compared to 15 years of data. Variance is determined to increase north of 60 deg latitude.

North, G. R.↗

NASA NTRS1982

Cloud fluctuation statistics

A space-time statistical analysis of total outgoing infrared radiation is used to determine the gross features of day-to-day cloudiness fluctuations over the Pacific Ocean in summer and winter. IR fluctuations arise from the passage of cloudiness systems through a grid box as well as the creation and destruction of cloudiness in the box. Which process dominates depends upon the size of the box relative to the size, speed and persistence time of a typical cloudiness system. In most regions the statistical analysis yields advection speeds characteristic of 700 mb mean flow with spatial dependence resembling the 300 mb mean flow. Spatial scales less than 2000 km predominate, smaller scales having less persistence. Characteristic time scales are on the order of one or two days, even for a grid box spanning the entire North Pacific storm track. This result is remarkable in view of the much longer time scales commonly associated with atmospheric disturbances. Apparently many cloudiness systems are created and destroyed during the lifetime of a single disturbance.

Cahalan, R. F.↗