Path-tracing Monte Carlo Library for 3D Radiative Transfer in Highly Resolved Cloudy Atmospheres

Interactions between clouds and radiation are at the root of many difficulties in numerically predicting future weather and climate and in retrieving the state of the atmosphere from remote sensing observations. The large range of issues related to these interactions, and in particular to three-dimensional interactions, motivated the development of accurate radiative tools able to compute all types of radiative metrics, from monochromatic, local and directional observables, to integrated energetic quantities. In the continuity of this community effort, we propose here an open-source library for general use in Monte Carlo algorithms. This library is devoted to the acceleration of path-tracing in complex data, typically high-resolution large-domain grounds and clouds. The main algorithmic advances embedded in the library are those related to the construction and traversal of hierarchical grids accelerating the tracing of paths through heterogeneous fields in null-collision (maximum cross-section) algorithms. We show that with these hierarchical grids, the computing time is only weakly sensitivive to the refinement of the volumetric data. The library is tested with a rendering algorithm that produces synthetic images of cloud radiances. Two other examples are given as illustrations, that are respectively used to analyse the transmission of solar radiation under a cloud together with its sensitivity to an optical parameter, and to assess a parametrization of 3D radiative effects of clouds.Comment: Submitted to JAMES, revised and submitted again (this is v2

Combining visible and infrared radiometry and lidar data to test simulations in clear and ice cloud conditions

Measurements taken during the 2003 Pacific THORPEX Observing System Test (P-TOST) by the MODIS Airborne Simulator (MAS), the Scanning High-resolution Interferometer Sounder (S-HIS) and the Cloud Physics Lidar (CPL) are compared to simulations performed with a line-by-line and multiple scattering modeling methodology (LBLMS). Formerly used for infrared hyper-spectral data analysis, LBLMS has been extended to the visible and near infrared with the inclusion of surface bi-directional reflectance properties. A number of scenes are evaluated: two clear scenes, one with nadir geometry and one cross-track encompassing sun glint, and three cloudy scenes, all with nadir geometry. <br><br> CPL data is used to estimate the particulate optical depth at 532 nm for the clear and cloudy scenes and cloud upper and lower boundaries. Cloud optical depth is retrieved from S-HIS infrared window radiances, and it agrees with CPL values, to within natural variability. MAS data are simulated convolving high resolution radiances. The paper discusses the results of the comparisons for the clear and cloudy cases. LBLMS clear simulations agree with MAS data to within 20% in the shortwave (SW) and near infrared (NIR) spectrum and within 2 K in the infrared (IR) range. It is shown that cloudy sky simulations using cloud parameters retrieved from IR radiances systematically underestimate the measured radiance in the SW and NIR by nearly 50%, although the IR retrieved optical thickness agree with same measured by CPL. <br><br> MODIS radiances measured from Terra are also compared to LBLMS simulations in cloudy conditions, using retrieved cloud optical depth and effective radius from MODIS, to understand the origin for the observed discrepancies. It is shown that the simulations agree, to within natural variability, with measurements in selected MODIS SW bands. <br><br> The impact of the assumed particles size distribution and vertical profile of ice content on results is evaluated. Sensitivity is much smaller than differences between measured and simulated radiances in the SW and NIR. <br><br> The paper dwells on a possible explanation of these contradictory results, involving the phase function of ice particles in the shortwave

Retrieving Decadal Climate Change from Satellite Radiance Observations-A 100-year CO2 Doubling OSSE Demonstration.

Preparing for climate change depends on the observation and prediction of decadal trends of the environmental variables, which have a direct impact on the sustainability of resources affecting the quality of life on our planet. The NASA Climate Absolute Radiance and Refractivity Observatory (CLARREO) mission is proposed to provide climate quality benchmark spectral radiance observations for the purpose of determining the decadal trends of climate variables, and validating and improving the long-range climate model forecasts needed to prepare for the changing climate of the Earth. The CLARREO will serve as an in-orbit, absolute, radiometric standard for the cross-calibration of hyperspectral radiance spectra observed by the international system of polar operational satellite sounding sensors. Here, we demonstrate that the resulting accurately cross-calibrated polar satellite global infrared spectral radiance trends (e.g., from the Metop IASI instrument considered here) can be interpreted in terms of temperature and water vapor profile trends. This demonstration is performed using atmospheric state data generated for a 100-year period from 2000-2099, produced by a numerical climate model prediction that was forced by the doubling of the global average atmospheric CO2 over the 100-year period. The vertical profiles and spatial distribution of temperature decadal trends were successfully diagnosed by applying a linear regression profile retrieval algorithm to the simulated hyperspectral radiance spectra for the 100-year period. These results indicate that it is possible to detect decadal trends in atmospheric climate variables from high accuracy all-sky satellite infrared radiance spectra using the linear regression retrieval technique

Spectral Detection of Human Skin in VIS-SWIR Hyperspectral Imagery without Radiometric Calibration

Many spectral detection algorithms require precise ground truth measurements that are hand-selected in the image to apply radiometric calibration, converting image pixels into estimated reflectance vectors. That process is impractical for mobile, real-time hyperspectral target detection systems, which cannot empirically derive a pixel-to-reflectance relationship from objects in the image. Implementing automatic target recognition on high-speed snapshot hyperspectral cameras requires the ability to spectrally detect targets without performing radiometric calibration. This thesis demonstrates human skin detection on hyperspectral data collected at a high frame rate without using calibration panels, even as the illumination in the scene changes. Compared to an established skin detection method that requires calibration panels, the illumination-invariant methods in this thesis achieve nearly as good detection performance in sunny scenes and superior detection performance in cloudy scenes

The FORUM end-to-end simulator project: architecture and results

FORUM (Far-infrared Outgoing Radiation Understanding and Monitoring) will fly as the ninth ESA's Earth Explorer mission, and an end-to-end simulator (E2ES) has been developed as a support tool for the mission selection process and the subsequent development phases. The current status of the FORUM E2ES project is presented together with the characterization of the capabilities of a full physics retrieval code applied to FORUM data. We show how the instrument characteristics and the observed scene conditions impact on the spectrum measured by the instrument, accounting for the main sources of error related to the entire acquisition process, and the consequences on the retrieval algorithm. Both homogeneous and heterogeneous case studies are simulated in clear and cloudy conditions, validating the E2ES against appropriate well-established correlative codes. The performed tests show that the performance of the retrieval algorithm is compliant with the project requirements both in clear and cloudy conditions. The far-infrared (FIR) part of the FORUM spectrum is shown to be sensitive to surface emissivity, in dry atmospheric conditions, and to cirrus clouds, resulting in improved performance of the retrieval algorithm in these conditions. The retrieval errors increase with increasing the scene heterogeneity, both in terms of surface characteristics and in terms of fractional cloud cover of the scene

An improved retrieval of tropospheric nitrogen dioxide from GOME

We present a retrieval of tropospheric nitrogen dioxide (NO2) columns from the Global Ozone Monitoring Experiment (GOME) satellite instrument that improves in several ways over previous retrievals, especially in the accounting of Rayleigh and cloud scattering. Slant columns, which are directly fitted without low-pass filtering or spectral smoothing, are corrected for an artificial offset likely induced by spectral structure on the diffuser plate of the GOME instrument. The stratospheric column is determined from NO2 columns over the remote Pacific Ocean to minimize contamination from tropospheric NO2. The air mass factor (AMF) used to convert slant columns to vertical columns is calculated from the integral of the relative vertical NO2 distribution from a global 3-D model of tropospheric chemistry driven by assimilated meteorological data (Global Earth Observing System (GEOS)-CHEM), weighted by altitude-dependent scattering weights computed with a radiative transfer model (Linearized Discrete Ordinate Radiative Transfer), using local surface albedos determined from GOME observations at NO2 wavelengths. The AMF calculation accounts for cloud scattering using cloud fraction, cloud top pressure, and cloud optical thickness from a cloud retrieval algorithm (GOME Cloud Retrieval Algorithm). Over continental regions with high surface emissions, clouds decrease the AMF by 20–30% relative to clear sky. GOME is almost twice as sensitive to tropospheric NO2 columns over ocean than over land. Comparison of the retrieved tropospheric NO2 columns for July 1996 with GEOS-CHEM values tests both the retrieval and the nitrogen oxide radical (NOx) emissions inventories used in GEOS-CHEM. Retrieved tropospheric NO2 columns over the United States, where NOx emissions are particularly well known, are within 18% of GEOS-CHEM columns and are strongly spatially correlated (r = 0.78, n = 288, p < 0.005). Retrieved columns show more NO2 than GEOS-CHEM columns over the Transvaal region of South Africa and industrial regions of the northeast United States and Europe. They are lower over Houston, India, eastern Asia, and the biomass burning region of central Africa, possibly because of biases from absorbing aerosols