Three-dimensional Radiative Transfer Effects Research Articles

Investigation of three-dimensional radiative transfer effects for UV–Vis satellite and ground-based observations of volcanic plumes

Abstract. We investigate effects of the three-dimensional (3D) structure of volcanic plumes on the retrieval results of satellite and ground-based UV–Vis observations. For the analysis of such measurements, 1D scenarios are usually assumed (the atmospheric properties only depend on altitude). While 1D assumptions are well suited for the analysis of many atmospheric phenomena, they are usually less appropriate for narrow trace gas plumes. For UV–Vis satellite instruments with large ground pixel sizes like the Global Ozone Monitoring Experiment-2 (GOME-2), the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) or the Ozone Monitoring Instrument (OMI), 3D effects are of minor importance, but usually these observations are not sensitive to small volcanic plumes. In contrast, observations of the TROPOspheric Monitoring Instrument (TROPOMI) on board Sentinel-5P have a much smaller ground pixel size (3.5 × 5.5 km2). Thus, on the one hand, TROPOMI can detect much smaller plumes than previous instruments. On the other hand, 3D effects become more important, because the TROPOMI ground pixel size is smaller than the height of the troposphere and also smaller than horizontal atmospheric photon path lengths in the UV–Vis spectral range. In this study we investigate the following 3D effects using Monte Carlo radiative transfer simulations: (1) the light-mixing effect caused by horizontal photon paths, (2) the saturation effect for strong SO2 absorption, (3) geometric effects related to slant illumination and viewing angles and (4) plume side-effects related to slant illumination angles and photons reaching the sensor from the sides of volcanic plumes. The first two effects especially can lead to a strong and systematic underestimation of the true trace gas content if 1D retrievals are applied (more than 50 % for the light-mixing effect and up to 100 % for the saturation effect). Besides the atmospheric radiative transfer, the saturation effect also affects the spectral retrievals. Geometric effects have a weaker influence on the quantitative analyses but can lead to a spatial smearing of elevated plumes or even to virtual double plumes. Plume side-effects are small for short wavelengths but can become large for longer wavelengths (up to 100 % for slant viewing and illumination angles). For ground-based observations, most of the above-mentioned 3D effects are not important because of the narrow field of view (FOV) and the closer distance between the instrument and the volcanic plume. However, the light-mixing effect shows a similar strong dependence on the horizontal plume extension as for satellite observations and should be taken into account for the analysis of ground-based observations.

Three-dimensional radiative transfer effects on airborne and ground-based trace gas remote sensing

Abstract. Air mass factors (AMFs) are used in passive trace gas remote sensing for converting slant column densities (SCDs) to vertical column densities (VCDs). AMFs are traditionally computed with 1D radiative transfer models assuming horizontally homogeneous conditions. However, when observations are made with high spatial resolution in a heterogeneous atmosphere or above a heterogeneous surface, 3D effects may not be negligible. To study the importance of 3D effects on AMFs for different types of trace gas remote sensing, we implemented 1D-layer and 3D-box AMFs into the Monte carlo code for the phYSically correct Tracing of photons In Cloudy atmospheres (MYSTIC), a solver of the libRadtran radiative transfer model (RTM). The 3D-box AMF implementation is fully consistent with 1D-layer AMFs under horizontally homogeneous conditions and agrees very well (<5 % relative error) with 1D-layer AMFs computed by other RTMs for a wide range of scenarios. The 3D-box AMFs make it possible to visualize the 3D spatial distribution of the sensitivity of a trace gas observation, which we demonstrate with two examples. First, we computed 3D-box AMFs for ground-based multi-axis spectrometer (MAX-DOAS) observations for different viewing geometry and aerosol scenarios. The results illustrate how the sensitivity reduces with distance from the instrument and that a non-negligible part of the signal originates from outside the line of sight. Such information is invaluable for interpreting MAX-DOAS observations in heterogeneous environments such as urban areas. Second, 3D-box AMFs were used to generate synthetic nitrogen dioxide (NO2) SCDs for an airborne imaging spectrometer observing the NO2 plume emitted from a tall stack. The plume was imaged under different solar zenith angles and solar azimuth angles. To demonstrate the limitations of classical 1D-layer AMFs, VCDs were then computed assuming horizontal homogeneity. As a result, the imaged NO2 plume was shifted in space, which led to a strong underestimation of the total VCDs in the plume maximum and an underestimation of the integrated line densities that can be used for estimating emissions from NO2 images. The two examples demonstrate the importance of 3D effects for several types of ground-based and airborne remote sensing when the atmosphere cannot be assumed to be horizontally homogeneous, which is typically the case in the vicinity of emission sources or in cities.

