Radiative Transfer Model Calculations Research Articles

An improved OMI ozone profile research product version 2.0 with collection 4 L1b data and algorithm updates

Abstract. We describe the new and improved version 2 of the ozone profile research product from the Ozone Monitoring Instrument (OMI) on the Aura satellite. One of the major changes is to switch the OMI L1b data from collection 3 to the recent collection 4 as well as the accompanying auxiliary datasets. The algorithm details are updated on radiative transfer model calculation and measurement calibrations, along with the input changes in meteorological data, and with the use of a tropopause-based ozone profile climatology, an improved high-resolution solar reference spectrum, and a recent ozone absorption cross-section dataset. A super Gaussian is applied to better represent OMI slit functions instead of a normal Gaussian. The effect of slit function errors on the spectral residuals is further accounted for as pseudo-absorbers in the iterative fit process. The OMI irradiances are averaged into monthly composites to reduce noise uncertainties in OMI daily measurements and to cancel out the temporal variations of instrument characteristics that are common in both radiance and irradiance measurements, which was previously neglected due to use of climatological composites. The empirical soft calibration spectra are re-derived to be consistent with the updated implementations and derived annually to remove the time-varying systematic biases between measured and simulated radiances. The “common mode” correction spectra are derived from remaining residual spectra after soft calibration as a function of solar zenith angle. The common mode is included as a pseudo-absorber in the iterative fit process, which helps to reduce the discrepancies of ozone retrieval accuracy between lower and higher solar zenith angles and between nadir and off-nadir pixels. Validation with ozonesonde measurements demonstrates the improvements of ozone profile retrievals in the troposphere, especially around the tropopause. The retrieval quality of tropospheric column ozone is improved with respect to the seasonal consistency between winter and summer as well as the long-term consistency before and after the row-anomaly occurrence.

Tropospheric bromine monoxide vertical profiles retrieved across the Alaskan Arctic in springtime

Abstract. Reactive halogen chemistry in the springtime Arctic causes ozone depletion events and alters the rate of pollution processing. There are still many uncertainties regarding this chemistry, including the multiphase recycling of halogens and how sea ice impacts the source strength of reactive bromine. Adding to these uncertainties are the impacts of a rapidly warming Arctic. We present observations from the CHACHA (CHemistry in the Arctic: Clouds, Halogens, and Aerosols) field campaign based out of Utqiaġvik, Alaska, from mid-February to mid-April of 2022 to provide information on the vertical distribution of bromine monoxide (BrO), which is a tracer for reactive bromine chemistry. Data were gathered using the Heidelberg Airborne Imaging DOAS (differential optical absorption spectroscopy) Instrument (HAIDI) on the Purdue University Airborne Laboratory for Atmospheric Research (ALAR) and employing a unique sampling technique of vertically profiling the lower atmosphere with the aircraft via “porpoising” maneuvers. Observations from HAIDI were coupled to radiative transfer model calculations to retrieve mixing ratio profiles throughout the lower atmosphere (below 1000 m), with unprecedented vertical resolution (50 m) and total information gathered (average of 17.5 degrees of freedom) for this region. A cluster analysis was used to categorize 245 retrieved BrO mixing ratio vertical profiles into four common profile shapes. We often found the highest BrO mixing ratios at the Earth's surface with a mean of nearly 30 pmol mol−1 in the lowest 50 m, indicating an important role for multiphase chemistry on the snowpack in reactive bromine production. Most lofted-BrO profiles corresponded with an aerosol profile that peaked at the same altitude (225 m above the ground), suggesting that BrO was maintained due to heterogeneous reactions on particle surfaces aloft during these profiles. A majority (11 of 15) of the identified lofted-BrO profiles occurred on a single day, 19 March 2022, over an area covering more than 24 000 km2, indicating that this was a large-scale lofted-BrO event. The clustered BrO mixing ratio profiles should be particularly useful for some MAX-DOAS (multi-axis DOAS) studies, where a priori BrO profiles and their uncertainties, used in optimal estimation inversion algorithms, are not often based on previous observations. Future MAX-DOAS studies (and past reanalyses) could rely on the profiles provided in this work to improve BrO retrievals.

