Top-of-atmosphere albedo bias from neglecting three-dimensional cloud radiative effects

Author(s):  
Clare E. Singer ◽  
Ignacio Lopez-Gomez ◽  
Xiyue Zhang ◽  
Tapio Schneider

AbstractClouds cover on average nearly 70% of Earth’s surface and regulate the global albedo. The magnitude of the shortwave reflection by clouds depends on their location, optical properties, and three-dimensional (3D) structure. Due to computational limitations, Earth system models are unable to perform 3D radiative transfer calculations. Instead they make assumptions, including the independent column approximation (ICA), that neglect effects of 3D cloud morphology on albedo. We show how the resulting radiative flux bias (ICA-3D) depends on cloud morphology and solar zenith angle. We use high-resolution (20–100 m horizontal resolution) large-eddy simulations to produce realistic 3D cloud fields covering three dominant regimes of low-latitude clouds: shallow cumulus, marine stratocumulus, and deep convective cumulonimbus. A Monte Carlo code is used to run 3D and ICA broadband radiative transfer calculations; we calculate the top-of-atmosphere (TOA) reflected flux and surface irradiance biases as functions of solar zenith angle for these three cloud regimes. Finally, we use satellite observations of cloud water path (CWP) climatology, and the robust correlation between CWP and TOA flux bias in our LES sample, to roughly estimate the impact of neglecting 3D cloud radiative effects on a global scale. We find that the flux bias is largest at small zenith angles and for deeper clouds, while the albedo bias is most prominent for large zenith angles. In the tropics, the annual-mean shortwave radiative flux bias is estimated to be 3.1±1.6 W m−2, reaching as much as 6.5 W m−2 locally.

2021 ◽  
Author(s):  
Claudia Emde ◽  
Huan Yu ◽  
Arve Kylling ◽  
Michel van Roozendael ◽  
Kerstin Stebel ◽  
...  

Abstract. Retrievals of trace gas concentrations from satellite observations are mostly performed for clear regions or regions with low cloud coverage. However, even fully clear pixels can be affected by clouds in the vicinity, either by shadowing or by scattering of radiation from clouds in the clear region. Quantifying the error of retrieved trace gas concentrations due to cloud scattering is a difficult task. One possibility is to generate synthetic data by three-dimensional (3D) radiative transfer simulations using realistic 3D atmospheric input data, including 3D cloud structures. Retrieval algorithms may be applied on the synthetic data and comparison to the known input trace gas concentrations yields the retrieval error due to cloud scattering. In this paper we present a comprehensive synthetic dataset which has been generated using the Monte Carlo radiative transfer model MYSTIC. The dataset includes simulated spectra in two spectral ranges (400–500 nm and the O2A-band from 755–775 nm). Moreover it includes layer air mass factors (layer-AMF) calculated at 460 nm. All simulations are performed for a fixed background atmosphere for various sun positions, viewing directions and surface albedos. Two cloud setups are considered: The first includes simple box-clouds with various geometrical and optical thicknesses. This can be used to systematically investigate the sensitivity of the retrieval error on solar zenith angle, surface albedo and cloud parameters. Corresponding 1D simulations are also provided. The second includes realistic three-dimensional clouds from an ICON large eddy simulation (LES) for a region covering Germany and parts of surrounding countries. The scene includes cloud types typical for central Europe such as shallow cumulus, convective cloud cells, cirrus, and stratocumulus. This large dataset can be used to quantify the trace gas concentration retrieval error statistically. Along with the dataset the impact of horizontal photon transport on reflectance spectra and layer-AMFs is analyzed for the box-cloud scenarios. Moreover, the impact of 3D cloud scattering on the NO2 vertical column density (VCD) retrieval is presented for a specific LES case. We find that the retrieval error is largest in cloud shadow regions, where the NO2 VCD is underestimated by more than 20 %. The dataset is available for the scientific community to assess the behavior of trace gas retrieval algorithms and cloud correction schemes in cloud conditions with 3D structure.


2007 ◽  
Vol 85 (11) ◽  
pp. 1159-1176
Author(s):  
D A Degenstein ◽  
A E Bourassa ◽  
E J Llewellyn ◽  
N D Lloyd

