Atmospheric radiative transfer in multiple scattering conditions. Application to NWP models.

High spectral resolution radiative transfer (RT) codes are essential tools in the study of the radiative energy transfer in the Earth atmosphere and a support for the development of parameterizations for fast RT codes used in climate and weather prediction models. Cirrus clouds cover permanently 30% of the Earth's surface, representing an important contribution to the Earth-atmosphere radiation balance. The work has been focussed on the development of the RT model LBLMS. The model, widely tested in the infra-red spectral range, has been extended to the short wave spectrum and it has been used in comparison with airborne and satellite measurements to study the optical properties of cirrus clouds. A new database of single scattering properties has been developed for mid latitude cirrus clouds. Ice clouds are treated as a mixture of ice crystals with various habits. The optical properties of the mixture are tested in comparison to radiometric measurements in selected case studies. Finally, a parameterization of the mixture for application to weather prediction and global circulation models has been developed. The bulk optical properties of ice crystals are parameterized as functions of the effective dimension of measured particle size distributions that are representative of mid latitude cirrus clouds. Tests with the Limited Area Weather Prediction model COSMO have shown the impact of the new parameterization with respect to cirrus cloud optical properties based on ice spheres.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

A new, fast and flexible radiative transfer method for Venus general circulation models

We present a new radiation scheme for the Oxford Planetary Unified Model System for Venus, suitable for the solar and thermal bands. This new and fast radiative parameterization uses a different approach in the two main radiative wavelength bands: solar radiation (0.1-5.5 mu m) and thermal radiation (1.7-260 mu m). The solar radiation calculation is based on the delta-Eddington approximation (two-stream-type) with an adding layer method. For the thermal radiation case, a code based on an absorptivity/emissivity formulation is used. The new radiative transfer formulation implemented is intended to be computationally light, to allow its incorporation in 3D global circulation models, but still allowing for the calculation of the effect of atmospheric conditions on radiative fluxes. This will allow us to investigate the dynamical-radiative-microphysical feedbacks. The model flexibility can be also used to explore the uncertainties in the Venus atmosphere such as the optical properties in the deep atmosphere or cloud amount. The results of radiative cooling and heating rates and the global-mean radiative-convective equilibrium temperature profiles for different atmospheric conditions are presented and discussed. This new scheme works in an atmospheric column and can be easily implemented in 3D Venus global circulation models. (C) 2014 Elsevier Ltd. All rights reserved.

Radiative effects of aerosols at an urban location in southern India: Observations versus model

The radiative impact of aerosols is one of the largest sources of uncertainty in estimating anthropogenic climate perturbations. Here we have used independent ground-based radiometer measurements made simultaneously with comprehensive measurements of aerosol microphysical and optical properties at a highly populated urban site, Bangalore (13.02 degrees N, 77.6 degrees E) in southern India during a dedicated campaign during winter of 2004 and summer and pre-monsoon season of 2005. We have also used longer term measurements carried out at this site to present general features of aerosols over this region. The aerosol radiative impact assessments were made from direct measurements of ground reaching irradiance as well as by incorporating measured aerosol properties into a radiative transfer model. Large discrepancies were observed between measured and modeled (using radiative transfer models, which employed measured aerosol properties) radiative impacts. It appears that the presence of elevated aerosol layers and (or) inappropriate description of aerosol state of mixing are (is) responsible for the discrepancies. On a monthly scale reduction of surface irradiance due to the presence of aerosols (estimated using radiative flux measurements) varies from 30 to 65 W m(-2). The lowest values in surface radiative impact were observed during June when there is large reduction in aerosol as a consequence of monsoon rainfall. Large increase in aerosol-induced surface radiative impact was observed from winter to summer. Our investigations re-iterate the inadequacy of aerosol measurements at the surface alone and importance of representing column properties (using vertical profiles) accurately in order to assess aerosol-induced climate changes accurately. (C) 2010 Elsevier Ltd. All rights reserved.

Effect of wall emissivities on radiation heat transfer in glass tank forehearths

In an earlier work, we had proposed a two-band, non-grey radiative transfer model for heat transfer in forehearths with simultaneous optically thick and thin approximations for molten glass interiors and at boundaries. Here using the same model, the radiative interaction of the top-crown and bottom-refractory walls with interior layers of shallow molten glass is studied by varying the wall emissivities. The forehearth exit temperature profiles for higher wall emissivities (0.9) show better conditioning of the glass for white flint glasses (optically thin).

Longwave radiative forcing of Indian Ocean tropospheric aerosol