Fast Calculation on Atmospheric Rayleigh Scattering Extinction Optical Depth for Day/Night Band Radiative Transfer Model With Periodic Lunar Illumination

AbstractRayleigh scattering effect inspired by tiny atmospheric molecules generates significant extinction at relatively short wavelength, such as Day/Night Band (DNB) covering visible (VIS, 0.4–0.76 μm) and near‐infrared (NIR, 0.76–1.0 μm) wavelengths. This salient extinction effect coupled with periodic changes in lunar illumination and wide DNB will further lead to increased complexity in nocturnal fast radiative transfer model (RTM) calculation of DNB. For these reasons, this study is aimed at developing a precise and robust fast calculation method on Rayleigh scattering extinction optical depth (RayOD) for pertinent fast RTM calculation of DNB. Firstly, we coupled the lunar spectral irradiance model and the parameterized Rayleigh scattering extinction calculation module into a RTM for forward simulations. Moreover, we find that the errors in DNB‐averaged RayOD calculated from the classical fast algorithm from RTTOV (Radiative Transfer for TOVS) of training the DNB instrument coefficients (about 20% of relative error in reflectance) and the parameterization method using the central wavelength of DNB (about 5%) are too large to be tolerable in fast RTM. Finally, to address this issue, we skillfully conceptualize a new equivalent optimal central wavelength algorithm to retrieve a dynamically wavelength in DNB using observation scattering and lunar phase angles for rapidly calculating DNB‐averaged RayOD. The relative errors in simulated reflectance from coupled RTM aforementioned are subtle nuances, which are about 10−5∼10−6 order of magnitude. Besides, this robust algorithm can also correct the relative error of about 1% at dark surface caused by periodic changes in lunar illumination.

The impact of cloudiness and cloud type on the atmospheric heating rate of black and brown carbon in the Po Valley

Abstract. We experimentally quantified the impact of cloud fraction and cloud type on the heating rate (HR) of black and brown carbon (HRBC and HRBrC). In particular, we examined in more detail the cloud effect on the HR detected in a previous study (Ferrero et al., 2018). High-time-resolution measurements of the aerosol absorption coefficient at multiple wavelengths were coupled with spectral measurements of the direct, diffuse and surface reflected irradiance and with lidar–ceilometer data during a field campaign in Milan, Po Valley (Italy). The experimental set-up allowed for a direct determination of the total HR (and its speciation: HRBC and HRBrC) in all-sky conditions (from clear-sky conditions to cloudy). The highest total HR values were found in the middle of winter (1.43 ± 0.05 K d−1), and the lowest were in spring (0.54 ± 0.02 K d−1). Overall, the HRBrC accounted for 13.7 ± 0.2 % of the total HR, with the BrC being characterized by an absorption Ångström exponent (AAE) of 3.49 ± 0.01. To investigate the role of clouds, sky conditions were classified in terms of cloudiness (fraction of the sky covered by clouds: oktas) and cloud type (stratus, St; cumulus, Cu; stratocumulus, Sc; altostratus, As; altocumulus, Ac; cirrus, Ci; and cirrocumulus–cirrostratus, Cc–Cs). During the campaign, clear-sky conditions were present 23 % of the time, with the remaining time (77 %) being characterized by cloudy conditions. The average cloudiness was 3.58 ± 0.04 oktas (highest in February at 4.56 ± 0.07 oktas and lowest in November at 2.91 ± 0.06 oktas). St clouds were mostly responsible for overcast conditions (7–8 oktas, frequency of 87 % and 96 %); Sc clouds dominated the intermediate cloudiness conditions (5–6 oktas, frequency of 47 % and 66 %); and the transition from Cc–Cs to Sc determined moderate cloudiness (3–4 oktas); finally, low cloudiness (1–2 oktas) was mostly dominated by Ci and Cu (frequency of 59 % and 40 %, respectively). HR measurements showed a constant decrease with increasing cloudiness of the atmosphere, enabling us to quantify for the first time the bias (in %) of the aerosol HR introduced by the simplified assumption of clear-sky conditions in radiative-transfer model calculations. Our results showed that the HR of light-absorbing aerosol was ∼ 20 %–30 % lower in low cloudiness (1–2 oktas) and up to 80 % lower in completely overcast conditions (i.e. 7–8 oktas) compared to clear-sky ones. This means that, in the simplified assumption of clear-sky conditions, the HR of light-absorbing aerosol can be largely overestimated (by 50 % in low cloudiness, 1–2 oktas, and up to 500 % in completely overcast conditions, 7–8 oktas). The impact of different cloud types on the HR was also investigated. Cirrus clouds were found to have a modest impact, decreasing the HRBC and HRBrC by −5 % at most. Cumulus clouds decreased the HRBC and HRBrC by −31 ± 12 % and −26 ± 7 %, respectively; cirrocumulus–cirrostratus clouds decreased the HRBC and HRBrC by −60 ± 8 % and −54 ± 4 %, which was comparable to the impact of altocumulus (−60 ± 6 % and −46 ± 4 %). A higher impact on the HRBC and HRBrC suppression was found for stratocumulus (−63 ± 6 % and −58 ± 4 %, respectively) and altostratus (−78 ± 5 % and −73 ± 4 %, respectively). The highest impact was associated with stratus, suppressing the HRBC and HRBrC by −85 ± 5 % and −83 ± 3 %, respectively. The presence of clouds caused a decrease of both the HRBC and HRBrC (normalized to the absorption coefficient of the respective species) of −11.8 ± 1.2 % and −12.6 ± 1.4 % per okta. This study highlights the need to take into account the role of both cloudiness and different cloud types when estimating the HR caused by both BC and BrC and in turn decrease the uncertainties associated with the quantification of their impact on the climate.