A simple radiative transfer model is developed to calculate the contribution of sea-glint to limb radiance. It is shown that the absolute sea-glint signal peaks between 70° and 80° solar zenith angle. Sea-glint can contribute 10–15% of the total limb radiance at wavelengths greater than 600 nm, which is several times brighter than an equivalent 5% reflecting Lambertian ocean surface. A test case was identified over the Arabian Sea in October 2002 and the model results compared to limb observations from the Optical Spectrograph and Infra-Red Imaging System (OSIRIS) on-board the Odin satellite. PACS Nos.: 94.10.Gb, 93.85.+q, 42.68.Ay, 42.68.Mj


2015 ◽  
Vol 15 (6) ◽  
pp. 8609-8646
Author(s):  
H. G. Chan ◽  
M. D. King ◽  
M. M. Frey

Abstract. Snow photochemical processes drive production of chemical trace gases, including nitrogen oxides (NO and NO2), and HOx radicals in snowpacks which are then released to the lower atmosphere. Coupled atmosphere–snow modelling on global scales requires simple parameterisations of actinic flux in snow to reduce computational cost. The disagreement between a physical radiative transfer method and a method based upon the e-folding depth of light-in snow is evaluated. In particular for the photolysis of the nitrate anion (NO3-), the nitrite anion (NO2-) and hydrogen peroxide (H2O2) within snow and photolysis of gas-phase nitrogen dioxide (NO2) within the snowpack interstitial air are considered. The emission flux from the snowpack is estimated as the depth-integrated photolysis rate, v, calculated (a) explicitly with a physical radiative transfer model (TUV), vTUV and (b) with a simple parameterisation based on e-folding depth, vze. The evaluation is based upon the deviation of the ratio of depth-integrated photolysis rate determined by the two methods,vTUV/vze, from unity. The disagreement in depth-integrated photolysis rate between the RT model and e-folding depth parameterisation depends primarily on the photolysis action spectrum of chemical species, solar zenith angle and optical properties of the snowpack, (scattering cross-section and a weak dependence on light absorbing impurity (black carbon) and density). For photolysis of NO2, the NO2- anion, the NO3- anion and H2O2 the ratio vTUV/vze varies within the range of 0.82–1.35, 0.88–1.28 and 0.92–1.27 respectively. The e-folding depth parameterisation underestimates for small solar zenith angles and overestimates at solar zenith angles around 60°. A simple algorithm has been developed to improve the parameterisation which reduced the ratio vTUV/vze to 0.97–1.02, 0.99–1.02 and 0.99–1.03 for photolysis of NO2, the NO2- anion, the NO3- anion and H2O2 respectively. The e-folding depth parameterisation may give acceptable results for the photolysis of the NO3- anion and H2O2 in cold polar snow with large solar zenith angles, but can be improved by a correction based on solar zenith angle.


2016 ◽  
Vol 55 (1) ◽  
pp. 93-117 ◽  
Author(s):  
Ehrhard Raschke ◽  
Stefan Kinne ◽  
William B. Rossow ◽  
Paul W. Stackhouse ◽  
Martin Wild

AbstractThis study examines radiative flux distributions and local spread of values from three major observational datasets (CERES, ISCCP, and SRB) and compares them with results from climate modeling (CMIP3). Examinations of the spread and differences also differentiate among contributions from cloudy and clear-sky conditions. The spread among observational datasets is in large part caused by noncloud ancillary data. Average differences of at least 10 W m−2 each for clear-sky downward solar, upward solar, and upward infrared fluxes at the surface demonstrate via spatial difference patterns major differences in assumptions for atmospheric aerosol, solar surface albedo and surface temperature, and/or emittance in observational datasets. At the top of the atmosphere (TOA), observational datasets are less influenced by the ancillary data errors than at the surface. Comparisons of spatial radiative flux distributions at the TOA between observations and climate modeling indicate large deficiencies in the strength and distribution of model-simulated cloud radiative effects. Differences are largest for lower-altitude clouds over low-latitude oceans. Global modeling simulates stronger cloud radiative effects (CRE) by +30 W m−2 over trade wind cumulus regions, yet smaller CRE by about −30 W m−2 over (smaller in area) stratocumulus regions. At the surface, climate modeling simulates on average about 15 W m−2 smaller radiative net flux imbalances, as if climate modeling underestimates latent heat release (and precipitation). Relative to observational datasets, simulated surface net fluxes are particularly lower over oceanic trade wind regions (where global modeling tends to overestimate the radiative impact of clouds). Still, with the uncertainty in noncloud ancillary data, observational data do not establish a reliable reference.


