Engineering PapersSearch

Engineering topics

Yang, Ping

Publications and source records attributed to Yang, Ping.

At least 37 records · Page 2

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing

Influence of Ice Particle Surface Roughening on the Global Cloud Radiative Effect

Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for severely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and small/large aggregates of plates). The scattering properties for these individual habits incorporate recent advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10 W m(exp −2) over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface roughness is estimated to be roughly 1-2 W m(exp −2). The CAM results indicate that ice particle roughening can result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud radiative effect.

clouds

A Numerical Combination of Extended Boundary Condition Method and Invariant Imbedding Method Applied to Light Scattering by Large Spheroids and Cylinders

The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.

invariant imbeddings

Impact of Radiatively Interactive Dust Aerosols in the NASA GEOS-5 Climate Model: Sensitivity to Dust Particle Shape and Refractive Index

We investigate the radiative effects of dust aerosols in the NASA GEOS-5 atmospheric general circulation model. GEOS-5 is improved with the inclusion of a sectional aerosol and cloud microphysics module, the Community Aerosol and Radiation Model for Atmospheres (CARMA). Into CARMA we introduce treatment of the dust and sea salt aerosol lifecycle, including sources, transport evolution, and sinks. The aerosols are radiatively coupled to GEOS-5, and we perform a series of multi-decade AMIP-style simulations in which dust optical properties (spectral refractive index and particle shape distribution) are varied. Optical properties assuming spherical dust particles are from Mie theory, while those for non-spherical shape distributions are drawn from a recently available database for tri-axial ellipsoids. The climatologies of the various simulations generally compare well to data from the MODIS, MISR, and CALIOP space-based sensors, the ground-based AERONET, and surface measurements of dust deposition and concentration. Focusing on the summertime Saharan dust cycle we show significant variability in our simulations resulting from different choices of dust optical properties. Atmospheric heating due to dust enhances surface winds over important Saharan dust sources, and we find a positive feedback where increased dust absorption leads to increased dust emissions. We further find that increased dust absorption leads to a strengthening of the summertime Hadley cell circulation, increasing dust lofting to higher altitudes and strengthening the African Easterly Jet. This leads to a longer atmospheric residence time, higher altitude, and generally more northward transport of dust in simulations with the most absorbing dust optical properties. We find that particle shape, although important for radiance simulations, is a minor effect compared to choices of refractive index, although total atmospheric forcing is enhanced by greater than 10 percent for simulations incorporating a spheroidal shape distribution versus ellipsoidal or spherical shapes.

radiative forcing

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties

Efficient Implementation of the Invariant Imbedding T-Matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.

Invariant imbedding method

Simulations of Infrared Radiances Over a Deep Convective Cloud System Observed During TC4: Potential for Enhancing Nocturnal Ice Cloud Retrievals

Retrievals of ice cloud properties using infrared measurements at 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 microns can provide consistent results regardless of solar illumination, but are limited to cloud optical thicknesses tau < approx.6. This paper investigates the variations in radiances at these wavelengths over a deep convective cloud system for their potential to extend retrievals of tau and ice particle size D(sub e) to optically thick clouds. Measurements from the Moderate Resolution Imaging Spectroradiometer Airborne Simulator--ASTER, the Scanning High-resolution Interferometer Sounder, the Cloud Physics Lidar (CPL), and the Cloud Radar System (CRS) aboard the NASA ER-2 aircraft during the NASA TC4 (Tropical Composition, Cloud and Climate Coupling) experiment flight during 5 August 2007, are used to examine the retrieval capabilities of infrared radiances over optically thick ice clouds. Simulations based on coincident in-situ measurements and combined cloud tau from CRS and CPL measurements are comparable to the observations. They reveal that brightness temperatures at these bands and their differences (BTD) are sensitive to tau up to approx.20 and that for ice clouds having tau > 20, the 3.7 - 10.8 microns and 3.7 - 6.7 microns BTDs are the most sensitive to D(sub e). Satellite imagery appears consistent with these results. Keywords: clouds; optical depth; particle size; satellite; TC4; multispectral thermal infrared

Minnis, Patrick

Simulations of Infrared Radiances Over a Deep Convective Cloud System Observed During TC4- Potential for Enhancing Nocturnal Ice Cloud Retrievals

Retrievals of ice cloud properties using infrared measurements at 3.7, 6.7, 7.3, 8.5, 10.8, and 12.0 microns can provide consistent results regardless of solar illumination, but are limited to cloud optical thicknesses tau < approx.6. This paper investigates the variations in radiances at these wavelengths over a deep convective cloud system for their potential to extend retrievals of tau and ice particle size D(sub e) to optically thick clouds. Measurements from the Moderate Resolution Imaging Spectroradiometer Airborne Simulator--ASTER, the Scanning High-resolution Interferometer Sounder, the Cloud Physics Lidar (CPL), and the Cloud Radar System (CRS) aboard the NASA ER-2 aircraft during the NASA TC4 (Tropical Composition, Cloud and Climate Coupling) experiment flight during 5 August 2007, are used to examine the retrieval capabilities of infrared radiances over optically thick ice clouds. Simulations based on coincident in-situ measurements and combined cloud tau from CRS and CPL measurements are comparable to the observations. They reveal that brightness temperatures at these bands and their differences (BTD) are sensitive to tau up to approx.20 and that for ice clouds having tau > 20, the 3.7 - 10.8 microns and 3.7 - 6.7 microns BTDs are the most sensitive to D(sub e). Satellite imagery appears consistent with these results. Keywords: clouds; optical depth; particle size; satellite; TC4; multispectral thermal infrared