The Combined Effect of Ozone and Aerosols on Erythemal Irradiance in an Extremely Low Ozone Event during May 2020

In this study we focus on measurements and modeled UV index in the region of Athens, Greece, during a low ozone event. During the period of 12–19 May 2020, total ozone column (TOC) showed extremely low values, 35–55 Dobson Units (up to 15%) decrease from the climatic mean (being lower than the −2σ). This condition favors the increase of UV erythemal irradiance, since stratospheric ozone is the most important attenuator at the UVB spectral region. Simultaneously, an intrusion of Saharan dust aerosols in the region has masked a large part of the low ozone effect on UV irradiance. In order to investigate the event, we have used spectral solar irradiance measurements from the Precision Solar Radiometer (PSR), TOC from the Brewer spectrophotometer, and Radiative Transfer Model (RTM) calculations. Model calculations of the UV Index (UVI) showed an increase of ~30% compared to the long-term normal UVI due to the low TOC while at the same time and for particular days, aerosols masked this effect by ~20%. The RTM has been used to investigate the response in the UV spectral region of these variations at different solar zenith angles (SZAs). Spectra simulated with the RTM have been compared to measured ones and an average difference of ~2% was found. The study points out the importance of accurate measurements or forecasts of both ozone and aerosols when deriving UVI under unusual low ozone–high aerosol conditions.

A first comparison of TROPOMI aerosol layer height (ALH) to CALIOP data

Abstract. The TROPOspheric Monitoring Instrument (TROPOMI) level-2 aerosol layer height (ALH) product has now been released to the general public. This product is retrieved using TROPOMI's measurements of the oxygen A-band, radiative transfer model (RTM) calculations augmented by neural networks and an iterative optimal estimation technique. The TROPOMI ALH product will deliver ALH estimates over cloud-free scenes over the ocean and land that contain aerosols above a certain threshold of the measured UV aerosol index (UVAI) in the ultraviolet region. This paper provides background for the ALH product and explores its quality by comparing ALH estimates to similar quantities derived from spaceborne lidars observing the same scene. The spaceborne lidar chosen for this study is the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) mission, which flies in formation with NASA's A-train constellation since 2006 and is a proven source of data for studying ALHs. The influence of the surface and clouds is discussed, and the aspects of the TROPOMI ALH algorithm that will require future development efforts are highlighted. A case-by-case analysis of the data from the four selected cases (mostly around the Saharan region with approximately 800 co-located TROPOMI pixels and CALIOP profiles in June and December 2018) shows that ALHs retrieved from TROPOMI using the operational Sentinel-5 Precursor Level-2 ALH algorithm is lower than CALIOP aerosol extinction heights by approximately 0.5 km. Looking at data beyond these cases, it is clear that there is a significant difference when it comes to retrievals over land, where these differences can easily go over 1 km on average.