2018 ◽  
Vol 18 (16) ◽  
pp. 12105-12121 ◽  
Author(s):  
Thomas Fauchez ◽  
Steven Platnick ◽  
Tamás Várnai ◽  
Kerry Meyer ◽  
Céline Cornet ◽  
...  

Abstract. In a context of global climate change, the understanding of the radiative role of clouds is crucial. On average, ice clouds such as cirrus have a significant positive radiative effect, but under some conditions the effect may be negative. However, many uncertainties remain regarding the role of ice clouds on Earth's radiative budget and in a changing climate. Global satellite observations are particularly well suited to monitoring clouds, retrieving their characteristics and inferring their radiative impact. To retrieve ice cloud properties (optical thickness and ice crystal effective size), current operational algorithms assume that each pixel of the observed scene is plane-parallel and homogeneous, and that there is no radiative connection between neighboring pixels. Yet these retrieval assumptions are far from accurate, as real radiative transfer is 3-D. This leads to the plane-parallel and homogeneous bias (PPHB) plus the independent pixel approximation bias (IPAB), which impacts both the estimation of top-of-the-atmosphere (TOA) radiation and the retrievals. An important factor that determines the impact of these assumptions is the sensor spatial resolution. High-spatial-resolution pixels can better represent cloud variability (low PPHB), but the radiative path through the cloud can involve many pixels (high IPAB). In contrast, low-spatial-resolution pixels poorly represent the cloud variability (high PPHB), but the radiation is better contained within the pixel field of view (low IPAB). In addition, the solar and viewing geometry (as well as cloud optical properties) can modulate the magnitude of the PPHB and IPAB. In this, Part II of our study, we simulate TOA 0.86 and 2.13 µm solar reflectances over a cirrus uncinus scene produced by the 3DCLOUD model. Then, 3-D radiative transfer simulations are performed with the 3DMCPOL code at spatial resolutions ranging from 50 m to 10 km, for 12 viewing geometries and nine solar geometries. It is found that, for simulated nadir observations taken at resolution higher than 2.5 km, horizontal radiation transport (HRT) dominates biases between 3-D and 1-D reflectance calculations, but these biases are mitigated by the side illumination and shadowing effects for off-zenith solar geometries. At resolutions coarser than 2.5 km, PPHB dominates. For off-nadir observations at resolutions higher than 2.5 km, the effect that we call THEAB (tilted and homogeneous extinction approximation bias) due to the oblique line of sight passing through many cloud columns contributes to a large increase of the reflectances, but 3-D radiative effects such as shadowing and side illumination for oblique Sun are also important. At resolutions coarser than 2.5 km, the PPHB is again the dominant effect. The magnitude and resolution dependence of PPHB and IPAB is very different for visible, near-infrared and shortwave infrared channels compared with the thermal infrared channels discussed in Part I of this study. The contrast of 3-D radiative effects between solar and thermal infrared channels may be a significant issue for retrieval techniques that simultaneously use radiative measurements across a wide range of solar reflectance and infrared wavelengths.


2015 ◽  
Vol 15 (14) ◽  
pp. 8147-8163 ◽  
Author(s):  
M. Schäfer ◽  
E. Bierwirth ◽  
A. Ehrlich ◽  
E. Jäkel ◽  
M. Wendisch