CRS and CPL

Determination of Ice Cloud Models Using MODIS and MISR Data

Representation of ice clouds in radiative transfer simulations is subject to uncertainties associated with the shapes and sizes of ice crystals within cirrus clouds. In this study, we examined several ice cloud models consisting of smooth, roughened, homogeneous and inhomogeneous hexagonal ice crystals with various aspect ratios. The sensitivity of the bulk scattering properties and solar reflectances of cirrus clouds to specific ice cloud models is investigated using the improved geometric optics method (IGOM) and the discrete ordinates radiative transfer (DISORT) model. The ice crystal habit fractions in the ice cloud model may significantly affect the simulations of cloud reflectances. A new algorithm was developed to help determine an appropriate ice cloud model for application to the satellite-based retrieval of ice cloud properties. The ice cloud particle size retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) data, collocated with Multi-angle Imaging Spectroradiometer (MISR) observations, is used to infer the optical thicknesses of ice clouds for nine MISR viewing angles. The relative differences between view-dependent cloud optical thickness and the averaged value over the nine MISR viewing angles can vary from −0.5 to 0.5 and are used to evaluate the ice cloud models. In the case for 2 July 2009, the ice cloud model with mixed ice crystal habits is the best fit to the observations (the root mean square (RMS) error of cloud optical thickness reaches 0.365). This ice cloud model also produces consistent cloud property retrievals for the nine MISR viewing configurations within the measurement uncertainties.

MODIS

Comparison of PARASOL Observations with Polarized Reflectances Simulated Using Different Ice Habit Mixtures

Insufficient knowledge of the habit distribution and the degree of surface roughness of ice crystals within ice clouds is a source of uncertainty in the forward light scattering and radiative transfer simulations required in downstream applications involving these clouds. The widely used MODerate Resolution Imaging Spectroradiometer (MODIS) Collection 5 ice microphysical model assumes a mixture of various ice crystal shapes with smooth-facets except aggregates of columns for which a moderately rough condition is assumed. When compared with PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) polarized reflection data, simulations of polarized reflectance using smooth particles show a poor fit to the measurements, whereas very rough-faceted particles provide an improved fit to the polarized reflectance. In this study a new microphysical model based on a mixture of 9 different ice crystal habits with severely roughened facets is developed. Simulated polarized reflectance using the new ice habit distribution is calculated using a vector adding-doubling radiative transfer model, and the simulations closely agree with the polarized reflectance observed by PARASOL. The new general habit mixture is also tested using a spherical albedo differences analysis, and surface roughening is found to improve the consistency of multi-angular observations. It is suggested that an ice model incorporating an ensemble of different habits with severely roughened surfaces would potentially be an adequate choice for global ice cloud retrievals.

Cole, Benjamin H.

Ultraspectral Sounding Retrieval Error Budget and Estimation

The ultraspectral infrared radiances obtained from satellite observations provide atmospheric, surface, and/or cloud information. The intent of the measurement of the thermodynamic state is the initialization of weather and climate models. Great effort has been given to retrieving and validating these atmospheric, surface, and/or cloud properties. Error Consistency Analysis Scheme (ECAS), through fast radiative transfer model (RTM) forward and inverse calculations, has been developed to estimate the error budget in terms of absolute and standard deviation of differences in both spectral radiance and retrieved geophysical parameter domains. The retrieval error is assessed through ECAS without assistance of other independent measurements such as radiosonde data. ECAS re-evaluates instrument random noise, and establishes the link between radiometric accuracy and retrieved geophysical parameter accuracy. ECAS can be applied to measurements of any ultraspectral instrument and any retrieval scheme with associated RTM. In this paper, ECAS is described and demonstration is made with the measurements of the METOP-A satellite Infrared Atmospheric Sounding Interferometer (IASI)..

Zhou, Daniel K.

Effect of Thin Cirrus Clouds on Dust Optical Depth Retrievals From MODIS Observations

The effect of thin cirrus clouds in retrieving the dust optical depth from MODIS observations is investigated by using a simplified aerosol retrieval algorithm based on the principles of the Deep Blue aerosol property retrieval method. Specifically, the errors of the retrieved dust optical depth due to thin cirrus contamination are quantified through the comparison of two retrievals by assuming dust-only atmospheres and the counterparts with overlapping mineral dust and thin cirrus clouds. To account for the effect of the polarization state of radiation field on radiance simulation, a vector radiative transfer model is used to generate the lookup tables. In the forward radiative transfer simulations involved in generating the lookup tables, the Rayleigh scattering by atmospheric gaseous molecules and the reflection of the surface assumed to be Lambertian are fully taken into account. Additionally, the spheroid model is utilized to account for the nonsphericity of dust particles In computing their optical properties. For simplicity, the single-scattering albedo, scattering phase matrix, and optical depth are specified a priori for thin cirrus clouds assumed to consist of droxtal ice crystals. The present results indicate that the errors in the retrieved dust optical depths due to the contamination of thin cirrus clouds depend on the scattering angle, underlying surface reflectance, and dust optical depth. Under heavy dusty conditions, the absolute errors are comparable to the predescribed optical depths of thin cirrus clouds.

