scispace - formally typeset
Search or ask a question

Showing papers on "Radiative transfer published in 1997"


Journal ArticleDOI
TL;DR: A rapid and accurate radiative transfer model (RRTM) for climate applications has been developed and the results extensively evaluated as discussed by the authors, which is performed using the correlated-k method: the k distributions are attained directly from the LBLRTM line-byline model, which connects the absorption coefficients used by RRTM to high-resolution radiance validations done with observations.
Abstract: A rapid and accurate radiative transfer model (RRTM) for climate applications has been developed and the results extensively evaluated. The current version of RRTM calculates fluxes and cooling rates for the longwave spectral region (10–3000 cm−1) for an arbitrary clear atmosphere. The molecular species treated in the model are water vapor, carbon dioxide, ozone, methane, nitrous oxide, and the common halocarbons. The radiative transfer in RRTM is performed using the correlated-k method: the k distributions are attained directly from the LBLRTM line-by-line model, which connects the absorption coefficients used by RRTM to high-resolution radiance validations done with observations. Refined methods have been developed for treating bands containing gases with overlapping absorption, for the determination of values of the Planck function appropriate for use in the correlated-k approach, and for the inclusion of minor absorbing species in a band. The flux and cooling rate results of RRTM are linked to measurement through the use of LBLRTM, which has been substantially validated with observations. Validations of RRTM using LBLRTM have been performed for the midlatitude summer, tropical, midlatitude winter, subarctic winter, and four atmospheres from the Spectral Radiance Experiment campaign. On the basis of these validations the longwave accuracy of RRTM for any atmosphere is as follows: 0.6 W m−2 (relative to LBLRTM) for net flux in each band at all altitudes, with a total (10–3000 cm−1) error of less than 1.0 W m−2 at any altitude; 0.07 K d−1 for total cooling rate error in the troposphere and lower stratosphere, and 0.75 K d−1 in the upper stratosphere and above. Other comparisons have been performed on RRTM using LBLRTM to gauge its sensitivity to changes in the abundance of specific species, including the halocarbons and carbon dioxide. The radiative forcing due to doubling the concentration of carbon dioxide is attained with an accuracy of 0.24 W m−2, an error of less than 5%. The speed of execution of RRTM compares favorably with that of other rapid radiation models, indicating that the model is suitable for use in general circulation models.

6,861 citations


Journal ArticleDOI
TL;DR: The 6S code has still limitations; it cannot handle spherical atmosphere and as a result, it cannot be used for limb observations, and the decoupling the authors are using for absorption and scattering effects does not allow to use the code in presence of strong absorption bands.
Abstract: Remote sensing from satellite or airborne platforms of land or sea surfaces in the visible and near infrared is strongly affected by the presence of the atmosphere along the path from Sun to target (surface) to sensor. This paper presents 6S (Second Simulation of the Satellite Signal in the Solar Spectrum), a computer code which can accurately simulate the above problems. The 6S code is an improved version of 5S (Simulation of the Satellite Signal in the Solar Spectrum), developed by the Laboratoire d'Optique Atmospherique ten years ago. The new version now permits calculations of near-nadir (down-looking) aircraft observations, accounting for target elevation, non lambertian surface conditions, and new absorbing species (CH/sub 4/, N/sub 2/O, CO). The computational accuracy for Rayleigh and aerosol scattering effects has been improved by the use of state-of-the-art approximations and implementation of the successive order of scattering (SOS) algorithm. The step size (resolution) used for spectral integration has been improved to 2.5 nm. The goal of this paper is not to provide a complete description of the methods used as that information is detailed in the 6S manual, but rather to illustrate the impact of the improvements between 5S and 6S by examining some typical remote sensing situations. Nevertheless, the 6S code has still limitations. It cannot handle spherical atmosphere and as a result, it cannot be used for limb observations. In addition, the decoupling the authors are using for absorption and scattering effects does not allow to use the code in presence of strong absorption bands.

2,955 citations