Abstract. Based on airborne spectral imaging observations, three-dimensional (3-D) radiative effects between Arctic boundary layer clouds and highly variable Arctic surfaces were identified and quantified. A method is presented to discriminate between sea ice and open water under cloudy conditions based on airborne nadir reflectivity γλ measurements in the visible spectral range. In cloudy cases the transition of γλ from open water to sea ice is not instantaneous but horizontally smoothed. In general, clouds reduce γλ above bright surfaces in the vicinity of open water, while γλ above open sea is enhanced. With the help of observations and 3-D radiative transfer simulations, this effect was quantified to range between 0 and 2200 m distance to the sea ice edge (for a dark-ocean albedo of αwater = 0.042 and a sea-ice albedo of αice = 0.91 at 645 nm wavelength). The affected distance Δ L was found to depend on both cloud and sea ice properties. For a low-level cloud at 0–200 m altitude, as observed during the Arctic field campaign VERtical Distribution of Ice in Arctic clouds (VERDI) in 2012, an increase in the cloud optical thickness τ from 1 to 10 leads to a decrease in Δ L from 600 to 250 m. An increase in the cloud base altitude or cloud geometrical thickness results in an increase in Δ L; for τ = 1/10 Δ L = 2200 m/1250 m in case of a cloud at 500–1000 m altitude. To quantify the effect for different shapes and sizes of ice floes, radiative transfer simulations were performed with various albedo fields (infinitely long straight ice edge, circular ice floes, squares, realistic ice floe field). The simulations show that Δ L increases with increasing radius of the ice floe and reaches maximum values for ice floes with radii larger than 6 km (500–1000 m cloud altitude), which matches the results found for an infinitely long, straight ice edge. Furthermore, the influence of these 3-D radiative effects on the retrieved cloud optical properties was investigated. The enhanced brightness of a dark pixel next to an ice edge results in uncertainties of up to 90 and 30 % in retrievals of τ and effective radius reff, respectively. With the help of Δ L, an estimate of the distance to the ice edge is given, where the retrieval uncertainties due to 3-D radiative effects are negligible.


2013 ◽  
Vol 13 (2) ◽  
pp. 5477-5507
Author(s):  
J. Tonttila ◽  
P. Räisänen ◽  
H. Järvinen

Abstract. A new method for parameterizing the subgrid variations of vertical velocity and cloud droplet number concentration (CDNC) is presented for GCMs. These parameterizations build on top of existing parameterizations that create stochastic subgrid cloud columns inside the GCM grid-cells, which can be employed by the Monte Carlo independent column approximation approach for radiative transfer. The new model version adds a description for vertical velocity in individual subgrid columns, which can be used to compute cloud activation and the subgrid distribution of the number of cloud droplets explicitly. This provides a consistent way for simulating the cloud radiative effects with two-moment cloud microphysical properties defined in subgrid-scale. The primary impact of the new parameterizations is to decrease the CDNC over polluted continents, while over the oceans the impact is smaller. This promotes changes in the global distribution of the cloud radiative effects and might thus have implications on model estimation of the indirect radiative effect of aerosols.


2005 ◽  
Vol 62 (6) ◽  
pp. 1989-2000 ◽  
Author(s):  
F. Di Giuseppe ◽  
A. M. Tompkins

Abstract Conflicting claims have been made concerning the magnitude of the bias in solar radiative transfer calculations when horizontal photon transport is neglected for deep convective scenarios. The difficulty of obtaining a realistic set of cloud scenes for situations of complex cloud geometry, while certain characteristics such as total cloud cover are systematically controlled, has hindered the attempt to reach a consensus. Here, a simple alternative approach is adopted. An ensemble of cloud scenes generated by a cloud resolving model are modified by an idealized function that progressively alters the cirrus anvil coverage without affecting the realism of the scene produced. Comparing three-dimensional radiative calculations with the independent column approximation for all cloud scenes, it is found that the bias in scene albedo can reach as much as 22% when the sun is overhead and 46% at low sun angles. The bias is an asymmetrical function of cloud cover with a maximum attained at cirrus anvil cloud cover of approximately 30%–40%. With a cloud cover of 15%, the bias is half its maximum value, while it is limited for coverage exceeding 80%. The position of the peak occurs at the cloud cover coinciding with the maximum number of independent clouds present in the scene. Increasing the cloud cover past this point produces a decrease in the number of isolated clouds because of cloud merging, with a consequential bias reduction. With this systematic documentation of the biases as a function of total cloud cover, it is possible to identify two contributions to the total error: the geometrical consequences of the effective cloud cover increase at low sun angles and the true 3D scattering effect of photons deviating from the original path direction. An attempt to account for the former geometrical contribution to the 1D bias is made by performing a simple correction technique, whereby the field is sheared by the tangent of the solar zenith angle. It is found that this greatly reduces the 1D biases at low sun angles. Because of the small aspect ratio of the cirrus cloud deck, the remaining bias contribution is small in magnitude and almost independent of solar zenith angle.


2005 ◽  
Vol 44 (6) ◽  
pp. 789-803 ◽  
Author(s):  
Jordi Badosa ◽  
Josep-Abel González ◽  
Josep Calbó ◽  
Michiel van Weele ◽  
Richard L. McKenzie

