Wildfire smoke is known as a highly absorptive aerosol type in the shortwave wavelength range. The absorption of sunlight by optically thick smoke layers results in heating of the ambient air. This heating is translated into self-lofting of the smoke up to more than 1 km in altitude per day. This study aims for a detailed analysis of tropospheric and stratospheric smoke lofting rates based on simulations and observations. The main goal is to demonstrate that radiative heating of intense smoke plumes is capable of lofting them from the lower and middle free troposphere (injection heights) up to the tropopause without the need of pyrocumulonimbus (pyroCb) convection. The further subsequent ascent within the lower stratosphere (caused by self-lofting) is already well documented in the literature. Simulations of absorbed solar radiation by smoke particles and resulting heating rates, which are then converted into lofting rates, are conducted by using the ECRAD (European Centre for Medium-Range Weather Forecasts Radiation) scheme. As input parameters thermodynamic profiles from CAMS (Copernicus Atmosphere Monitoring Service) reanalysis data, aerosol profiles from ground-based lidar observations, radiosonde potential temperature profiles, CALIOP (Cloud–Aerosol Lidar with Orthogonal Polarization) aerosol measurements, and MODIS (Moderate Resolution Imaging Spectroradiometer) aerosol optical depth retrievals were used. The sensitivity analysis revealed that the lofting rate strongly depends on aerosol optical thickness (AOT), layer depth, layer height, and black carbon (BC) fraction. We also looked at the influence of different meteorological parameters such as cloudiness, relative humidity, and potential temperature gradient. To demonstrate the applicability of our self-lofting model, we compared our simulations with the lofting processes in the stratosphere observed with CALIOP after major pyroCb events (Canadian fires in 2017, Australian fires in 2019–2020). We analyzed long-term CALIOP observations of smoke layers and plumes evolving in the UTLS (upper troposphere and lower stratosphere) height region over Siberia and the adjacent Arctic Ocean during the summer season of 2019. Our results indicate that self-lofting contributed to the vertical transport of smoke. We hypothesize that the formation of a near-tropopause aerosol layer, observed with CALIOP, was the result of self-lofting processes because this is in line with the simulations. Furthermore, Raman-lidar-based aerosol typing (in Leipzig and the High Arctic) clearly indicated the dominance of smoke in the UTLS aerosol layer since August 2019, most probably also the result of smoke self-lofting.

Cyprus plans to drastically increase the share of renewable energy sources from 13.9% in 2020 to 22.9% in 2030. Solar energy can play a key role in the effort to fulfil this goal. The potential for production of solar energy over the island is much higher than most of European territory because of the low latitude of the island and the nearly cloudless summers. In this study, high quality and fine resolution satellite retrievals of aerosols and dust, from the newly developed MIDAS climatology, and information for clouds from CM SAF are used in order to quantify the effects of aerosols, dust, and clouds on the levels of surface solar radiation for 2004–2017 and the corresponding financial loss for different types of installations for the production of solar energy. Surface solar radiation climatology has also been developed based on the above information. Ground-based measurements were also incorporated to study the contribution of different species to the aerosol mixture and the effects of day-to-day variability of aerosols on SSR. Aerosols attenuate 5–10% of the annual global horizontal irradiation and 15–35% of the annual direct normal irradiation, while clouds attenuate 25–30% and 35–50% respectively. Dust is responsible for 30–50% of the overall attenuation by aerosols and is the main regulator of the variability of total aerosol. All-sky annual global horizontal irradiation increased significantly in the period of study by 2%, which was mainly attributed to changes in cloudiness.

The time series of global radiation observed by a dense network of 99 autonomous pyranometers during the HOPE campaign around Jülich, Germany, are investigated with a multiresolution analysis based on the maximum overlap discrete wavelet transform and the Haar wavelet. For different sky conditions, typical wavelet power spectra are calculated to quantify the timescale dependence of variability in global transmittance. Distinctly higher variability is observed at all frequencies in the power spectra of global transmittance under broken-cloud conditions compared to clear, cirrus, or overcast skies. The spatial autocorrelation function including its frequency dependence is determined to quantify the degree of similarity of two time series measurements as a function of their spatial separation. Distances ranging from 100 m to 10 km are considered, and a rapid decrease of the autocorrelation function is found with increasing frequency and distance. For frequencies above 1∕3 min⁻¹ and points separated by more than 1 km, variations in transmittance become completely uncorrelated. A method is introduced to estimate the deviation between a point measurement and a spatially averaged value for a surrounding domain, which takes into account domain size and averaging period, and is used to explore the representativeness of a single pyranometer observation for its surrounding region. Two distinct mechanisms are identified, which limit the representativeness; on the one hand, spatial averaging reduces variability and thus modifies the shape of the power spectrum. On the other hand, the correlation of variations of the spatially averaged field and a point measurement decreases rapidly with increasing temporal frequency. For a grid box of 10 km × 10 km and averaging periods of 1.5–3 h, the deviation of global transmittance between a point measurement and an area-averaged value depends on the prevailing sky conditions: 2.8 (clear), 1.8 (cirrus), 1.5 (overcast), and 4.2 % (broken clouds). The solar global radiation observed at a single station is found to deviate from the spatial average by as much as 14–23 (clear), 8–26 (cirrus), 4–23 (overcast), and 31–79 W m⁻² (broken clouds) from domain averages ranging from 1 km × 1 km to 10 km × 10 km in area.