Journal ArticleDOI
TL;DR: This paper examined the sensitivity of a climate model to a wide range of radiative forcings, including changes of solar irradiance, atmospheric CO2, O3, CFCs, clouds, aerosols, surface albedo, and a "ghost" forcing introduced at arbitrary heights, latitudes, longitudes, seasons, and times of day.
Abstract: We examine the sensitivity of a climate model to a wide range of radiative forcings, including changes of solar irradiance, atmospheric CO2, O3, CFCs, clouds, aerosols, surface albedo, and a “ghost” forcing introduced at arbitrary heights, latitudes, longitudes, seasons, and times of day. We show that, in general, the climate response, specifically the global mean temperature change, is sensitive to the altitude, latitude, and nature of the forcing; that is, the response to a given forcing can vary by 50% or more depending upon characteristics of the forcing other than its magnitude measured in watts per square meter. The consistency of the response among different forcings is higher, within 20% or better, for most of the globally distributed forcings suspected of influencing global mean temperature in the past century, but exceptions occur for certain changes of ozone or absorbing aerosols, for which the climate response is less well behaved. In all cases the physical basis for the variations of the response can be understood. The principal mechanisms involve alterations of lapse rate and decrease (increase) of large-scale cloud cover in layers that are preferentially heated (cooled). Although the magnitude of these effects must be model-dependent, the existence and sense of the mechanisms appear to be reasonable. Overall, we reaffirm the value of the radiative forcing concept for predicting climate response and for comparative studies of different forcings; indeed, the present results can help improve the accuracy of such analyses and define error estimates. Our results also emphasize the need for measurements having the specificity and precision needed to define poorly known forcings such as absorbing aerosols and ozone change. Available data on aerosol single scatter albedo imply that anthropogenic aerosols cause less cooling than has commonly been assumed. However, negative forcing due to the net ozone change since 1979 appears to have counterbalanced 30–50% of the positive forcing due to the increase of well-mixed greenhouse gases in the same period. As the net ozone change includes halogen-driven ozone depletion with negative radiative forcing and a tropospheric ozone increase with positive radiative forcing, it is possible that the halogen-driven ozone depletion has counterbalanced more than half of the radiative forcing due to well-mixed greenhouse gases since 1979.

2,044 citations


Journal ArticleDOI
TL;DR: In this paper, the authors derive hydrostatic, radiative equilibrium models for passive disks surrounding T Tauri stars, where each disk is encased by an optically thin layer of superheated dust grains.
Abstract: We derive hydrostatic, radiative equilibrium models for passive disks surrounding T Tauri stars. Each disk is encased by an optically thin layer of superheated dust grains. This layer reemits directly to space about half the stellar energy it absorbs. The other half is emitted inward and regulates the interior temperature of the disk. The heated disk flares. As a consequence, it absorbs more stellar radiation, especially at large radii, than a flat disk would. The portion of the spectral energy distribution contributed by the disk is fairly flat throughout the thermal infrared. At fixed frequency, the contribution from the surface layer exceeds that from the interior by about a factor 3 and is emitted at more than an order of magnitude greater radius. Spectral features from dust grains in the superheated layer appear in emission if the disk is viewed nearly face-on.

1,642 citations


Journal ArticleDOI
TL;DR: A simple method for the estimation of global leaf area index (LAI) and fraction of photosynthetically active radiation absorbed by the vegetation (FAPAR) from atmospherically corrected Normalized Difference Vegetation Index (NDVI) observations is described.
Abstract: A simple method for the estimation of global leaf area index (LAI) and fraction of photosynthetically active radiation absorbed by the vegetation (FAPAR) from atmospherically corrected Normalized Difference Vegetation Index (NDVI) observations is described. Recent improvements to the authors' three dimensional radiative transfer model of a vegetated surface are described. Example simulation results and a validation exercise are discussed. The model was utilized to derive land cover specific NDVI-LAI and NDVI-FAPAR relations. The method therefore requires stratification of global vegetation into cover types that are compatible with the radiative transfer model. Such a classification based on vegetation structure is proposed and a simple method for its derivation is presented. Proof-of-concept results are given to illustrate the feasibility of the proposed method.

979 citations