Using polarimetric observations to detect and quantify the three-dimensional radiative transfer effects in passive satellite cloud property retrievals: Theoretical framework and feasibility study

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions. This work is supported by the grant CyberTraining: DSE: Cross-Training of Researchers in Computing, Applied Mathematics and Atmospheric Sciences using Advanced Cyber infrastructure Resources from the US National Science Foundation (grant no. OAC–1730250). The computations in this study were performed at the UMBC High Performance Computing Facility (HPCF). The facility is supported by the US National Science Foundation through the MRI program (grant nos. CNS0821258 and CNS-1228778) and the SCREMS program (grant no. DMS-0821311), with substantial support from UMBC. ZZ’s research is partly supported by NASA grant (80NSSC20K0130) to UMBC.

A three-dimensional parallel radiative transfer model for atmospheric heating rates for use in cloud resolving models—The TenStream solver

This paper presents a new method to compute three-dimensional heating rates in atmospheric models, in particular numerical weather prediction models and large eddy simulations. The radiative transfer in such models is usually calculated for each vertical column independent of its neighbouring columns. Earlier studies showed that the neglect of horizontal energy transport introduces significant errors at model grid spacings below 1km. To date, there is no method to calculate 3D heating rates which is fast enough to systematically study the effect of radiation on cloud evolution. Here, we present a new algorithm which provides a fast yet accurate approximation for realistic three-dimensional heating rates. The method extends the well-known one-dimensional two-stream theory to 10 streams in three dimensions. Special emphasis is laid on scalable parallelism and speed.It is found that the new solver significantly reduces the root mean square error for atmospheric heating and surface heating rates when compared to traditionally employed one-dimensional solvers. The TenStream solver reduces the relative root mean square error of heating rates by a factor of five when compared to the independent column approximation. In the case of a strato-cumulus cloud field and the solar zenith angle being 60°, the error was reduced from 178% to 31% and for a deep-convective cumulus cloud from 138% to 28%. The model described here will open the way to answer the question, if and how much three-dimensional radiative transfer effects indeed affect cloud development and precipitation.

PaNTICA: A Fast 3D Radiative Transfer Scheme to Calculate Surface Solar Irradiance for NWP and LES Models

AbstractThe resolution of numerical weather prediction models is constantly increasing, making it necessary to consider three-dimensional radiative transfer effects such as cloud shadows cast into neighboring grid cells and thus affecting radiative heating. For that purpose, fast approximations are needed since three-dimensional radiative transfer solvers are computationally far too expensive. For the solar spectral range, different approaches of how to consider three-dimensional effects were presented in the past—in particular, the tilted independent column approximation (TICA), which aims at improving the calculation of the direct radiation, and the nonlocal tilted independent column approximation (NTICA), which is used to additionally correct the diffuse radiation. Here a new version of NTICA is presented that—in contrast to earlier approaches—is applicable for a variety of cloud scenes and grid resolutions and for arbitrary solar zenith angles. This new parameterization for the diffuse irradiance is then applied to the two different TICA approaches and the results are compared with a full 3D Monte Carlo calculation. It is shown that both approaches strongly improve the calculation of radiation fluxes if the new parameterization for the diffuse irradiance—what the authors call “parameterized NTICA (paNTICA)”—is applied. It is found that the method in which TICA is only applied to direct radiation yields the better results. The studies show that consideration of three-dimensional effects is inevitable if higher model resolutions are used in the future. This paper proposes ways to consider these effects and, thus, to substantially reduce the errors made with one-dimensional radiative transfer solvers.

Shortwave Radiative Impacts from Aerosol Effects on Marine Shallow Cumuli