The clear-sky radiative effect of aerosol–radiation interactions is of relevance for our understanding of the climate system. The influence of aerosol on the surface energy budget is of high interest for the renewable energy sector. In this study, the radiative effect is investigated in particular with respect to seasonal and regional variations for the region of Germany and the year 2015 at the surface and top of atmosphere using two complementary approaches. First, an ensemble of clear-sky models which explicitly consider aerosols is utilized to retrieve the aerosol optical depth and the surface direct radiative effect of aerosols by means of a clear-sky fitting technique. For this, short-wave broadband irradiance measurements in the absence of clouds are used as a basis. A clear-sky detection algorithm is used to identify cloud-free observations. Considered are measurements of the short-wave broadband global and diffuse horizontal irradiance with shaded and unshaded pyranometers at 25 stations across Germany within the observational network of the German Weather Service (DWD). The clear-sky models used are the Modified MAC model (MMAC), the Meteorological Radiation Model (MRM) v6.1, the Meteorological–Statistical solar radiation model (METSTAT), the European Solar Radiation Atlas (ESRA), Heliosat-1, the Center for Environment and Man solar radiation model (CEM), and the simplified Solis model. The definition of aerosol and atmospheric characteristics of the models are examined in detail for their suitability for this approach. Second, the radiative effect is estimated using explicit radiative transfer simulations with inputs on the meteorological state of the atmosphere, trace gases and aerosol from the Copernicus Atmosphere Monitoring Service (CAMS) reanalysis. The aerosol optical properties (aerosol optical depth, Ångström exponent, single scattering albedo and asymmetry parameter) are first evaluated with AERONET direct sun and inversion products. The largest inconsistency is found for the aerosol absorption, which is overestimated by about 0.03 or about 30 % by the CAMS reanalysis. Compared to the DWD observational network, the simulated global, direct and diffuse irradiances show reasonable agreement within the measurement uncertainty. The radiative kernel method is used to estimate the resulting uncertainty and bias of the simulated direct radiative effect. The uncertainty is estimated to −1.5 ± 7.7 and 0.6 ± 3.5 W m⁻² at the surface and top of atmosphere, respectively, while the annual-mean biases at the surface, top of atmosphere and total atmosphere are −10.6, −6.5 and 4.1 W m⁻², respectively. The retrieval of the aerosol radiative effect with the clear-sky models shows a high level of agreement with the radiative transfer simulations, with an RMSE of 5.8 W m⁻² and a correlation of 0.75. The annual mean of the REari at the surface for the 25 DWD stations shows a value of −12.8 ± 5 W m⁻² as the average over the clear-sky models, compared to −11 W m⁻² from the radiative transfer simulations. Since all models assume a fixed aerosol characterization, the annual cycle of the aerosol radiation effect cannot be reproduced. Out of this set of clear-sky models, the largest level of agreement is shown by the ESRA and MRM v6.1 models.

The modification of an existing cloud property retrieval scheme for the Spinning Enhanced Visible and Infrared Imager (SEVIRI) instrument on board the geostationary Meteosat satellites is described to utilize its high-resolution visible (HRV) channel for increasing the spatial resolution of its physical outputs. This results in products with a nadir spatial resolution of 1×1 km² compared to the standard 3×3 km² resolution offered by the narrowband channels. This improvement thus greatly reduces the resolution gap between current geostationary and polar-orbiting meteorological satellite imagers. In the first processing step, cloudiness is determined from the HRV observations by a threshold-based cloud masking algorithm. Subsequently, a linear model that links the 0.6 µm, 0.8 µm, and HRV reflectances provides a physical constraint to incorporate the spatial high-frequency component of the HRV observations into the retrieval of cloud optical depth. The implementation of the method is described, including the ancillary datasets used. It is demonstrated that the omission of high-frequency variations in the cloud-absorbing 1.6 µm channel results in comparatively large uncertainties in the retrieved cloud effective radius, likely due to the mismatch in channel resolutions. A newly developed downscaling scheme for the 1.6 µm reflectance is therefore applied to mitigate the effects of this scale mismatch. Benefits of the increased spatial resolution of the resulting SEVIRI products are demonstrated for three example applications: (i) for a convective cloud field, it is shown that significantly better agreement between the distributions of cloud optical depth retrieved from SEVIRI and from collocated MODIS observations is achieved. (ii) The temporal evolution of cloud properties for a growing convective storm at standard and HRV spatial resolutions are compared, illustrating an improved contrast in growth signatures resulting from the use of the HRV channel. (iii) An example of surface solar irradiance, determined from the retrieved cloud properties, is shown, for which the HRV channel helps to better capture the large spatiotemporal variability induced by convective clouds. These results suggest that incorporating the HRV channel into the retrieval has potential for improving Meteosat-based cloud products for several application domains.

The temperature of photovoltaic modules is modelled as a dynamic function of ambient temperature, shortwave and longwave irradiance and wind speed, in order to allow for a more accurate characterisation of their efficiency. A simple dynamic thermal model is developed by extending an existing parametric steady-state model using an exponential smoothing kernel to include the effect of the heat capacity of the system. The four parameters of the model are fitted to measured data from three photovoltaic systems in the Allgäu region in Germany using non-linear optimisation. The dynamic model reduces the root-mean-square error between measured and modelled module temperature to 1.58 K on average, compared to 3.03 K for the steady-state model, whereas the maximum instantaneous error is reduced from 20.02 to 6.58 K.