Search

An ice crystal single-scattering property database is developed in the microwave spectral region (1 to 874 GHz) to provide the scattering, absorption, and polarization properties of 12 ice crystal habits (10-plate aggregate, 5-plate aggregate, 8-column aggregate, solid hexagonal column, hollow hexagonal column, hexagonal plate, solid bullet rosette, hollow bullet rosette, droxtal, oblate spheroid, prolate spheroid, and sphere) with particle maximum dimensions from 2 µm to 10 mm. For each habit, four temperatures (160, 200, 230, and 270 K) are selected to account for temperature dependence of the ice refractive index. The microphysical and scattering properties include projected area, volume, extinction efficiency, single-scattering albedo, asymmetry factor, and six independent nonzero phase matrix elements (i.e. P 11 , P 12 , P 22 , P 33 , P 43 and P 44 ). The scattering properties are computed by the Invariant Imbedding T-Matrix (II-TM) method and the Improved Geometric Optics Method (IGOM). The computation results show that the temperature dependence of the ice single-scattering properties in the microwave region is significant, particularly at high frequencies. Potential active and passive remote sensing applications of the database are illustrated through radar reflectivity and radiative transfer calculations. For cloud radar applications, ignoring temperature dependence has little effect on ice water content measurements. For passive microwave remote sensing, ignoring temperature dependence may lead to brightness temperature biases up to 5 K in the case of a large ice water path. [ABSTRACT FROM AUTHOR]

To derive the bulk radiative properties of ice clouds, aircraft contrails and snow grains, which are fundamental to atmospheric radiative transfer calculations in downstream applications, it is necessary to accurately simulate the scattering of light by individual ice particles. An ice particle optical property database reported in 2013 (hereafter, TAMUice2013) is updated (hereafter, TAMUice2016) to incorporate recent advances in computation of the optical properties of nonspherical particles. Specifically, we employ the invariant imbedding T-matrix (II-TM) method to compute the optical properties of particles with small to moderate size parameters. Both versions use the Improved Geometric Optics Method (IGOM) to compute the optical properties of large ice crystals, but TAMUice2016 improves the treatment of inhomogeneous waves inside the scattering particles in the case where ice is absorptive such as at infrared wavelengths. To bridge the gap between the extinction efficiencies computed from the II-TM and the IGOM, TAMUice2016 includes spectrally dependent higher order terms of the edge effect in addition to the first order counterpart considered in TAMUice2013. Furthermore, the differences between TAMUice2013 and TAMUice2016 are quantified with respect to the computation of the bulk optical properties of ice clouds. [ABSTRACT FROM AUTHOR]

Constructing an appropriate particle morphology model is essential for realistic simulation of optical properties of atmospheric particles. This paper presents a model for generating surface roughness based on a combination of methods from discrete differential geometry combined with a stochastic partial differential equation for surface evolution introduced by Edwards and Wilkinson. Scattering of light by roughened particles is simulated using the Invariant Imbedding T-Matrix (II-TM) method. The effects of surface roughness on the single-scattering properties, namely, the phase matrix, asymmetry factor, and extinction efficiency, are investigated for a single wavelength in the visible range and for a range of size parameters up to x =50. Three different smooth shapes are considered: spherical, spheroidal, and hexagonal, the latter two in just the “compact particle” case of unit aspect ratio. It is shown that roughness has negligible effects on the optical scattering properties for size parameters less than 20. For size parameters ranging from 20 to 50, the phase matrix elements are more sensitive to the surface roughness than are two important integral optical properties, the extinction efficiency and asymmetry factor. As has been seen in studies using other forms of roughening, the phase function is progressively smoothed as roughness increases. The effect on extinction efficiency is to increase it, and on asymmetry factor is to decrease it. Each of these effects is relatively modest in the size range considered, but the trend of results suggests that greater effects will be seen for size parameters larger than ones considered here. [ABSTRACT FROM AUTHOR]