Journal ArticleDOI
TL;DR: In this paper, the spectral radiance of an aerosol with a bilognormal size distribution is simulated with a single lognormal aerosol, with an appropriate mean radius and width of distribution.
Abstract: Spectral radiances measured at the top of the atmosphere in a wide spectral range (0.55–2.13 μm) are used to monitor the aerosol optical thickness and the aerosol size distribution (integrated on the vertical column) of the ambient (undisturbed) aerosol over the oceans. Even for the moderate resolution imaging spectrometer (MODIS) wide spectral range, only three parameters that describe the aerosol loading and size distribution can be retrieved. These three parameters are not always unique. For instance, the spectral radiance of an aerosol with a bilognormal size distribution can be simulated very well with a single lognormal aerosol with an appropriate mean radius and width of distribution. Preassumptions on the general structure of the size distribution are therefore required in the inversion of MODIS data. The retrieval of the aerosol properties is performed using lookup table computations. The volume size distribution in the lookup table is described with two lognormal modes: a single mode to describe the accumulation mode particles (radius 1.0 μm). Note that two accumulation modes may be present, one dominated by gas phase processes and a second dominated by cloud phase processes. The coarse mode can also be split into several partially overlapping modes describing maritime salt particles and dust. The aerosol parameters we expect to retrieve are η, the fractional contribution of the accumulation mode to scattering; τ, the spectral optical thickness; and rm, the mean particle size of the dominant mode. Additional radiative quantities such as asymmetry parameter and effective radius are derived subsequently. The impact of the surface conditions, wind speed and chlorophyll content on the retrieval is estimated, the impact of potential sources of error like the calibration of the instrument is also tested. The algorithm has been applied successfully to actual data sets provided by the Thematic Mapper on Landsat 5 and by the MODIS airborne simulator on the ER-2 and tested against ground and airborne measurements. A first estimate of the general accuracy is Δτ = ±0.05±0.05τ (at 550 nm), Δrm = 0.3rm, Δη = ±0.25.

944 citations


Journal ArticleDOI
TL;DR: The diffusion approximation of the radiative transfer equation is a model used widely to describe photon migration in highly diffusing media and is an important matter in biological tissue optics as mentioned in this paper, however, it is not a suitable model for biological tissue imaging.
Abstract: The diffusion approximation of the radiative transfer equation is a model used widely to describe photon migration in highly diffusing media and is an important matter in biological tissue optics An analysis of the time-dependent diffusion equation together with its solutions for the slab geometry and for a semi-infinite diffusing medium are reported These solutions, presented for both the time-dependent and the continuous wave source, account for the refractive index mismatch between the turbid medium and the surrounding medium The results have been compared with those obtained when different boundary conditions were assumed The comparison has shown that the effect of the refractive index mismatch cannot be disregarded This effect is particularly important for the transmittance The discussion of results also provides an analysis of the role of the absorption coefficient in the expression of the diffusion coefficient

473 citations


Journal ArticleDOI
TL;DR: In this paper, the authors use numerically accurate solutions of the vector radiative transfer equation for a realistic atmosphere-ocean model to theoretically simulate several types of satellite aerosol retrievals over the ocean utilizing radiance measurements alone, polarization measurements alone and radiance and polarization measurements combined.
Abstract: Most current and proposed satellite remote sensing of tropospheric aerosols relies upon radiance measurements that are interpreted using algorithms that determine best fits to precalculated scattered sunlight for one or more "standard" aerosol models. However, the number of different types of aerosol and the substantial space and time variations typically encountered can pose a severe uniqueness problem even for the multiple constraints provided by multispectral radiances of a scene at a number of observation zenith angles. Experience with polarimetremote sensing on planetary missions has demonstrated that the measurement of polarization as well as the radiance can resolve such uniqueness problems. We use numerically accurate solutions of the vector radiative transfer equation for a realistic atmosphere-ocean model to theoretically simulate several types of satellite aerosol retrievals over the ocean utilizing radiance measurements alone, polarization measurements alone, and radiance and polarization measurements combined. We have restricted all simulations to a single near-infrared wavelength of 0.865 tm and assumed that aerosols are 'spherical, monomodal, and nonabsorbing. These simplifications permit a study of practical scope that tests the retrieval algorithms under exactly the same conditions, thus clearly demonstrating their relative capabilities. In agreement with previous analyses, we have found that radiance-only algorithms using multiple-viewing-angle observations perform far better than those based on single-viewing-angle measurements. However, even multiple-viewing-angle radiance measurements taken at a single wavelength are not always sufficient to determine the aerosol optical thickness, effective radius, and refractive index with high enough accuracy. In contrast, high-accuracy, single-wavelength, multiple-viewing-angle polarimetry alone is capable of uniquely retrieving all three aerosol characteristics with extremely high accuracy (+_0.015 in aerosol optical thickness, +_0.03 gm in effective radius, and +_0.01 in refractive index). Furthermore, the accuracy of the optical thickness retrieval can be slightly improved by simultaneously using radiance measurements: Our analysis demonstrates that algorithms utilizing high-accuracy polarization as well as radiance measurements are much less dependent on the availability and use of a priori information and can be expected to provide a physically based retrieval of aerosol characteristics (optical thickness, refractive index, and size) with accuracy needed for long-term monitoring of global climate forcings and feedbacks.

412 citations



