We present a new model based on a convolutional neural network (CNN) to predict daytime cloud cover (CC) from sky images captured by all‐sky cameras, which is called CNN‐CC. A total of 49,016 daytime sky images, recorded at different Spanish locations (Valladolid, La Palma, and Izaña) from two different all‐sky camera types, are manually classified into different CC (oktas) values by trained researchers. Subsequently, the images are randomly split into a training set and a test set to validate the model. The CC values predicted by the CNN‐CC model are compared with the observations made by trained people on the test set, which serve as reference. The predicted CC values closely match the reference values within 1 oktas in 99% of the cloud‐free and overcast cases. Moreover, this percentage is above 93% for the rest of partially cloudy cases. The mean bias error (MBE) and standard deviation (SD) of the differences between the predicted and reference CC values are calculated, resulting in oktas and oktas. The MBE and SD are also represented for different intervals of measured aerosol optical depth and Ångström exponent values, revealing that the performance of the CNN‐CC model does not depend on aerosol load or size. Once the model is validated, the CC obtained from a set of images captured every 5 min, from January 2018 to March 2022, at the Antarctic station of Marambio (Argentina) is compared against direct field observations of CC (not from images) taken at this location, which is not used in the training process. As a result, the model slightly underestimates the observations with an MBE of 0.3 oktas. The retrieved data are analyzed in detail. The monthly and annual CC values are calculated. Overcast conditions are the most frequent, accounting for 46.5% of all observations throughout the year, rising to 64.5% in January. The annual mean CC value at this location is 5.5 oktas, with a standard deviation of approximately 3.1 oktas. A similar analysis is conducted, separating data by hours, but no significant diurnal cycles are observed except for some isolated months.

Aerosol optical depth (AOD) at night-time has become a hot topic in recent years due to the development of new instruments recording accurate ground-based lunar irradiance measurements, and the development of calibration methods and extraterrestrial irradiance models adapted to lunar photometry. This study uses all daytime and night-time AOD data available at Valladolid (Spain) from October 2016 to March 2022 in order to analyze its behavior and the added contribution of night data. The annual, monthly and daily AOD evolution is studied comparing daytime and night-time values and checking the correlation between them. For this purpose, the daily averages are computed, showing an annual pattern, with low AOD values throughout the year (mean value of AOD at 440 nm: 0.122), where winter months have the lower and summer the higher values, as observed in previous studies. All these AOD values are modulated by frequent desert dust events over the Iberian Peninsula, with a strong influence on daily and monthly mean values in February and March, where the strongest desert outbreaks occurred. The added new data confirm these results and the good correlation between daytime and night-time data. Also, a complete daily evolution is shown, observing that AOD and Ångström exponent (AE) mean values vary by only ±0.02 in 24 h, with a maximum value at 06:00 UTC and minimum at 18:00 UTC for both parameters.

The lack of aerosol optical depth (AOD) data at night can be partially addressed through Moon photometer measurements or fully covered with star photometer observations. However, the limited availability and complexity of star photometers have motivated this study to use all-sky cameras to extract starlight signals and derive AOD at night using star photometry. For this purpose, eight all-sky cameras were configured and deployed in nine different locations to capture raw images with varying exposure times every 2 min during the night. This work proposes a novel methodology to extract the starlight signal from the raw data from all-sky cameras and convert it into AOD values. This process consists of the following steps: removing the background image, selecting the pixels, and extracting the signal for each star from a predefined list of 56 stars; performing in situ Langley calibration of the instruments and retrieving the total optical depth (TOD); calculating the effective wavelength for each camera channel; deriving the AOD by subtracting the gas contribution to TOD; and averaging, cloud-screening, and quality-assuring the AOD time series. The AOD time series obtained through this methodology are compared with independent AOD measurements from collocated Moon photometers in the nine locations. The obtained results show that the AOD values derived with the proposed method generally correlate with reference values, often achieving correlation coefficients (r) above 0.90. The AOD values retrieved using the cameras tend to overestimate the reference values by approximately 0.02 and exhibit a precision of around 0.03–0.04. The agreement between the two datasets varies with wavelength and decreases at high-latitude locations, likely due to the poorer performance of Langley calibration in these regions. AOD values align well with day-to-night transitions obtained by solar photometers, demonstrating their reliability. Despite the slight overestimation, the AOD values derived by this new method approximate the real values and provide coverage throughout the entire night, without requiring the presence of the Moon. Therefore, they serve to study and monitor the nocturnal evolution of AOD.