Feng, Qian

Global Land Surface Emissivity Retrieved From Satellite Ultraspectral IR Measurements

Ultraspectral resolution infrared (IR) radiances obtained from nadir observations provide information about the atmosphere, surface, aerosols, and clouds. Surface spectral emissivity (SSE) and surface skin temperature from current and future operational satellites can and will reveal critical information about the Earth s ecosystem and land-surface-type properties, which might be utilized as a means of long-term monitoring of the Earth s environment and global climate change. In this study, fast radiative transfer models applied to the atmosphere under all weather conditions are used for atmospheric profile and surface or cloud parameter retrieval from ultraspectral and/or hyperspectral spaceborne IR soundings. An inversion scheme, dealing with cloudy as well as cloud-free radiances observed with ultraspectral IR sounders, has been developed to simultaneously retrieve atmospheric thermodynamic and surface or cloud microphysical parameters. This inversion scheme has been applied to the Infrared Atmospheric Sounding Interferometer (IASI). Rapidly produced SSE is initially evaluated through quality control checks on the retrievals of other impacted surface and atmospheric parameters. Initial validation of retrieved emissivity spectra is conducted with Namib and Kalahari desert laboratory measurements. Seasonal products of global land SSE and surface skin temperature retrieved with IASI are presented to demonstrate seasonal variation of SSE.

Zhou, D. K.

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu

How Well Can Infrared Sounders Observe the Atmosphere and Surface Through Clouds?

Infrared sounders, such as the Atmospheric Infrared Sounder (AIRS), the Infrared Atmospheric Sounding Interferometer (IASI), and the Cross-track Infrared sounder (CrIS), have a cloud-impenetrable disadvantage in observing the atmosphere and surface under opaque cloudy conditions. However, recent studies indicate that hyperspectral, infrared sounders have the ability to detect cloud effective-optical and microphysical properties and to penetrate optically thin clouds in observing the atmosphere and surface to a certain degree. We have developed a retrieval scheme dealing with atmospheric conditions with cloud presence. This scheme can be used to analyze the retrieval accuracy of atmospheric and surface parameters under clear and cloudy conditions. In this paper, we present the surface emissivity results derived from IASI global measurements under both clear and cloudy conditions. The accuracy of surface emissivity derived under cloudy conditions is statistically estimated in comparison with those derived under clear sky conditions. The retrieval error caused by the clouds is shown as a function of cloud optical depth, which helps us to understand how well infrared sounders can observe the atmosphere and surface through clouds.

Zhou, Daniel K.

A Principal Component-Based Radiative Transfer Forward Model (PCRTM) for Vertically in Homogeneous Cloud

A principal-component based radiative transfer model (PCRTM) is developed for simulating the infrared spectral radiance at the top of the atmosphere (TOA). The PCRTM approach removes the redundancy in radiative transfer calculation in high resolution infrared spectra, and saves significant amount of computational time with great accuracy. In PCRTM, both ice and water clouds are treated as effective transmissivity and reflectivity stored in a pre-calculated lookup tables. These quantities are calculated using cloud single scattering properties such as cloud optical depth, cloud particle size, cloud phase, etc. The cloud can be inserted into any pressure layer in the PCRTM model (up to 100 layers). The effective temperature of each cloud layer is treated as a function of its optical depth. To test the accuracy of this approximation, the results are compared with the more rigorous DISORT model, which treats cloud as a plane parallel layer. The root-mean-square error of PCRTM, with respect to DISORT results, is generally less than 0.5 K in brightness temperature. However, the CPU time required by PCRTM was approximately two orders of magnitude less than that required by DISORT.

Li, Hui

Hyperspectrally-Resolved Surface Emissivity Derived Under Optically Thin Clouds

Surface spectral emissivity derived from current and future satellites can and will reveal critical information about the Earth s ecosystem and land surface type properties, which can be utilized as a means of long-term monitoring of global environment and climate change. Hyperspectrally-resolved surface emissivities are derived with an algorithm utilizes a combined fast radiative transfer model (RTM) with a molecular RTM and a cloud RTM accounting for both atmospheric absorption and cloud absorption/scattering. Clouds are automatically detected and cloud microphysical parameters are retrieved; and emissivity is retrieved under clear and optically thin cloud conditions. This technique separates surface emissivity from skin temperature by representing the emissivity spectrum with eigenvectors derived from a laboratory measured emissivity database; in other words, using the constraint as a means for the emissivity to vary smoothly across atmospheric absorption lines. Here we present the emissivity derived under optically thin clouds in comparison with that under clear conditions.

Zhou, Daniel K.