Infrared Properties of Asymptotic Giant Branch Stars in Our Galaxy and the Magellanic Clouds

Abstract We investigate the infrared properties of asymptotic giant branch (AGB) stars in our Galaxy and the Magellanic Clouds using various infrared observational data and theoretical models. We use catalogs for the sample of 4996 AGB stars in our Galaxy and about 39,000 AGB stars in the Magellanic Clouds from the available literature. For each object in the sample, we cross-identify the 2MASS, Wide-field Infrared Survey Explorer, and Spitzer counterparts. To compare the physical properties of O- and C-rich AGB stars in our Galaxy and the Magellanic Clouds, we present IR two-color diagrams (2CDs) using various photometric data. We perform radiative transfer model calculations for AGB stars using various possible parameters of central stars and dust shells. Using the dust opacity functions of amorphous silicate and carbon, the theoretical dust shell models can roughly reproduce the observations of AGB stars on various IR 2CDs. Compared with our Galaxy, we find that the Magellanic Clouds are deficient in AGB stars with thick dust shells. Compared with the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC) is more deficient in AGB stars with thick dust shells. This could be because the Magellanic Clouds are more metal-poor than our Galaxy and the LMC is more metal-rich than the SMC. We also present the IR properties of known pulsating variables. Investigating the magnitude distributions at mid-IR (MIR) bands for AGB stars in the Magellanic Clouds, we find that the SMC is more deficient in bright AGB stars at MIR bands compared with the LMC.

Fast scheme for determination of direct normal irradiance. Part I: New aerosol parameterization and performance assessment

A new aerosol direct transmittance parameterization is developed based on detailed radiative transfer model calculations for improving direct normal irradiance (DNI) for the solar energy forecast in the fast radiation scheme SUNFLUX. The modified SUNFLUX scheme is evaluated against observational datasets collected from the 14 Baseline Surface Radiation Network (BSRN) stations. The aerosol data from AERONET stations collocated with the BSRN network are used in the evaluations. The evaluation results have indicated that the accuracy of the calculated DNI using the new version of SUNFLUX has been improved significantly compared to the old version. The relative mean bias difference (MBD) and relative root mean square difference (RMSD) are −3.03% and 4.95%, respectively. The performance of the modified DNI model is also compared with those from the other five models. The results predicted by the current model are comparable with those from the best reference models. More importantly, the evaluation results show that the new aerosol scheme developed in this study can also be applied to other models to improve their performance.

A Deep Learning Trained Clear-Sky Mask Algorithm for VIIRS Radiometric Bias Assessment

Clear-sky mask (CSM) is a crucial influence on the calculating accuracy of the sensor radiometric biases for spectral bands of visible, infrared, and microwave regions. In this study, a fully connected deep neural network (FCDN) was proposed to generate CSM for the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard Suomi National Polar-Orbiting Partnership (S-NPP) and NOAA-20 satellites. The model, well-trained by S-NPP data, was used to generate both S-NPP and NOAA-20 CSMs for the independent data, and the results were validated against the biases between the sensor observations and Community Radiative Transfer Model (CRTM) calculations (O-M). The preliminary result shows that the FCDN-CSM model works well for identifying clear-sky pixels. Both O-M mean biases and standard deviations were comparable with the Advance Clear-Sky Processor over Ocean (ACSPO) and were significantly better than a prototype cloud mask (PCM) and the case without a clear-sky check. In addition, by replacing CRTM brightness temperatures (BTs) with the atmosphere air temperature and water vapor contents as input features, the FCDN-CSM exhibits its potential to generate fast and accurate VIIRS CSM onboard follow-up Joint Polar Satellite System (JPSS) satellites for sensor calibration and validation before the physics-based CSM is available.

Toward a Consistent Definition between Satellite and Model Clear-Sky Radiative Fluxes