Journal ArticleDOI
TL;DR: Improved parameters for the description of Rayleigh scattering in air and for the detailed rotational Raman scattering component for scattering by O (2) and N(2) are presented for the wavelength range 200-1000 nm to enable more accurate calculations of bulk molecular scattering and of the Ring effect for a variety of atmospheric radiative transfer and constituent retrieval applications.
Abstract: Improved parameters for the description of Rayleigh scattering in air and for the detailed rotational Raman scattering component for scattering by O2 and N2 are presented for the wavelength range 200 ‐1000 nm. These parameters enable more accurate calculations to be made of bulk molecular scattering and of the Ring effect for a variety of atmospheric radiative transfer and constituent retrieval applications. A solar reference spectrum with accurate absolute vacuum wavelength calibration, suitable for convolution with the rotational Raman spectrum for Ring effect calculations, has been produced at 0.01-nm resolution from several sources. It is convolved with the rotational Raman spectra of O2 and N2 to produce an atmospheric Ring effect source spectrum. © 1997 Optical Society of America

404 citations


Journal ArticleDOI
TL;DR: In this article, temperature-dependent fluorescence quantum yield and lifetime measurements on the LH-1 and LH-2 complexes of Rhodobacter sphaeroides and on the isolated B820 subunit of Rhodospirillum rubrum were presented.
Abstract: We present temperature-dependent fluorescence quantum yield and lifetime measurements on the LH-1 and LH-2 complexes of Rhodobacter sphaeroides and on the isolated B820 subunit of Rhodospirillum rubrum. From these measurements the superradiance is calculated, which is related to the delocalization of excitations in these complexes. In the B820 preparation we find a radiative rate that is 30% higher than that of monomeric bacteriochlorophyll, in agreement with a dimer model of this subunit. At room temperature both LH-1 and LH-2 are superradiant relative to monomeric Bchl-a with enhancement factors of 3.8 and 2.8, respectively. In LH-2 the radiative rate does not change significantly upon lowering the temperature to 4 K. LH-1 however exhibits a strong temperature dependence, giving rise to a 2.4 times higher radiative rate at 4 K relative to room temperature. From modeling of the superradiance using a Hamiltonian based on the LH-2 structure and including site inhomogeneity, we conclude that the ratio of in...

Journal ArticleDOI
TL;DR: In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models as discussed by the authors.
Abstract: In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W m22 in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (6 10 Wm 2 2). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models’ neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of

Journal ArticleDOI
TL;DR: In this paper, the authors developed a new, fully scale-free method for solving radiative transfer, present exact numerical results, and derive approximate analytical solutions for spherical geometry, covering the entire range of parameter space relevant to observations.
Abstract: Dust infrared emission possesses scaling properties. Overall luminosity is never an input parameter of the radiative transfer problem, spectral shape is the only relevant property of the heating radiation when the inner boundary of the dusty region is controlled by dust sublimation. Similarly, the absolute scales of densities and distances are irrelevant; the geometry enters only through angles, relative thicknesses and aspect ratios, and the actual magnitudes of densities and distances enter only through one independent parameter, the overall optical depth. Dust properties enter only through dimensionless, normalized distributions that describe the spatial variation of density and the wavelength dependence of scattering and absorption efficiencies. Scaling enables a systematic approach to modeling and classification of IR spectra. We develop a new, fully scale-free method for solving radiative transfer, present exact numerical results, and derive approximate analytical solutions for spherical geometry, covering the entire range of parameter space relevant to observations. Scaling implies tight correlations among the SEDs of various members of the same class of sources such as young stellar objects, late-type stars, etc. In particular, all members of the same class occupy common, well defined regions in color-color diagrams. The observational data corroborate the existence of these correlations.

Journal ArticleDOI
TL;DR: In this article, a new particle size bin structure that nearly eliminates numerical diffusion during growth but still treats nucleation, emissions, coagulation, reversible chemistry, condensation, dissolution, evaporation, irreversible chemistry, sedimentation, dry deposition, and radiative scattering and absorption by particles is discussed.

