Feasibility of Ground-Based Sky-Camera HDR Imagery to Determine Solar Irradiance and Sky Radiance over Different Geometries and Sky Conditions

We propose a methodological approach to provide the accurate and calibrated measurements of sky radiance and broadband solar irradiance using the High Dynamic Range (HDR) images of a sky-camera. This approach is based on a detailed instrumental characterization of a SONA sky-camera in terms of image acquisition and processing, as well as geometric and radiometric calibrations. As a result, a 1 min time resolution database of geometrically and radiometrically calibrated HDR images has been created and has been available since February 2020, with daily updates. An extensive validation of our radiometric retrievals has been performed in all sky conditions. Our results show a very good agreement with the independent measurements of the AERONET almucantar for sky radiance and pyranometers for broadband retrievals. The SONA sky radiance shows a difference of an RMBD < 10% while the broadband diffuse radiation shows differences of 2% and 5% over a horizontal plane and arbitrarily oriented surfaces, respectively. These results support the developed methodology and allow us to glimpse the great potential of sky-cameras to carry out accurate measurements of sky radiance and solar radiation components. Thus, the remote sensing techniques described here will undoubtedly be of great help for solar and atmospheric research.

Assessment of influence of urban aerosol vertical profile on clear-sky diffuse radiance pattern

Aerosol particles spread in the atmosphere play an important role in solar light scattering and thus co-determine the sky radiance/luminance pattern as well as diffuse irradiances/illuminances at the ground. The particular influence is given by their optical properties and by their distribution in the atmosphere. The dependence of the aerosol extinction coefficient on altitude is usually described by the exponential law, which results from averaging of a large amount of aerosol realizations. This is also frequently the case of simulating of the solar diffuse radiance/luminance distribution over the sky, when it is based on solving the radiative transfer problem. However, the aerosol vertical profile can sometimes be significantly different from the exponential one. Mainly in the urban environment, the aerosol is often well-mixed within the atmospheric boundary layer, so its volume extinction coefficient is almost constant there. This work investigates how such different profiles affect the clear sky radiance pattern and consequently also the ground-based horizontal diffuse irradiance. The numerical simulations reveal that the discrepancies are negligible in practice. So it appears that the aerosol vertical distribution does not play any important role in sky radiance calculations and the standard exponential law is general enough to cover also various specific aerosol conditions.

Estimation of Aerosol Columnar Size Distribution from Spectral Extinction Data in Coastal and Maritime Environment

Aerosol columnar size distributions (SDs) are commonly provided by aerosol inversions based on measurements of both spectral extinction and sky radiance. These inversions developed for a fully clear sky offer few SDs for areas with abundant clouds. Here, we estimate SDs from spectral extinction data alone for cloudy coastal and maritime regions using aerosol refractive index (RI) obtained from chemical composition data. Our estimation involves finding volume and mean radius of lognormally distributed modes of an assumed bimodal size distribution through fitting of the spectral extinction data. We demonstrate that vertically integrated SDs obtained from aircraft measurements over a coastal site have distinct seasonal changes, and these changes are captured reasonably well by the estimated columnar SDs. We also demonstrate that similar seasonal changes occur at a maritime site, and columnar SDs retrieved from the combined extinction and sky radiance measurements are approximated quite well by their extinction only counterparts (correlation exceeds 0.9) during a 7-year period (2013–2019). The level of agreement between the estimated and retrieved SDs depends weakly on wavelength selection within a given spectral interval (roughly 0.4–1 µm). Since the extinction-based estimations can be performed frequently for partly cloudy skies, the number of periods where SDs can be found is greatly increased.

Quality Scoring of the Fengyun 4A Clear Sky Radiance Product