AbstractA new method of determining clear-sky radiative fluxes from satellite observations for climate model evaluation is presented. The method consists of applying adjustment factors to existing satellite clear-sky broadband radiative fluxes that make the observed and simulated clear-sky flux definitions more consistent. The adjustment factors are determined from the difference between observation-based radiative transfer model calculations of monthly mean clear-sky fluxes obtained by ignoring clouds in the atmospheric column and by weighting hourly mean clear-sky fluxes with imager-based clear-area fractions. The global mean longwave (LW) adjustment factor is −2.2 W m−2at the top of the atmosphere and 2.7 W m−2at the surface. The LW adjustment factors are pronounced at high latitudes during winter and in regions with high upper-tropospheric humidity and cirrus cloud cover, such as over the west tropical Pacific, and the South Pacific and intertropical convergence zones. In the shortwave (SW), global mean adjustment is 0.5 W m−2at TOA and −1.9 W m−2at the surface. It is most pronounced over sea ice off of Antarctica and over heavy aerosol regions, such as eastern China. However, interannual variations in the regional SW and LW adjustment factors are small compared to those in cloud radiative effect. After applying the LW adjustment factors, differences in zonal mean cloud radiative effect between observations and climate models decrease markedly between 60°S and 60°N and poleward of 65°N. The largest regional improvements occur over the west tropical Pacific and Indian Oceans. In contrast, the impact of the SW adjustment factors is much smaller.

Light-absorbing impurities accelerating glacial melting in southeastern Tibetan Plateau

Deposition of light-absorbing particles on glacier surfaces poses a series of adverse impacts on the cryospheric environment, climate and human health. Broad attention of the scientific community has been paid on insoluble light-absorbing impurities (ILAIs) in snow and ice on glaciers over the Tibetan Plateau (TP). However, systematic investigation of ILAIs in snowpack of glaciers on the TP is scarce. In this study, the properties and darkening effect of ILAIs in snowpack on glaciers are extensively investigated in the southeast of TP. Results show that ILAIs concentrations in multiple types of snow and ice samples were significantly different. Snowpit depths varied substantially from one profile to another during May and June 2016. The average concentrations of ILAIs in snowpits increase as snow melting progresses. Black carbon (BC) and dust cause snow albedo reduction more in snow with larger grain size Re. Based on a radiative transfer model calculation, the average albedo reduction induced by BC in the snowpack was 0.141 ± 0.02, and associated daily maximum radiative forcing (RF) was 72.97 ± 12.7 W m−2. BC is a controlling light-absorbing factor in snowpack and causes substantial albedo reduction and thus the associated daily maximum RF. The maximum reduction of snow cover duration was 4.56 ± 0.71 days caused by BC and dust in snowpack in southeastern TP. The average mass absorption cross-section (MAC) of BC from multiple snowpits was 3.26 ± 0.46 m2 g−1, which represents a typical value of MAC in snow on glaciers, but it is type-dependent of snow/ice samples. Tropospheric aerosols vertically extended up to 8 km over the TP and its surrounding areas, which indicates the transport of aerosols from remote sources through elevated pathways. A large amount of carbon stored in the brittle glaciers can be potentially released with meltwater runoff under a warming climate. This study provides a new insight for investigating carbonaceous and light-absorbing particles in glacierization areas.

Analysis of Water Vapor Absorption in the Far‐Infrared and Submillimeter Regions Using Surface Radiometric Measurements From Extremely Dry Locations

AbstractThe second Radiative Heating in Underexplored Bands Campaign (RHUBC‐II) was conducted in 2009 by the U.S. Department of Energy Atmospheric Radiation Measurement program to improve water vapor spectroscopy in the far‐infrared spectral region. RHUBC‐II was located in an extremely dry region of Chile to ensure very low opacities in this spectral region. Spectrally resolved measurements by a far‐infrared spectrometer and a submillimeter interferometer from RHUBC‐II are compared with line‐by‐line radiative transfer model calculations. Water vapor amounts and temperatures used in the calculations come from collocated radiosondes, with extensive adjustments to correct for issues due to the campaign's dry conditions and mountainous terrain. A reanalysis is also performed of far‐infrared measurements taken at the Atmospheric Radiation Measurement North Slope of Alaska site before and during the first RHUBC campaign. These analyses determine that differences between the measurements and model calculations using existing spectroscopic parameters are significant in the far‐infrared and submillimeter regions, leading to the derivation of improved water vapor continuum absorption coefficients and air‐broadened widths of 74 water vapor lines. The foreign continuum is increased by more than 50% in part of the far‐infrared and the widths of more than 20 lines are changed by more than 10%. The uncertainty in the foreign continuum coefficients is estimated as greater than 20% in some spectral regions, primarily a consequence of the uncertainty in the specification of water vapor. The improved far‐infrared spectroscopic parameters have a notable impact on calculated spectral radiances and a modest impact on broadband radiative fluxes and heating rates.