Abstract The net shortwave radiative impact of aerosol on simulations of two shallow marine cloud cases is investigated using a Monte Carlo radiative transfer model. For a shallow cumulus case, increased aerosol concentrations are associated not only with smaller droplet sizes but also reduced cloud fractions and cloud dimensions, a result of evaporation-induced mixing and a lack of precipitation. Three-dimensional radiative transfer (3DRT) effects alter the fluxes by 10%–20% from values calculated using the independent column approximation for these simulations. The first (Twomey) aerosol indirect effect is dominant but the decreased cloud fraction reduces the magnitude of the shortwave cloud forcing substantially. The 3DRT effects slightly decrease the sensitivity of the cloud albedo to changes in droplet size under an overhead sun for the two ranges of cloud liquid water paths examined, but not strongly so. A popular two-stream radiative transfer approximation to the cloud susceptibility overestimates the more directly calculated values for the low liquid-water-path clouds within pristine aerosol conditions by a factor of 2 despite performing well otherwise, suggesting caution in its application to the cloud albedos within broken cloud fields. An evaluation of the influence of cloud susceptibility and cloud fraction changes to a “domain” area-weighted cloud susceptibility found that the domain cloud albedo is more likely to increase under aerosol loading at intermediate aerosol concentrations than under the most pristine conditions, contrary to traditional expectations. The second simulation (cumulus penetrating into stratus) is characterized by higher cloud fractions and more precipitation. This case has two regimes: a clean, precipitating regime where cloud fraction increases with increasing aerosol, and a more polluted regime where cloud fraction decreases with increasing aerosol. For this case the domain-mean cloud albedo increases steadily with aerosol loading under clean conditions, but increases only slightly after the cloud coverage decreases. Three-dimensional radiative transfer effects are mostly negligible for this case. Both sets of simulations suggest that aerosol-induced cloud fraction changes must be considered in tandem with the Twomey effect for clouds of small dimensions when assessing the net radiative impact, because both effects are drop size dependent and radiatively significant.

Twisting Motions of Sunspot Penumbral Filaments

The penumbra of a sunspot is composed of numerous thin, radially extended, bright and dark filaments carrying outward gas flows (the Evershed flow). Using high-resolution images obtained by the Solar Optical Telescope aboard the solar physics satellite Hinode, we discovered a number of penumbral bright filaments revealing twisting motions about their axes. These twisting motions are observed only in penumbrae located in the direction perpendicular to the symmetry line connecting the sunspot center and the solar disk center, and the direction of the twist (that is, lateral motions of intensity fluctuation across filaments) is always from limb side to disk-center side. Thus, the twisting feature is not an actual twist or turn of filaments but a manifestation of dynamics of penumbral filaments with three-dimensional radiative transfer effects.

The Accuracy of Determining Three-Dimensional Radiative Transfer Effects in Cumulus Clouds Using Ground-Based Profiling Instruments

Abstract Three-dimensional radiative transfer calculations are accurate, though computationally expensive, if the spatial distribution of cloud properties is known. The difference between these calculations and those using the much less expensive independent column approximation is called the 3D radiative transfer effect. Assessing the magnitude of this effect in the real atmosphere requires that many realistic cloud fields be obtained, and profiling instruments such as ground-based radars may provide the best long-term observations of cloud structure. Cloud morphology can be inferred from a time series of vertical profiles obtained from profilers by converting time to horizontal distance with an advection velocity, although this restricts variability to two dimensions. This paper assesses the accuracy of estimates of the 3D effect in shallow cumulus clouds when cloud structure is inferred in this way. Large-eddy simulations provide full three-dimensional, time-evolving cloud fields, which are sampled every 10 s to provide a “radar’s eye view” of the same cloud fields. The 3D effect for shortwave surface fluxes is computed for both sets of fields using a broadband Monte Carlo radiative transfer model, and intermediate calculations are made to identify reasons why estimates of the 3D effect differ in these fields. The magnitude of the 3D effect is systematically underestimated in the two-dimensional cloud fields because there are fewer cloud edges that cause the effect, while the random error in hourly estimates is driven by the limited sample observed by the profiling instrument.