The present work proposes a new model based on a convolutional neural network (CNN) to retrieve solar shortwave (SW) irradiance via the estimation of the cloud modification factor (CMF) from daytime sky images captured by all-sky cameras; this model is named CNN-CMF. To this end, a total of 237,669 sky images paired with SW irradiance measurements obtained by using pyranometers were selected at the following three sites: Valladolid and Izaña, Spain, and Lindenberg, Germany. This dataset was randomly split into training and testing sets, with the latter excluded from the training model in order to validate it using the same locations. Subsequently, the test dataset was compared with the corresponding SW irradiance measurements obtained by the pyranometers in scatter density plots. The linear fit shows a high determination coefficient ($R^2$) of $0.99$. Statistical analyses based on the mean bias error (MBE) values and the standard deviation (SD) of the SW irradiance differences yield results close to $-2\%$ and $9\%$, respectively. The MBE indicates a slight underestimation of the CNN-CMF model compared to the measurement values. After its validation, model performance was evaluated at the Antarctic station of Marambio (Argentina), a location not used in the training process. A similar comparison between the model-predicted SW irradiance and pyranometer measurements yielded $R^2=0.95$, with an MBE of around $2\%$ and an SD of approximately $26\%$. Although the precision provided by the SD at the Marambio station is lower, the MBE shows that the model’s accuracy is similar to previous results but with a slight overestimation of the SW irradiance. Finally, the determination coefficient improved to $0.99$, and the MBE and SD are about $3\%$ and $11\%$, respectively, when the CNN-CMF model is used to estimate daily SW irradiation values.

This work introduces CAECENET, a new system capable of automatically retrieving columnar and vertically-resolved aerosol properties running the GRASP (Generalized Retrieval of Atmosphere and Surface Properties) algorithm using sun-sky photometer (aerosol optical depth, AOD; and sky radiance measurements) and ceilometer (range corrected signal; RCS) data as input. This method, so called GRASP_pac, is implemented in CAECENET, which assimilates sun-sky photometers data from CÆLIS database and ceilometer data from ICENET database (Iberian Ceilometer Network). CAECENET allows for continuous and near-real-time monitoring of both vertical and columnar aerosol properties. The main characteristics and workflow of CAECENET are explained in detail. This work also explores the potential of CAECENET to monitor and analyze the evolution of transported aerosol events on a regional scale by means of the distribution of CAECENET stations across the Iberian Peninsula. As an example, this paper analyzes, using the CAECENET products, the case of a Saharan dust outbreak that occurred between the 3rd and 5th of October 2022. This was an intense event, with AOD at 440 nm values around 0.5 in Madrid and Valladolid, and reaching 1.55 in Granada. Transport from the Canadian wildfires at the end of June 2023 is also studied. Despite the long-range transport of the smoke particles in this event, measured volume concentrations reached and surpassed 80 μm³/cm³ in some stations. The results obtained point to the utility of this CAECENET tool for analyzing changes in the height and speed of the event propagation, in the aerosol concentration, and how this affects the optical properties.

Volcanic emissions from the Tajogaite volcano, located on the Cumbre Vieja edifice on the island of La Palma (Canary Islands, Spain), caused significant public health and aviation disruptions throughout the eruption (19 September–13 December 2021, officially declared over on 25 December). Nonetheless, it is considered the most significant volcanic event in Europe over the past 75 years due to the substantial amount of SO2 released into the atmosphere. The Instituto Geográfico Nacional (IGN), the authority responsible for volcano surveillance in Spain, implemented extensive operational monitoring to track volcanic activity and to provide a robust estimation of the volcanic plume height using a video-surveillance network. In parallel, the State Meteorological Agency of Spain (AEMET), in partnership with other Spanish ACTRIS (Aerosol, Clouds, and Trace Gases Research Infrastructure) members and collaborating institutions, conducted an unprecedented instrumental deployment to evaluate the impacts of this volcanic event on atmospheric composition. This effort included a network of aerosol profilers surrounding the volcano. A total of four profiling instruments were installed on La Palma: one MPL-4B lidar and three ceilometers. Additionally, a pre-existing Raman lidar on the island contributed valuable data to this study. These efforts are undertaken due to the importance of monitoring volcanic plume height in terms of air quality (necessary for the implementation of effective civil protection policies), volcanic activity surveillance (for tracking and forecasting eruptive behaviour), and, from a scientific perspective, for improving our understanding of the climatic and radiative impacts of this type of aerosol. In this study, the eruptive process was characterised in terms of the altitude of the dispersive volcanic plume (hd), measured by both IGN and AEMET-ACTRIS, and the altitude of the eruptive column (hec), measured by IGN. Modulating factors such as seismicity and meteorological conditions were also analysed. The consistency between the two independent and complementary datasets (hd,IGN and hd,AEMET) was assessed throughout the eruption (mean difference of 258.6 m). Our results confirmed the existence of three distinct eruptive phases, encompassing a range of styles from Strombolian explosive to effusive activity. While these phases have been characterised in previous studies, the results of the present work provide complementary information and novel insights from an alternative observational approach, which may be of use in future volcanic crises and will be applied to operational surveillance during such events. A subsequent comparison of hd,AEMET with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aerosol layer height product (ALHCALIOP) revealed a systematic underestimation by the satellite product, with a mean difference of 392.2 m. Finally, the impact of using hec in estimating SO2 emissions from the NASA MSVOLSO2L4 satellite-based product was evaluated. When a fixed (standard) plume altitude of 8 km was used instead of the observed hec, the total SO2 emission was significantly underestimated by an average of 56.2 %, and by up to 84.7 %. These findings underscore the importance of accurately determining the volcanic plume height when deriving SO2 emissions from satellite data.