Assessment of the role of the Pacific Ocean in present climate changes

Hydrothermodynamic processes in the atmosphere–ocean system played in favour of global warming slowdown in 1998–2014 were studied in this work. On the base of remote sensing and reanalysis data, close relationships between total global and regional column water vapour, terrestrial wind speed and temperature anomalies of upper layer water in tropical Pacific region were revealed. Increase of the wind speed in tropical Pacific has been observed since 1980 (linear trend ratio is –0.017 m · s–1/year). The most significant wind speed increase was in 1992–2013 (–0.025 m · s–1/year). During this period, the following phenomena were also observed: water temperature rise in upper layers of central and east equatorial Pacific regions by 0.024 K/year and accumulation of heat in the deeper layers of western Pacific north of the equator. These tendencies contributed to decrease in evaporation from the surface of the Pacific, which exerts considerable influence on the global mean water vapour content in the atmosphere with nearly 1-year lag (correlation coefficient is 0.88). Thus, average total column water vapour had been decreasing with average rate 0.12 mm/year until 2014. Atmospheric radiation transfer model calculations showed that decrease of water vapour content in atmospheric during 2001–2014 reduced the incoming part of Earth’s surface radiation balance by 0.93 W/m², which exceeds CO2-related increase in greenhouse warming by 11 times. Such behaviour of greenhouse gases concentrations could be the reason of decrease of winter temperature in Northern hemisphere. Summer temperatures continued to grow due to decrease in cloud optical depth in 35°N–70°N latitude zone and following radiation heating of the land surface.

No Surface Cooling over Antarctica from the Negative Greenhouse Effect Associated with Instantaneous Quadrupling of CO2 Concentrations

Abstract Over the highest elevations of Antarctica, during many months of the year, air near the surface is colder than in much of the overlying atmosphere. This unique feature of the Antarctic atmosphere has been shown to result in a negative greenhouse effect and a negative instantaneous radiative forcing at the top of the atmosphere , when carbon dioxide (CO2) concentrations are increased, and it has been suggested that this effect might play some role in te recent cooling trends observed over East Antarctica. Here, using fully coupled global climate model integrations, in addition to radiative transfer model calculations, the authors confirm the existence of such a negative over parts of Antarctica in response to an instantaneous quadrupling of CO2. However, it is also shown that the instantaneous radiative forcing at the tropopause is positive. Further, the negative lasts only a few days following the imposed perturbation, and rapidly disappears as the stratosphere cools in response to increased CO2. As a consequence, like the , the stratosphere-adjusted radiative forcing at the TOA is positive over all of Antarctica and, in the model presented herein, surface temperatures increase everywhere over that continent in response to quadrupled CO2. The results, therefore, clearly demonstrate that the curious negative instantaneous radiative forcing plays no role in the recently observed East Antarctic cooling.

Impact of Ice Cloud Microphysics on Satellite Cloud Retrievals and Broadband Flux Radiative Transfer Model Calculations