The accuracy of the physical-geometric optics (PG-O) approximation is examined for the simulation of electromagnetic scattering by nonspherical dielectric particles. This study seeks a better understanding of the tunneling effect on the phase matrix by employing the invariant imbedding method to rigorously compute the zeroth-order Debye series, from which the tunneling efficiency and the phase matrix corresponding to the diffraction and external reflection are obtained. The tunneling efficiency is shown to be a factor quantifying the relative importance of the tunneling effect over the Fraunhofer diffraction near the forward scattering direction. Due to the tunneling effect, different geometries with the same projected cross section might have different diffraction patterns, which are traditionally assumed to be identical according to the Babinet principle. For particles with a fixed orientation, the PG-O approximation yields the external reflection pattern with reasonable accuracy, but ordinarily fails to predict the locations of peaks and minima in the diffraction pattern. The larger the tunneling efficiency, the worse the PG-O accuracy is at scattering angles less than 90°. If the particles are assumed to be randomly oriented, the PG-O approximation yields the phase matrix close to the rigorous counterpart, primarily due to error cancellations in the orientation-average process. Furthermore, the PG-O approximation based on an electric field volume-integral equation is shown to usually be much more accurate than the Kirchhoff surface integral equation at side-scattering angles, particularly when the modulus of the complex refractive index is close to unity. Finally, tunneling efficiencies are tabulated for representative faceted particles. [ABSTRACT FROM AUTHOR]

• A new internally-mixed aerosol optics scheme was developed. • Artificial neural networks were used to parameterize aerosol optical properties. • The inhomogeneous aerosol model had lower absorption than the well-mixed model. • The optics scheme was evaluated using the Weather Research and Forecasting model. • The impact of particle inhomogeneity was obvious in the planetary boundary layer. A new internally-mixed aerosol optics scheme was developed. Aerosols (organic carbon, sea salt, aerosol water, sulfate, and nitrate) were assumed to be internally mixed with black carbon and soil dust (BC and SD, respectively), both of which are treated as cores. In the SD-containing particles, the particle non-sphericity was incorporated. Aerosol optical properties were computed from the invariant imbedding T-matrix method and the Lorenz–Mie theory and then parameterized using a novel method involving artificial neural networks (ANN). The new scheme was evaluated using the Weather Research and Forecasting (WRF) model. A case study involving East Asia was performed on Dec 3-4, 2010. Compared to uniformly mixed aerosols, inhomogeneous mixed aerosols had lower absorbability and stronger solar backscattering, causing a positive aerosol direct radiative effect (DRE) bias at the surface (3.0 and 2.6 W/m2 for the Sichuan Basin and Northern China, respectively) and negative DRE bias at the top of atmosphere ([TOA], −3.5 and −2.2 W/m2, respectively). Inhomogeneity effects caused a large reduction (up to 0.2K) in the atmospheric warming above the planetary boundary layer (PBL) and caused a northerly anomaly at 850 hPa on the order of 0.3 m/s. The non-sphericity effect was only significant over the Tarim Basin. DREs related to non-sphericity decreased at both the surface and TOA (−1.0 and −0.6 W/m2, respectively) due to stronger extinction. Non-sphericity caused more surface dimming and more solar heating, which had opposite effects on temperature perturbation, producing a more obvious impact within the PBL than in the free atmosphere. [ABSTRACT FROM AUTHOR]