Journal ArticleDOI
TL;DR: In this article, the authors adopt the view that the classes of active galactic nuclei (AGN) known variously as gigahertz peak spectrum (GPS) sources, compact steep-spectrum (CSS) sources and compact symmetric objects (CSO) generally represent the same sort of object and show that both their radio spectra and optical emission can be explained by a single model which incorporates the e†ects produced by the interaction of a jet-driven nonthermal lobe with a dense interstellar medium.
Abstract: We adopt the view that the classes of active galactic nuclei (AGN) known variously as gigahertz peak spectrum (GPS) sources, compact steep-spectrum (CSS) sources, and compact symmetric objects (CSO) generally represent the same sort of object and show that both their radio spectra and optical emission can be explained by a single model which incorporates the e†ects produced by the interaction of a jet- driven nonthermal lobe with a dense interstellar medium. Following Begelman, we assume that these sources are young AGNs yr) in which the jets are propagating through an interstellar (ages ( 106 medium in which the hydrogen number density, decreases as a power law with radius, with the index n H dB 1 E2 and cm~3 at 1 kpc. The bow shock preceding the radio lobe is radiative at early n H D 10E100 times in such a dense environment, and the optical line emission produced by the shocked ISM and the associated photoionized precursor is proportional to the monochromatic radio power, consistent with the observational data of Gelderman & Whittle. The ionized gas surrounding the lobes has a signi-cant emission measure and a correspondingly high free-free opacity which is responsible for the 0.1E1 GHz peaks in the radio spectra. For jet energy Nuxes D1045E1046 ergs s~1, consistent with the observed radio powers of these objects, the crucial observed anticorrelation between peak frequency and size is readily recovered. The form of the radio spectra (power laws at high and low frequencies) indicate that the absorption is due to a cloudy/-lamentary medium with an approximately uniform distribution of opa- cities resulting from a combination of a two-phase interstellar medium, shock shredding of clouds impacted by the bow shock and thermal instabilities in the shocked ISM. The ionized medium envelop- ing the radio source also forms a Faraday screen which produces high rotation measure and substantial depolarization, readily accounting for another key property of this class of AGNs. Subject headings: radio continuum: galaxies

Journal ArticleDOI
TL;DR: In this article, it was shown that pulsation modes are driven through an opacity bump due to a local enhancement of the iron abundance in the envelopes of sdB stars, which provides a natural explanation for the instabilities found in the newly discovered class of pulsating sdB star.
Abstract: We present new calculations that strongly reinforce the idea—originally proposed by Charpinet et al.—that pulsation modes are driven through an opacity bump due to a local enhancement of the iron abundance in the envelopes of sdB stars. Our improved models incorporate nonuniform iron abundance distributions obtained through the condition of diffusive equilibrium between gravitational settling and radiative levitation. They also include special Rosseland opacity tables that take into account the large variations of the iron abundance about the cosmic value that are predicted by equilibrium radiative levitation theory. For representative models with M = 0.48 M☉ and log g = 5.8, we find strong instabilities for low-order radial and nonradial (p and f) pulsation modes in the range 36,500 K Teff 29,000 K. The four pulsating sdB stars currently known all have effective temperatures in that range. In addition, one of our models with Teff = 34,000 K has a band of unstable modes with periods in the range 116-195 s, in excellent agreement with those of the known pulsators. We therefore claim that our proposed iron bump mechanism provides a natural explanation for the instabilities found in the newly discovered class of pulsating sdB stars.

Journal ArticleDOI
16 May 1997-Science
TL;DR: The sensitivity of predictions of climate variations to assumptions about the microphysical processes has led to new approaches to atmospheric measurements and to heightened interest and progress in understanding the physical chemistry, radiative properties, and kinetics of small solid and liquid aqueous particles as discussed by the authors.
Abstract: Recent studies have shown that global radiative and hydrological fluxes are strongly linked to microphysical processes in clouds. The sensitivity of predictions of climate variations to assumptions about the microphysical processes has led to new approaches to atmospheric measurements and to heightened interest and progress in understanding the physical chemistry, radiative properties, and kinetics of small solid and liquid aqueous particles.

Journal ArticleDOI
TL;DR: In this article, a systematic formulation of various radiative transfer parameterizations is presented, including the absorption approximation (AA), δ-two-stream approximation (D2S), and D4S, in a consistent manner for thermal infrared flux calculations.
Abstract: A systematic formulation of various radiative transfer parameterizations is presented, including the absorption approximation (AA), δ-two-stream approximation (D2S), δ-four-stream approximation (D4S), and δ-two- and four-stream combination approximation (D2/4S), in a consistent manner for thermal infrared flux calculations. The D2/4S scheme uses a source function from the δ-two-stream approximation and evaluates intensities in the four-stream directions. A wide range of accuracy checks for monochromatic emissivity of a homogeneous layer and broadband heating rates and fluxes in nonhomogeneous atmospheres is performed with respect to the “exact” results computed from the δ-128-stream scheme for radiative transfer. The computer time required for the calculations using different radiative transfer parameterizations is compared. The results pertaining to the accuracy and efficiency of various radiative transfer approximations can be utilized to decide which approximate method is most appropriate for ...


Journal ArticleDOI
TL;DR: In this paper, the authors reported the modification of photoluminescence from conjugated polymers due to the proximity of metal films and showed that the presence of a metal film alters the radiative decay rate of an emitter via interference effects and also opens up an efficient nonradiative decay channel via energy transfer to the metal film.
Abstract: We report the modification of photoluminescence (PL) and electroluminescence (EL) from conjugated polymers due to the proximity of metal films. The presence of a metal film alters the radiative decay rate of an emitter via interference effects, and also opens up an efficient nonradiative decay channel via energy transfer to the metal film. We show that these effects lead to substantial changes in the PL and EL quantum efficiencies and the emission spectra of the polymers studied here [cyano derivatives of poly($p$-phenylenevinylene), PPV] as a function of the distance of the emitting dipoles from the metal film. We have measured the PL quantum efficiency directly using an integrating sphere, and found its distance dependence to be in good agreement with earlier theoretical predictions. Using the spectral dependence of the emission, we have been able to investigate the effect of interference on the radiative rate as a function of the wavelength and the distance between the emitter and the mirror. We compare our results with simulations of the radiative power of an oscillating dipole in a similar system. From our results we can determine the orientation of the dipoles in the polymer film, and the branching ratio that gives the fraction of absorbed photons leading to singlet excitons. We propose design rules for light-emitting diodes (LED's) and photovoltaic cells that optimize the effects of the metal film. By making optimum use of above effects we have substantially increased the EL quantum efficiencies of PPV/cyano-PPV double-layer LED's.

Journal ArticleDOI
TL;DR: In this article, it was shown that radiative torques on irregular dust grains, in addition to producing superthermal rotation, play a direct dynamical role in the alignment of interstellar dust with the local magnetic field.
Abstract: Radiative torques on irregular dust grains, in addition to producing superthermal rotation, play a direct dynamical role in the alignment of interstellar dust with the local magnetic field. The equations governing the orientation of spinning, precessing grains are derived; H2 formation torques and paramagnetic dissipation are included in the dynamics. Stationary solutions (constant alignment angle and spin rate) are found; these solutions may be stable ("attractors") or unstable ("repellors"). The equations of motion are numerically integrated for three exemplary irregular grain geometries, exposed to anisotropic radiation with the spectrum of interstellar starlight. The resulting "trajectory maps" are classified as "noncyclic," "semicyclic," or "cyclic," with examples of each given. We find that radiative torques result in rapid grain alignment, even in the absence of paramagnetic dissipation. It appears that radiative torques due to starlight can account for the observed alignment of interstellar grains with the Galactic magnetic field.

Journal ArticleDOI
TL;DR: The parallel implementation of the generalized stellar atmosphere and non-LTE (NLTE) radiative transfer computer program PHOENIX is described and the parallel algorithms developed for radiativeTransfer, spectral line opacity, and NLTE opacity and rate calculations are discussed.
Abstract: We describe the parallel implementation of our generalized stellar atmosphere and non-LTE (NLTE) radiative transfer computer program PHOENIX. We discuss the parallel algorithms we have developed for radiative transfer, spectral line opacity, and NLTE opacity and rate calculations. Our implementation uses a multiple instruction-multiple data design based on a relatively small number of MPI library calls. We report the results of test calculations on a number of different parallel computers and discuss the results of scalability tests.

Journal ArticleDOI
TL;DR: In this article, a new radiation code within a general circulation model was used to assess the direct solar and thermal radiative forcing by sulfate aerosol of anthropogenic origin and soot aerosol from fossil-fuel burning.
Abstract: A new radiation code within a general circulation model is used to assess the direct solar and thermal radiative forcing by sulfate aerosol of anthropogenic origin and soot aerosol from fossil-fuel burning. The radiative effects of different aerosol profiles, relative humidity parameterizations, chemical compositions, and internal and external mixtures of the two aerosol types are investigated. The contribution to the radiative forcing from cloudy sky regions is found to be negligible for sulfate aerosol; this is in contrast to recent studies where the cloudy sky contribution was estimated using a method in which the spatial correlation between cloud amount and sulfate burden was ignored. However, the radiative forcing due to fossil-fuel soot aerosol is enhanced in cloudy regions if soot aerosol exists within or above the cloud. The global solar radiative forcing due to sulfate aerosol is estimated to be −0.38 W m−2 and the global thermal radiative forcing is estimated to be +0.01 W m−2. The hemi...

Journal ArticleDOI
TL;DR: In this article, a multi-spectral radiative transfer code was used to investigate the important features of the direct radiative forcing of sulphate and soot aerosol; the indirect effect of both aerosol species is not investigated in this study.
Abstract: A newly developed multi-spectral radiative-transfer code is used to investigate the important features of the direct radiative forcing of sulphate and soot aerosol; the indirect effect of both aerosol species is not investigated in this study. the direct radiative forcing is presented for different surface albedos, solar zenith angles, relative humidities and aerosol vertical profiles together with effects upon the surface irradiance and long-wave radiative forcing. the effect of subgrid-scale variations in relative humidity are examined using idealized relative humidity distributions. the results show that subgrid-scale variations in relative humidity and the spatial correlation between cloud and areas of high relative humidity should be considered in future general-circulation model calculations of the direct forcing due to sulphate aerosol. A comparison of the direct forcing obtained by adjusting the surface albedo to that using the full multi-spectral column calculation is performed; the results indicate that recent estimates of the climate response to the direct forcing of sulphate may be too large. the contribution to the direct forcing from cloudy-sky regions appears to be negligible for sulphate aerosol but there is a considerable enhancement of the forcing due to soot aerosol if soot exists within or above clouds. These calculations show that a small amount of soot, relative to the sulphate mass loading, can cause a significant positive direct forcing and emphasize that the vertical profile of soot aerosol relative to cloud must be established to enable accurate assessment of the direct radiative effects of anthropogenic emissions of soot aerosol.

Journal ArticleDOI
TL;DR: In this paper, the average ice crystal size distributions are parameterized as functions of temperature and ice water content (IWC), based on observations in cirrus produced as outflows of deep convection made during the Central Equatorial Pacific Experiment (CEPEX), as the sum of a first-order gamma function describing ice crystals with melted equivalent diameters (Dm) less than 100 mm, and a lognormal function, describing larger ice crystals.
Abstract: Average ice crystal size distributions are parameterized as functions of temperature and ice water content (IWC), based on observations in cirrus produced as outflows of deep convection made during the Central Equatorial Pacific Experiment (CEPEX), as the sum of a first-order gamma function, describing ice crystals with melted equivalent diameters (Dm) less than 100 mm, and a lognormal function, describing larger ice crystals. The fit parameters are chosen to minimize the chi-squared function describing the difference between observed and parameterized distribution functions. The parameterization is mass conserving, accurately represents small ice crystals, and is easily integrable. The parameterization gives accurate estimates of mass, area, and number contained in different size ranges. The radiative properties estimated from midlatitude parameterizations are compared with those estimated from this parameterization using anomalous diffraction theory. As opposed to some previous studies, small crystals do not dominate the mass and radiative properties of cirrus. Comparison with midlatitude observations shows that size distribution shape can vary substantially depending on where, when, and how the cirrus is measured and on how it forms.

Journal ArticleDOI
TL;DR: It is proved formally that the extended KM model under perpendicular illumination is a good approximation for the cases of optically thick films that contain weakly or nonabsorbing particles.
Abstract: The description of optical properties of light-scattering materials has made extensive use of radiative transfer models. One of the most successful and simplest models is that of Kubelka and Munk (KM). With this model, optical properties of particulate films under diffuse illumination can be predicted from effective absorption and scattering coefficients of the material. We consider the applicability conditions of this kind of model. An extended KM model for the case of perpendicular collimated illumination is compared with results from a more general four-flux approach, and the differences between them are characterized in terms of a correction factor that depends on particle scattering and absorption, concentration of the scatterers, and film thickness. It is proved formally that the extended KM model under perpendicular illumination is a good approximation for the cases of optically thick films that contain weakly or nonabsorbing particles.

Journal ArticleDOI
TL;DR: In this article, the authors reexamine the empirical evidence for this scale break and explain it using fractal cloud models, Monte Carlo simulations, and a Green function approach to multiple scattering theory.
Abstract: Several studies have uncovered a break in the scaling properties of Landsat cloud scenes at nonabsorbing wavelengths. For scales greater than 200‐400 m, the wavenumber spectrum is approximately power law in k25/3, but from there down to the smallest observable scales (50‐100 m) follows another k2b law with b . 3. This implies very smooth radiance fields. The authors reexamine the empirical evidence for this scale break and explain it using fractal cloud models, Monte Carlo simulations, and a Green function approach to multiple scattering theory. In particular, the authors define the ‘‘radiative smoothing scale’’ and relate it to the characteristic scale of horizontal photon transport. The scale break was originally thought to occur at a scale commensurate with either the geometrical thickness Dz of the cloud, or with the ‘‘transport’’ mean free path lt 5 [(1 2 g)s]21, which incorporates the effect of forward scattering (s is extinction and g the asymmetry factor of the phase function). The smoothing scale is found to be approximately ltDz at cloud top; this is the prediction of diffusion ˇ theory which applies when (1 2 g)t 5D z / l t * 1( tis optical thickness). Since the scale break is a tangible effect of net horizontal radiative fluxes excited by the fluctuations of t, the smoothing scale sets an absolute lower bound on the range where one can neglect these fluxes and use plane-parallel theory locally, even for stratiform clouds. In particular, this constrains the retrieval of cloud properties from remotely sensed data. Finally, the characterization of horizontal photon transport suggests a new lidar technique for joint measurements of optical and geometrical thicknesses at about 0.5-km resolution.

Journal ArticleDOI
TL;DR: In this paper, the infrared absorption cross sections of CFC-11 and 16 other halogenated compounds have been measured and used in detailed line-by-line calculations to derive radiative forcing values.
Abstract: The infrared absorption cross sections of CFC-11 and 16 other halogenated compounds have been measured. These spectra were used in detailed line-by-line calculations to derive radiative forcing values. The radiative forcing values for 14 of these gases have not, to our knowledge, been previously reported in the literature. The accuracy of a computationally inexpensive narrowband scheme, which included the effect of clouds and stratospheric adjustment, was investigated. Global warming potentials are presented where atmospheric lifetimes are available. In light of the substantial disagreement in values for the forcing due to CFC-11 reported in the literature and its use as a standard to which other halogenated gases are often compared, we examined the sensitivity of CFC-11 forcing to a number of assumptions. We find that the uncertainties in the calculated value of the radiative forcing caused by neglect of the temperature and pressure dependence of the IR absorption spectrum are much smaller than those resulting from uncertainties in the absolute absorption cross sections or the vertical profile of CFC-11. Our best estimate is 0.285 W m−2 ppbv−1, which is 30% higher than the value adopted by the Intergovernmental Panel on Climate Change and is believed to be accurate to within about 10%. For the other gases represented here the lack of detailed knowledge of the likely vertical and horizontal distribution probably represents the most significant uncertainty in evaluating their radiative forcing.

Journal ArticleDOI
TL;DR: In this article, it was shown that the bulk motion of the converging flow is more efficient in upscattering photons than thermal Comptonization, provided that the electron temperature in the flow is of order a few keV or less.
Abstract: We study Compton upscattering of low-frequency photons in a converging flow of thermal plasma. The photons escape diffusively, and electron scattering is the dominant source of opacity. We solve the equation of radiative transfer in the case of spherical, steady state accretion into black holes numerically and approximately analytically. Unlike previous work on this subject, we consider the inner boundary at a finite radius, and this has a significant effect on the emergent spectrum. It is shown that the bulk motion of the converging flow is more efficient in upscattering photons than thermal Comptonization, provided that the electron temperature in the flow is of order a few keV or less. In this case, the spectrum observed at infinity consists of a soft component coming from input photons that escaped after a few scatterings without any significant energy change and of a power law that extends to high energies and is made of those photons that underwent significant upscattering. The luminosity of the power law is relatively small compared to that of the soft component. The more reflective the inner boundary is, the flatter the power-law spectrum becomes. The spectral energy power-law index for black hole accretion is always higher than 1, and it is approximately 1.5 for high accretion rates. This result tempts us to say that bulk motion Comptonization might be the mechanism behind the power-law spectra seen in black hole X-ray sources.

Journal ArticleDOI
TL;DR: In this paper, the lifetime of excess charge carriers in GaN epitaxial layers grown on sapphire at temperatures up to 300 K was measured by using picosecond time-resolved photoluminescence.
Abstract: By using picosecond time-resolved photoluminescence we have measured the lifetime of excess charge carriers in GaN epitaxial layers grown on sapphire at temperatures up to 300 K. The decay time turns out to be dominated by trapping processes at low excitation levels. The radiative lifetime derived from our data is dominated by free excitons at temperatures below 150 K, but also clearly shows the gradual thermal dissociation of excitons at higher temperatures. From our data, we are able to determine the free exciton binding energy and the free carrier radiative recombination coefficient. By combining these data with optical absorption data, we find the interband momentum matrix element and an estimate for the hole effective mass, which is much larger than previously thought.