Ice cloud particles exhibit a range of shapes and sizes affecting a cloud’s single-scattering properties. Because they cannot be inferred from passive visible/infrared imager measurements, assumptions about the bulk single-scattering properties of ice clouds are fundamental to satellite cloud retrievals and broadband radiative flux calculations. To examine the sensitivity to ice particle model assumptions, three sets of models are used in satellite imager retrievals of ice cloud fraction, thermodynamic phase, optical depth, effective height, and particle size, and in top-of-atmosphere (TOA) and surface broadband radiative flux calculations. The three ice particle models include smooth hexagonal ice columns (SMOOTH), roughened hexagonal ice columns, and a two-habit model (THM) comprising an ensemble of hexagonal columns and 20-element aggregates. While the choice of ice particle model has a negligible impact on daytime cloud fraction and thermodynamic phase, the global mean ice cloud optical depth retrieved from THM is smaller than from SMOOTH by 2.3 (28%), and the regional root-mean-square difference (RMSD) is 2.8 (32%). Effective radii derived from THM are 3.9 μm (16%) smaller than SMOOTH values and the RMSD is 5.2 μm (21%). In contrast, the regional RMSD in TOA and surface flux between THM and SMOOTH is only 1% in the shortwave and 0.3% in the longwave when a consistent ice particle model is assumed in the cloud property retrievals and forward radiative transfer model calculations. Consequently, radiative fluxes derived using a consistent ice particle model assumption throughout provide a more robust reference for climate model evaluation compared to ice cloud property retrievals.

Process-Based Decomposition of the Decadal Climate Difference between 2002–13 and 1984–95

This study examines at the process level the climate difference between 2002–13 and 1984–95 in ERA-Interim. A linearized radiative transfer model is used to calculate the temperature change such that its thermal radiative cooling would balance the energy flux perturbation associated with the change of an individual process, without regard to what causes the change of the process in the first place. The global mean error of the offline radiative transfer model calculations is 0.09 K, which corresponds to the upper limit of the uncertainties from a single term in the decomposition analysis. The process-based decomposition indicates that the direct effect of the increase of CO2 (0.15 K) is the largest contributor to the global warming between the two periods (about 0.27 K). The second and third largest contributors are the cloud feedback (0.14 K) and the combined effect of the oceanic heat storage and evaporation terms (0.11 K), respectively. The largest warming associated with the oceanic heat storage term is found in the tropical Pacific and Indian Oceans, with relatively weaker warming over the tropical Atlantic Ocean. The increase in atmospheric moisture adds another 0.1 K to the global surface warming, but the enhancement in tropical convections acts to reduce the surface warming by 0.17 K. The ice-albedo and atmospheric dynamical feedbacks are the two leading factors responsible for the Arctic polar warming amplification (PWA). The increase of atmospheric water vapor over the Arctic region also contributes substantially to the Arctic PWA pattern.

Effect of temperature-dependent cross sections on O4 slant column density estimation by a space-borne UV–visible hyperspectral sensor

The sensitivities of oxygen dimer (O4) slant column densities (SCDs) were examined by applying temperature-dependent O4 cross sections using the radiative transfer model (RTM) calculation with the linearized pseudo-spherical vector discrete ordinate radiative transfer model. For the sensitivity study, we used a newly developed cross section database in place of the database used in the operational algorithm. Newly investigated O4 cross section databases for 203 K and 293 K were used for the radiance simulation by interpolating temperature for each atmospheric layer based on the vertical profile of standard atmosphere in the RTM. The effect of the temperature-dependent cross sections was a significant O4 SCD increase of 8.3% with dependence on satellite and solar viewing geometries. Furthermore, the O4 SCD generally increased by an estimated 3.9% based on the observation geometries of the Ozone Monitoring Instrument. For the long-term comparison, the O4 SCD estimated from the temperature-dependent cross sections corrects 20% of the total underestimation of O4 SCD between the observation and simulation. Although the surface pressure variation and background aerosol effect also correct the O4 SCD discrepancy, the effect of temperature-dependent cross sections was more important than the effects of surface pressure variation and background aerosols. Therefore, temperature dependence of the cross section in the RTM calculation is essential for the accurate simulation of O4 SCD.

The Zugspitze radiative closure experiment for quantifying water vapor absorption over the terrestrial and solar infrared – Part 2: Accurate calibration of high spectral-resolution infrared measurements of surface solar radiation