Understanding the inherent optical properties (IOPs) of coccoliths and coccolithophores is important in oceanic radiative transfer simulations and remote sensing implementations. In this study, the invariant imbedding T-matrix method (II-TM) is employed to investigate the IOPs of coccoliths and coccolithophores. The Emiliania huxleyi ( Ehux ) coccolith and coccolithophore models are built based on observed biometric parameters including the eccentricity, the number of slits, and the rim width of detached coccoliths. The calcification state that specifies the amount of calcium of a single coccolith is critical in the determination of the size–volume/mass relationship (note, the volume/mass of coccoltihs at different calcification states are different although the diameters are the same). The present results show that the calcification state, namely, under-calcification, normal-calcification, or over-calcification, significantly influences the backscattering cross section and the phase matrix. Furthermore, the linear depolarization ratio of the light scattered by coccoliths is sensitive to the degree of calcification, and provides a potentially valuable parameter for interpreting oceanic remote sensing data. The phase function of an ensemble of randomly oriented coccolithophores has a similar pattern to that of individual coccoliths, but the forward scattering is dominant in the coccolithophores due to the large geometric cross sections. The linear depolarization ratio associated with coccolithophores is found to be larger than that for coccoliths as polarization is more sensitive to multiple scattering than the phase function. The simulated coccolithophore phase matrix numerical results are compared with laboratory measurements. For scattering angles larger than 100°, an increase of the phase function with respect to the scattering angle is confirmed based on the present coccolithophore model while the spherical approximation fails. [ABSTRACT FROM AUTHOR]

Abstract: 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 II-TM 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μm) 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μm, 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−2. [Copyright &y& Elsevier]

• The optical properties of dust particles associated with three sets of refractive indices were compared in a spectrum region (650–1135 cm−1). • The refractive indices have a clear impact on the spectrum slope of brightness temperature. • The refractive index impacts on bright temperature are essentially associated with the height of dust layer and the aerosol optical thickness. In this paper, we modeled dust impacts on atmospheric radiative transfer in the thermal infrared region of 650–1135 cm−1 by using the radiative transfer for TOVs model (RTTOV). We assessed three sets of Asian dust refractive indices (the Volz refractive index implemented in the RTTOV as default, and the Gobi and Taklamakan refractive indices) for brightness temperature simulations. The US 76 standard profile was employed to represent the atmospheric status. Three scenarios of high (532.6~575.2 hPa), middle (661~703.6 hPa) and low (786.3~825.8 hPa) dust layers and four cases of aerosol optical thickness (AOT) (0.1, 0.5, 1.0, and 1.5 at 1000 cm−1) were considered in the simulations. It was found that the refractive indices have a clear impact on simulating the values and the spectral slope of brightness temperature in the mid-thermal infrared band. Brightness temperature difference (BTD) for different refractive indices (Gobi/Taklamakan–Volz) increases with the height of dust layer and the AOT. The BTD reaches ~4 K for the high dust layer (532.6~575.2 hPa) and the AOT of 1.5. Most BTDs (0.5–4 K) are much larger than the noise-equivalent temperature difference of high-spectral resolution infrared sounders. Even for the AOT of 0.1, BTDs for different refractive indices could be close to or larger than 0.2 K. This study emphasizes the requirement of accurate refractive indices according to dust origins for improving brightness temperature simulations and provides valuable references for potential improvement of dust retrievals. [ABSTRACT FROM AUTHOR]

Inclusions of air bubbles or soot particles have significant effects on the single-scattering properties of ice crystals, effects that in turn have significant impacts on the radiation budget of an atmosphere containing the crystals. This study investigates some of the single-scattering effects in the case of hexagonal ice crystals, including effects on the backscattering depolarization ratio, a quantity of practical importance in the interpretation of lidar observations. One distinguishing feature of the study is an investigation of scattering properties at a visible wavelength for a crystal with size parameter ( x ) above 100, a size regime where one expects some agreement between exact methods and geometrical optics methods. This expectation is generally borne out in a test comparison of how the sensitivity of scattering properties to the distribution of a given volume fraction of included air is represented using (i) an approximate Monte Carlo Ray Tracing (MCRT) method and (ii) a numerically exact pseudo-spectral time-domain (PSTD) method. Another distinguishing feature of the study is a close examination, using the numerically exact Invariant-Imbedding T-Matrix (II-TM) method, of how some optical properties of importance to satellite remote sensing vary as the volume fraction of inclusions and size of crystal are varied. Although such an investigation of properties in the x >100 regime faces serious computational burdens that force a large number of idealizations and simplifications in the study, the results nevertheless provide an intriguing glimpse of what is evidently a quite complex sensitivity of optical scattering properties to inclusions of air or soot as volume fraction and size parameter are varied. [ABSTRACT FROM AUTHOR]