The Clear Sky Radiance (CSR) product has been widely used instead of Level 1 (L1) geostationary imager data in data assimilation for numerical weather prediction due to its many advantages concerning superobservation methodology. In this study, CSR was produced in two water vapor channels (channels 9 and channel 10, with wavelengths at 5.8–6.7 μm and 6.9–7.3 μm) of the Advanced Geostationary Radiation Imager aboard Fengyun 4A. The root mean square error (RMSE) between CSR observations and backgrounds was used as a quality flag and was predicted by cloud cover, standard deviation (STD), surface type, and elevation of a CSR field of view (FOV). Then, a centesimal scoring system based on the predicted RMSE was set to a CSR FOV that indicates its percentile point in the quality distribution of the whole FOV. Validations of the scoring system demonstrated that the biases of the predicted RMSE were small for all FOVs and that the score was consistent with the predicted RMSE, especially for FOVs with high scores. We suggest using this score for quality control (QC) to replace the QC of cloud cover, STD, and elevation of CSR, and we propose 40 points as the QC threshold for the two channels, above which the predicted RMSE of a CSR is superior to the RMSE of averaged clear-sky L1 data.

Exploring the Assimilation of ATMS Cloud Retrieval Products and its Benefits for CrIS All-sky Radiance Simulations

Aimed at improving all-sky Cross-track Infrared Sounder (CrIS) radiance assimilation, this study explores the benefits for CrIS all-sky radiance simulations, focusing on the accuracy of background cloud information, through assimilating cloud liquid water path (LWP), ice water path (IWP), and rain water path (RWP) data retrieved from the Advanced Technology Microwave Sounder (ATMS). The Community Radiative Transfer Model (CRTM), which considers cloud scattering and absorption processes, is applied to simulate CrIS radiances. The Gridpoint Statistical Interpolation ensemble-variational data assimilation (DA) is updated by incorporating ensemble covariances of hydrometeor variables and observation operators of LWP, IWP, and RWP. First, two DA experiments named DActrl and DAcwp are conducted with (DAcwp) and without (DActrl) assimilating ATMS LWP, IWP, and RWP data. Assimilating ATMS cloud retrieval data results in better spatial distributions of hydrometers for both a Meiyu rainfall case and a typhoon case. Analyses of DActrl and DAcwp are then used as input to the CRTM to generate CrIS all-sky radiance simulations SMallsky_DActrl and SMallsky_DAcwp, respectively. Improvements in the DAcwp analyses of hydrometeor variables are found to benefit CrIS radiance simulations, especially in cloudy regions. A long period of statistics reveals that the biases and standard deviations of all-sky observations minus simulations from SMallsky_DAcwp are notably smaller than those from SMallsky_DActrl. This pilot study suggests the potential benefit of combining the use of microwave cloud retrieval products for all-sky infrared DA.

Computing light pollution indicators for environmental assessment

Light pollution modelling and monitoring has traditionally used zenith sky brightness as its main indicator. Several other indicators (e.g. average sky radiance, horizontal irradiance, average sky radiance at given interval of zenith distances) may be more useful, both for general and for specific purposes of ecology studies, night sky and environmental monitoring. These indicators can be calculated after the whole sky radiance is known with sufficient angular detail. This means, for each site, to integrate the contribution in each direction of the sky of each light source in the radius of hundreds of km. This approach is extremely high time consuming if the mapping is desired for a large territory. Here we present a way to obtain maps of large territories for a large subset of useful indicators, bypassing the need to calculate first the radiance map of the whole sky in each site to obtain from it the desired indicator in that site. For each indicator, a point spread function (PSF) is calculated from the whole sky radiance maps generated by a single source at sufficiently dense number of distances from the observing site. If the PSF is transversally shift-invariant, i.e. if it depends only on the relative position of source and observer, then we can further speed up the map calculation via the use of fast Fourier-transform (FFT). We present here examples of maps for different indicators. Precise results can be calculated for any single site, taking into account the site and light sources altitudes, by means of specific inhomogeneous (spatially-variant) and anisotropic (non rotationally symmetric) PSFs.