Abstract. Quantitative knowledge of water vapor absorption is crucial for accurate climate simulations. An open science question in this context concerns the strength of the water vapor continuum in the near infrared (NIR) at atmospheric temperatures, which is still to be quantified by measurements. This issue can be addressed with radiative closure experiments using solar absorption spectra. However, the spectra used for water vapor continuum quantification have to be radiometrically calibrated. We present for the first time a method that yields sufficient calibration accuracy for NIR water vapor continuum quantification in an atmospheric closure experiment. Our method combines the Langley method with spectral radiance measurements of a high-temperature blackbody calibration source (< 2000 K). The calibration scheme is demonstrated in the spectral range 2500 to 7800 cm−1, but minor modifications to the method enable calibration also throughout the remainder of the NIR spectral range. The resulting uncertainty (2σ) excluding the contribution due to inaccuracies in the extra-atmospheric solar spectrum (ESS) is below 1 % in window regions and up to 1.7 % within absorption bands. The overall radiometric accuracy of the calibration depends on the ESS uncertainty, on which at present no firm consensus has been reached in the NIR. However, as is shown in the companion publication Reichert and Sussmann (2016), ESS uncertainty is only of minor importance for the specific aim of this study, i.e., the quantification of the water vapor continuum in a closure experiment. The calibration uncertainty estimate is substantiated by the investigation of calibration self-consistency, which yields compatible results within the estimated errors for 91.1 % of the 2500 to 7800 cm−1 range. Additionally, a comparison of a set of calibrated spectra to radiative transfer model calculations yields consistent results within the estimated errors for 97.7 % of the spectral range.

A multi-wavelength classification method for polar stratospheric cloud types using infrared limb spectra

Abstract. The Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument on board the ESA Envisat satellite operated from July 2002 until April 2012. The infrared limb emission measurements represent a unique dataset of daytime and night-time observations of polar stratospheric clouds (PSCs) up to both poles. Cloud detection sensitivity is comparable to space-borne lidars, and it is possible to classify different cloud types from the spectral measurements in different atmospheric windows regions. Here we present a new infrared PSC classification scheme based on the combination of a well-established two-colour ratio method and multiple 2-D brightness temperature difference probability density functions. The method is a simple probabilistic classifier based on Bayes' theorem with a strong independence assumption. The method has been tested in conjunction with a database of radiative transfer model calculations of realistic PSC particle size distributions, geometries, and composition. The Bayesian classifier distinguishes between solid particles of ice and nitric acid trihydrate (NAT), as well as liquid droplets of super-cooled ternary solution (STS). The classification results are compared to coincident measurements from the space-borne lidar Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument over the temporal overlap of both satellite missions (June 2006–March 2012). Both datasets show a good agreement for the specific PSC classes, although the viewing geometries and the vertical and horizontal resolution are quite different. Discrepancies are observed between the CALIOP and the MIPAS ice class. The Bayesian classifier for MIPAS identifies substantially more ice clouds in the Southern Hemisphere polar vortex than CALIOP. This disagreement is attributed in part to the difference in the sensitivity on mixed-type clouds. Ice seems to dominate the spectral behaviour in the limb infrared spectra and may cause an overestimation in ice occurrence compared to the real fraction of ice within the PSC area in the polar vortex. The entire MIPAS measurement period was processed with the new classification approach. Examples like the detection of the Antarctic NAT belt during early winter, and its possible link to mountain wave events over the Antarctic Peninsula, which are observed by the Atmospheric Infrared Sounder (AIRS) instrument, highlight the importance of a climatology of 9 Southern Hemisphere and 10 Northern Hemisphere winters in total. The new dataset is valuable both for detailed process studies, and for comparisons with and improvements of the PSC parameterizations used in chemistry transport and climate models.

Dust Disks Around Young Stellar Objects

To reproduce the spectral energy distributions (SEDs) of young stellar objects (YSOs), we perform radiative transfer model calculations for the circumstellar dust disks with various shapes and many dust species. For eight sample objects of T Tauri and Herbig Ae/Be stars, we compare the theoretical model SEDs with the observed SEDs described by the infrared space observatory and Spitzer space telescope spectral data. We use the model, CGPLUS, for a passive irradiated circumstellar dust disk with an inner hole and an inner rim for the eight sample YSOs. We present model parameters for the dust disk, which reproduce the observed SEDs. We find that the model requires a higher mass, luminosity, and temperature for the central star for the Herbig Ae/Be stars than those for the T Tauri stars. Generally, the outer radius, total mass, thickness, and rim height of the theoretical dust disk for the Herbig Ae/Be stars are larger than those for the T Tauri stars.