• Impacts of cloud scattering properties on FY-3D HIRAS simulations are assessed in ARMS and CRTM. • ARMS and CRTM simulated brightness temperatures are compared with HIRAS observations to provide quantitative evaluation. • Spherical ice cloud approximation has been found to be not appropriate for ARMS model simulations. Cloudy sky spectral radiance at the top of the atmosphere has always been an important while difficult variable to simulate for fast radiative transfer models. In this paper, we focus on examining the impacts of cloud scattering properties on the spectral radiance signature of the High-spectral Resolution Infrared Atmospheric Sounder (HIRAS) onboard the Fengyun-3D (FY-3D) satellite by using the Advanced Radiative transfer Modeling System (ARMS) and the Community Radiative Transfer Model (CRTM). Cloud scattering properties used in the radiative transfer models are critical for modeling the spectral radiance under cloudy sky, which involves choices of appropriate cloud particle models and particle size distributions, etc. Multiple FY-3D HIRAS observations over Southern China and Southeast Asia with ice or liquid water cloud cover on 6 May 2018 are examined, respectively. Vertical atmospheric profiles are derived from the Modern-Era Retrospective analysis for Research and Applications, Version 2 reanalysis product. Cloud property retrievals from the Moderate Resolution Imaging Spectroradiometer are used. Cloud scattering property parameterization schemes based on spherical and nonspherical cloud particle shapes are implemented for liquid water and ice clouds in ARMS and CRTM, respectively. Results show that both ARMS and CRTM can well simulate the radiance at the HIRAS spectral ranges under liquid water cloud condition as compared with the HIRAS observation with mean absolute error (MAE) of brightness temperature of less than 1 K. However, for ice cloud conditions, ARMS model using assumed spherical ice properties exhibits large biases between simulation and observation. CRTM with nonspherical ice properties using 16-stream approximation shows MAE less than 1 K and MAE of about 1 K using 2-stream approximation. [ABSTRACT FROM AUTHOR]

Abstract: The many-body iterative T-matrix (MBIT) method is based on the T-matrix method and focuses on scatterers with large aspect ratios. In this study, the MBIT method is extended so that it can be applied to an infinite scatterer with a 1-D periodic finite structure. A semi-analytical solution of infinite scatterers is obtained. In addition, the MBIT method is theoretically extended to include scatterers without axial rotational symmetry. The direct validation of the present method is performed by comparing the analytic solution for an infinite cylinder and the numerically accurate discrete dipole approximation method whereas indirect validation is conducted by using asymptotic solutions. Furthermore, the study also demonstrates that scattering properties of infinite scatterers can be used as a substitute for finite scatterers with large aspect ratios. [Copyright &y& Elsevier]

• An advanced radiative transfer modeling system has been developed for satellite data assimilation and remote sensing applications. • It focuses on the radiative transfer components for applications to Fengyun satellites. • The accuracy is validated against other model results and the measurements. An advanced radiative transfer modeling system (ARMS) has been developed for satellite data assimilation and remote sensing applications. It inherits the capability of other fast radiative transfer models developed mainly for US and European satellite programs but focuses on the radiative transfer components that are specific for the assimilation of the Fengyun satellites and those sensors not included in existing models. The ARMS forward operator includes a fast transmittance module, a new particle absorption and scattering look-up table, surface emissivity and polarimetric radiative transfer solver. The particle optical properties of aerosols and ice clouds are calculated based on a super-spheroidal T-matrix model. The polarimetric two-stream approximation is adopted as the preferred solver to simulate the Stokes vector and their Jacobians. The simulations of ARMS are compared with those of other models and observations of Fengyun satellites, respectively. The comparison shows close agreement between ARMS and other existing fast radiative transfer models. Therefore, ARMS has robust simulation capabilities for existing satellite sensors. [ABSTRACT FROM AUTHOR]