A spectrally resolved discrete-ordinates radiative transfer model is used to calculate the change in downwelling surface and top-of-the-atmosphere (TOA) outgoing longwave (3.9-500 mum) radiative fluxes induced by tropospheric aerosols of the type observed over the Indian Ocean during the Indian Ocean Experiment (INDOEX). Both external and internal aerosol mixtures were considered. Throughout the longwave, the aerosol volume extinction depends more strongly on relative humidity than in most of the shortwave (0.28-3.9 mum), implying that particle growth factors and realistic relative humidity profiles must be taken into account when modeling the longwave radiative effects of aerosols. A typical boundary layer aerosol loading, with a 500-nm optical depth of 0.3, will increase the downwelling longwave flux at the surface by 7.7 W m(-2) over the clean air case while decreasing the outgoing longwave radiation by 1.3 W m(-2). A more vertically extended aerosol loading, exhibiting a high opacity plume between 2 and 3 km above the surface and having a typical 500-nm optical depth of 0.7, will increase the downwelling longwave flux at the surface by 11.2 W m(-2) over the clean air case while decreasing the outgoing longwave radiation by 2.7 W m(-2). For a vertically extended aerosol profile, approximately 30% of the TOA radiative forcing comes from sea salt and approximately 60% of the forcing comes from the combination of sea salt and dust. The remaining forcing is from anthropogenic constituents. These results are for the external mixture. For an internal mixture, TOA longwave forcings can be up to a factor of two larger. Therefore, to complete our understanding of this region's longwave aerosol radiative properties, more detailed information is needed about aerosol mixing states. These longwave radiative effects partially offset the large shortwave aerosol radiative forcing and should be included in regional and global climate modeling simulations.

Implications of multiple scattering on the assessment of black carbon aerosol radiative forcing

The effects of radiative coupling between scattering and absorbing aerosols, in an external mixture, on the aerosol radiative forcing (ARF) due to black carbon (BC), its sensitivity to the composite aerosol loading and composition, and surface reflectance are investigated using radiative transfer model simulations. The ARF due to BC is found to depend significantly on the optical properties of the `neighboring' (non-BC) aerosol species. The scattering due to these species significantly increases the top of the atmospheric warming due to black carbon aerosols, and significant changes in the radiative forcing efficiency of BC. This is especially significant over dark surfaces (such as oceans), despite the ARF due to BC being higher over snow and land-surfaces. The spatial heterogeneity of this effect (coupling or multiple scattering by neighboring aerosol species) imposes large uncertainty in the estimation ARF due to BC aerosols, especially over the oceans. (C) 2014 Elsevier Ltd. All rights reserved.

Scaling of magneto-quantum-radiative hydrodynamic equations: From laser-produced plasmas to astrophysics

We introduce the equations of magneto-quantum-radiative hydrodynamics. By rewriting them in a dimensionless form, we obtain a set of parameters that describe scale-dependent ratios of characteristic hydrodynamic quantities. We discuss how these dimensionless parameters relate to the scaling between astrophysical observations and laboratory experiments.

Radiative hydrodynamic modelling and observations of the X-class solar flare on 2011 March 9

Aims. We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters.

Methods. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 s, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300 s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy.

Results. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density (ne ≤ 1015 cm-3) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

Intercomparison of radiative forcing calculations of stratospheric water vapour and contrails

Seven groups have participated in an intercomparison study of calculations of radiative forcing (RF) due to stratospheric water vapour (SWV) and contrails. A combination of detailed radiative transfer schemes and codes for global-scale calculations have been used, as well as a combination of idealized simulations and more realistic global-scale changes in stratospheric water vapour and contrails. Detailed line-by-line codes agree within about 15 % for longwave (LW) and shortwave (SW) RF, except in one case where the difference is 30 %. Since the LW and SW RF due to contrails and SWV changes are of opposite sign, the differences between the models seen in the individual LW and SW components can be either compensated or strengthened in the net RF, and thus in relative terms uncertainties are much larger for the net RF. Some of the models used for global-scale simulations of changes in SWV and contrails differ substantially in RF from the more detailed radiative transfer schemes. For the global-scale calculations we use a method of weighting the results to calculate a best estimate based on their performance compared to the more detailed radiative transfer schemes in the idealized simulations.

The importance of the diurnal and annual cycle of air traffic for contrail radiative forcing

Air traffic condensation trails, or contrails, are believed to have a net atmospheric warming effect(1), although one that is currently small compared to that induced by other sources of human emissions. However, the comparably large growth rate of air traffic requires an improved understanding of the resulting impact of aircraft radiative forcing on climate(2). Contrails have an effect on the Earth's energy balance similar to that of high thin ice clouds(3). Their trapping of outgoing longwave radiation emitted by the Earth and atmosphere (positive radiative forcing) is partly compensated by their reflection of incoming solar radiation (negative radiative forcing). On average, the longwave effect dominates and the net contrail radiative forcing is believed to be positive(1,2,4). Over daily and annual timescales, varying levels of air traffic, meteorological conditions, and solar insolation influence the net forcing effect of contrails. Here we determine the factors most important for contrail climate forcing using a sophisticated radiative transfer model(5,6) for a site in southeast England, located in the entrance to the North Atlantic flight corridor. We find that night-time flights during winter (December to February) are responsible for most of the contrail radiative forcing. Night flights account for only 25 per cent of daily air traffic, but contribute 60 to 80 per cent of the contrail forcing. Further, winter flights account for only 22 per cent of annual air traffic, but contribute half of the annual mean forcing. These results suggest that flight rescheduling could help to minimize the climate impact of aviation.