Abstract To perform a climatic analysis of the annual UV index (UVI) variations in Catalonia, Spain (northeast of the Iberian Peninsula), a new simple parameterization scheme is presented based on a multilayer radiative transfer model. The parameterization performs fast UVI calculations for a wide range of cloudless and snow-free situations and can be applied anywhere. The following parameters are considered: solar zenith angle, total ozone column, altitude, aerosol optical depth, and single-scattering albedo. A sensitivity analysis is presented to justify this choice with special attention to aerosol information. Comparisons with the base model show good agreement, most of all for the most common cases, giving an absolute error within ±0.2 in the UVI for a wide range of cases considered. Two tests are done to show the performance of the parameterization against UVI measurements. One uses data from a high-quality spectroradiometer from Lauder, New Zealand [45.04°S, 169.684°E, 370 m above mean sea level (MSL)], where there is a low presence of aerosols. The other uses data from a Robertson–Berger-type meter from Girona, Spain (41.97°N, 2.82°E, 100 m MSL), where there is more aerosol load and where it has been possible to study the effect of aerosol information on the model versus measurement comparison. The parameterization is applied to a climatic analysis of the annual UVI variation in Catalonia, showing the contributions of solar zenith angle, ozone, and aerosols. High-resolution seasonal maps of typical UV index values in Catalonia are presented.


2020 ◽  
Author(s):  
Christoph Braun ◽  
Aiko Voigt ◽  
Johannes Hörner ◽  
Joaquim G. Pinto

<p>Atmospheric general circulation models developed for the Earth system include comprehensive parameterizations of clouds. Applying them to exoplanet atmospheres provides an opportunity to advance understanding of clouds, atmosphere dynamics, and their coupling in the context of planetary climate dynamics and habitability.</p><p>Here, we study a deep-time extreme climate of Earth as an example of the cold limit of the habitable zone. Geological evidence indicates near-global ice cover during the Neoproterozoic (1000 – 541 Million years ago) associated with considerable hysteresis of atmospheric CO<sub>2</sub>. The Snowball Earth hypothesis provides a straightforward interpretation of Neoproterozoic proxies based on a runaway of the sea-ice albedo feedback. However, the Snowball Earth hypothesis relies on the existence of local habitats to explain the survival of photosynthetic marine species on an entirely ice-covered planet. The Jormungand hypothesis may resolve this issue by considering a weakening of the sea-ice albedo feedback by exposure of dark bare sea ice when sea ice enters the subtropics. This potentially allows the Earth system to stabilize in a climate state - the Jormungand state - with near-global ice cover. Around the equator, a narrow strip of ocean remains ice-free, where life would have easily survived during the pan-glaciations.</p><p>The weakening of the sea-ice albedo feedback is based on the change of the meridional structure of planetary albedo with a moving sea-ice edge. While previous work focused on the contribution of surface albedo to planetary albedo, we here focus on the impact of subtropical and tropical cloudiness on planetary albedo. Enhanced cloudiness generally weakens the sea-ice albedo feedback and thus decreases the climate sensitivity of the Jormungand state, i.e. it stabilizes the Jormungand state. We analyze the impact of cloudiness on the stability of the Jormungand state in the general circulation models CAM3 and ICON-AES with idealized aquaplanet setups. While CAM3 shows significant CO<sub>2</sub>-hysteresis of the Jormungand state, ICON-AES exhibits no stable Jormungand state. Consistently, CAM3 exhibits stronger cloudiness than ICON-AES, especially in the subtropics. An analysis with a one-dimensional energy balance model shows that the Jormungand hysteresis strongly depends on the sensitivity of the planetary albedo to an advance of sea ice into the subtropics. Accordingly, we demonstrate that the absence of cloud-radiative effects within vertical columns in the subtropics drastically decreases the Jormungand hysteresis in CAM3.</p><p>Overall, the magnitude of the Jormungand hysteresis is tightly linked to the representation of cloud-radiative effects in general circulation models. Our results highlight the important role of uncertainties associated with cloud-radiative effects for climate feedbacks on planet Earth in the context of extreme climates, such as they have occurred in Earth’s deep past or might be found on Earth-like planets. In consequence, this also stresses the need and challenges of accounting for adequate cloud modeling for planetary climates.</p>


Sign in / Sign up

Export Citation Format

Share Document