Accurate assessment of atmospheric nitrogen dioxide (NO$_2$) and sulfur dioxide (SO$_2$) is essential for understanding climate-air quality interactions, supporting environmental policy, and protecting public health. Traditional monitoring approaches face limitations: satellite observations provide broad spatial coverage but suffer from data gaps, while ground-based sensors offer high temporal resolution but limited spatial extent. To address these challenges, we propose PollutionNet, a Vision Transformer-based framework that integrates Sentinel-5P TROPOMI vertical column density (VCD) data with ground-level observations. By leveraging self-attention mechanisms, PollutionNet captures complex spatiotemporal dependencies that are often missed by conventional CNN and RNN models. Applied to Ireland (2020-2021), our case study demonstrates that PollutionNet achieves state-of-the-art performance (RMSE: 6.89 $μ$g/m$^3$ for NO$_2$, 4.49 $μ$g/m$^3$ for SO$_2$), reducing prediction errors by up to 14% compared to baseline models. Beyond accuracy gains, PollutionNet provides a scalable and data-efficient tool for applied climatology, enabling robust pollution assessments in regions with spars
Tornadoes have caused billions of dollars in damage and are one of the leading causes of weather-related deaths in the United States each year. Tornadoes are known to frequently form in the warm sector of extratropical cyclones (ETCs), yet relatively little research exists quantifying the climatological relationship between the two phenomena and their covariability. This work analyzes the climatology of F/EF1+ tornadoes relative to ETCs (ETCTORs) using historical databases of tornadoes and hourly ETC centers for 1980-2022. Most tornadoes (72\%) occurred broadly to the southeast of an ETC center within 2000km, with a median distance of approximately 500km. Of those tornadoes, 69\% occurred in outbreaks of 6+ tornadoes. The spatial and ETC-relative distributions are similar across all intensity levels. Through the seasonal cycle, tornadoes shift north and south along with ETCs and the jet stream, and they are much more strongly (weakly) associated with ETCs in the winter (76\%; summer 38\%) when the jet stream and ETCs are strongest (weakest). Finally, median tornado and tornado-producing ETC locations covary strongly with one another interannually. A long-term eastward trend in torn
This study presents a comprehensive climatological benchmarking of tropical cyclones (TCs) generated by AI-based global weather prediction models. Using all TC events from the North Atlantic and Western Pacific basins between 2020 and 2025, we assess the ability of two AI models (Pangu-Weather and Aurora) to reproduce observed TC track density, climatology of storm characteristics, and physical consistency with TC theory. By comparing AI-simulated TCs with ERA5 reanalysis, we benchmark the distributions of intensity, size, forward speed, and evaluate the model's ability to credibly simulate extratropical transition. Results show that both Pangu and Aurora perform well in reproducing storm track density, forward speed distribution, and outer size distribution. Aurora shows an improved performance in simulating storm intensity compared to Pangu, with less bias in the distribution of minimum central pressure and maximum wind speed. However, both models overestimate the distribution of storm inner size (radius of maximum winds), especially for extreme events. AI models capture the relative frequency and temporal evolution of extratropical transition patterns with reasonable accuracy. T
We propose an index to quantify and analyse the impact of climatological variability on the energy system at different timescales. We define the Climatological Renewable Energy Deviation Index (CREDI) as the cumulative anomaly of a renewable resource with respect to its climate over a specific time period of interest. For this we introduce the smooth, yet physical, hourly rolling window climatology that captures the expected hourly to yearly behaviour of renewable resources. We analyse the presented index at decadal, annual and (sub-)seasonal timescales for a sample region and discuss scientific and practical implications. CREDI is meant as an analytical tool for researchers and stakeholders to help them quantify, understand, and explain, the impact of energy-meteorological variability on future energy system. Improved understanding translates to better assessments of how renewable resources, and the associated risks for energy security, may fare in current and future climatological settings. The practical use of the index is in resource planning. For example transmission system operators may be able to adjust short-term planning to reduce adequacy issues before they occur or combi
Planetary surface habitability has so far been considered, in the main, upon a global scale. The increasing number of 3D modelling studies of (exo)planetary climate has highlighted the need for a more nuanced understanding of surface habitability. Using satellite-derived data of photosynthetic life to represent the observed surface habitability of modern Earth, we validate a set of climatologically-defined metrics previously used in habitability studies. The comparison finds that the metrics defined by surface temperature alone show spatial patterns of habitability distinct to those defined by aridity or water availability, with no metric able to completely replicate the observed habitability. We build upon these results to introduce a new metric defined by the observed thermal limits of modern Earth-based life, along with surface water fluxes as an analogue for water and nutrient availability. Furthermore, we pay attention to not only the thermal bounds of macroscopic complex life, but additionally the limits of microbial life which have been vital to the generation of Earth's biosignatures, thus expanding considerations of climatic habitability out of a historically binary defini
In this short review, we discuss how Earth's climatological and geological history and also how the shadows of galactic black holes might reveal our Universe's past evolution. Specifically we point out that a pressure singularity that occurred in our Universe's past might have left its imprint on Earth's geological and climatological history and on the shadows of cosmological black holes. Our approach is based on the fact that the $H_0$ tension problem may be resolved if some sort of abrupt physics change occurred in our Universe $70-150\,$Myrs ago, an abrupt change that deeply affected the Cepheid parameters. We review how such an abrupt physics change might have been caused in our Universe by a smooth passage of it through a pressure finite-time singularity. Such finite-time singularities might occur in modified gravity and specifically in $F(R)$ gravity, so we show how modified gravity might drive this type of evolution, without resorting to peculiar cosmic fluids or scalar fields. The presence of such a pressure singularity can distort the elliptic trajectories of bound objects in the Universe, causing possible geological and climatological changes on Earth, if its elliptic tra
Eigen techniques such as empirical orthogonal function (EOF) or coupled pattern (CP) / maximum covariance analysis have been frequently used for detecting patterns in multivariate climatological data sets. Recently, statistical methods originating from the theory of complex networks have been employed for the very same purpose of spatio-temporal analysis. This climate network (CN) analysis is usually based on the same set of similarity matrices as is used in classical EOF or CP analysis, e.g., the correlation matrix of a single climatological field or the cross-correlation matrix between two distinct climatological fields. In this study, formal relationships as well as conceptual differences between both eigen and network approaches are derived and illustrated using exemplary global precipitation, evaporation and surface air temperature data sets. These results allow to pinpoint that CN analysis can complement classical eigen techniques and provides additional information on the higher-order structure of statistical interrelationships in climatological data. Hence, CNs are a valuable supplement to the statistical toolbox of the climatologist, particularly for making sense out of ve
Over the past decades, temperature and density of the upper atmosphere show negative trends and decrease of the upper atmospheric temperature is attributed to the declining neutral density. Specifically, nitric oxide (NO) and carbon dioxide (CO2) govern thermospheric cooling at 5.3 and 15 micron, respectively. While a lot of efforts have focused on the CO2 effects on the long-term trends, relatively less attention has been paid to the impacts by NO, which responds to solar and geomagnetic activities dynamically. In this study, we investigate the role of NO in climatological global energy budget for the recent three solar cycles using the Global Ionosphere Thermosphere Model. From 1982 to 2013, the F10.7 and Ap indices showed a decadal decrease of ~8% and ~20%, respectively. By imposing temporal-varying F10.7 and Ap values in the simulations, we find a decadal change of -0.28x1011 W or -17.3% in total NO cooling power, which agrees well with that (-0.34x1011 W or -17.2%) from the empirical Thermosphere Climate Index derived from the TIMED/SABER data. Neutral density decreases by 10-20% at 200-450km and Tex decreases by 25.3 K per decade. The deduced-decadal change of NO cooling reac
The proliferation of data-driven models in weather and climate sciences has marked a significant paradigm shift, with advanced models demonstrating exceptional skill in medium-range forecasting. However, these models are often limited by long-term instabilities, climatological drift, and substantial computational costs during training and inference, restricting their broader application for climate studies. Addressing these limitations, Guan et al. (2024) introduced LUCIE, a lightweight, physically consistent climate emulator utilizing a Spherical Fourier Neural Operator (SFNO) architecture. This model is able to reproduce accurate long-term statistics including climatological mean and seasonal variability. However, LUCIE's native resolution (~300 km) is inadequate for detailed regional impact assessments. To overcome this limitation, we introduce a deep learning-based downscaling framework, leveraging probabilistic diffusion-based generative models with conditional and posterior sampling frameworks. These models downscale coarse LUCIE outputs to 25 km resolution. They are trained on approximately 14,000 ERA5 timesteps spanning 2000-2009 and evaluated on LUCIE predictions from 2010
MJO and SPV are prominent sources of subseasonal predictability in the Extratropics. With relevance for European weather it has been shown that the joint interaction of MJO and the SPV can modulate the preferred phase of the NAO and the occurrence of weather regimes. However, improving extended-range NWP at three-week lead times remain under-explored. This study investigates how MJO and SPV phases affect Greenland Blocking (GL) activity and integrates atmospheric state information into a neural network to enhance week-three weather regime activity forecasts. We define a weather regime activity metric using ECMWF reanalysis and reforecasts. In reanalyses we find increased GL activity following MJO phases 7,8 and 1, as well as weak SPV phases, indicating climatological windows of opportunity in line with previous studies. However, ECMWF forecast skill improves only in MJO phases 8 and 1 and weak SPV phases, identifying somewhat different model windows of opportunities. Next we explore using these findings in post-processing tools. Climatological forecasts based on MJO/SPV-NAO relationships provide a purely statistical approach to extended-range GL activity forecasting, independent of
The frequency of extratropical cyclones in East Asia, including those traveling along the Kuroshio off the south coast of Japan, maximizes climatologically in spring in harmony with local enhancement of precipitation. The springtime cyclone activity is of great socioeconomic importance for East Asian countries. However, mechanisms for the spring peak in the East Asian cyclone activity have been poorly understood. This study aims to unravel the mechanisms, focusing particularly on favorable conditions for relevant cyclogenesis. Through a composite analysis based on atmospheric reanalysis data, we show that cyclogenesis enhanced around the East China Sea under anomalously strengthened cyclonic wind shear and temperature gradient, in addition to enhanced moisture flux from the south, is important for the spring peak in the cyclone activity in East Asia. In spring, climatologically strengthened cyclonic shear north of the low-level jet axis and associated frequent atmospheric frontogenesis in southern China and the East China Sea serve as favorable background conditions for low-level cyclogenesis. We also demonstrate that climatologically enhanced diabatic heating around East Asia is p
An Ensemble of Data Assimilations (EDA) can provide valuable information on the analysis and short-range forecast uncertainties. The present ECMWF operational ocean analysis and reanalysis system, called ORAS5, produces an ensemble but does not exploit it for the specification of the background-error covariance matrix $\mathbf{B}$, a key component of the data assimilation system. In this article, we describe EDA developments for the ocean, which take advantage of the short-range forecast ensemble for specifying, in two distinct ways, parameters of a covariance model representation of $\mathbf{B}$. First, we generate a climatological ensemble over an extended period to produce seasonally varying climatological estimates of background-error variances and horizontal correlation length-scales. Second, on each assimilation cycle, we diagnose flow-dependent variances from the ensemble and blend them with the climatological estimates to form hybrid variances. We also use the ensemble to diagnose flow-dependent vertical correlation length-scales. We demonstrate for the Argo-rich period that this new, hybrid formulation of $\mathbf{B}$ results in a significant reduction of background errors
Seasonal forecast of Arctic sea ice concentration is key to mitigate the negative impact and assess potential opportunities posed by the rapid decline of sea ice coverage. Seasonal prediction systems based on climate models often show systematic biases and complex spatio-temporal errors that grow with the forecasts. Consequently, operational predictions are routinely bias corrected and calibrated using retrospective forecasts. For predictions of Arctic sea ice concentration, error corrections are mainly based on one-to-one post-processing methods including climatological mean or linear regression correction and, more recently, machine learning. Such deterministic adjustments are confined at best to the limited number of costly-to-run ensemble members of the raw forecast. However, decision-making requires proper quantification of uncertainty and likelihood of events, particularly of extremes. We introduce a probabilistic error correction framework based on a conditional Variational Autoencoder model to map the conditional distribution of observations given the biased model prediction. This method naturally allows for generating large ensembles of adjusted forecasts. We evaluate our
Long-range ensemble forecasts are typically verified as anomalies with respect to a lead-time dependent climatological mean to remove the influence of systematic biases. However, common methods for calculating anomalies result in statistical inconsistencies between forecast and verification anomalies, even for a perfectly reliable ensemble. It is important to account for these systematic effects when evaluating ensemble forecast systems, particularly when tuning a model to improve the reliability of forecast anomalies or when comparing spread-error diagnostics between systems with different reforecast periods. Here, we show that unbiased variances and spread-error ratios can be recovered by deriving estimators that are consistent with the values that would be achieved when calculating anomalies relative to the true, but unknown, climatological mean. An elegant alternative is to construct forecast climatologies separately for each member, which ensures that forecast and verification anomalies are defined relative to reference climatological means with the same sampling uncertainty. This alternative approach has no impact on forecast ensemble means but systematically modifies the tot
Accurate Subseasonal-to-Seasonal (S2S) forecasting is vital for decision-making in agriculture, energy production, and emergency management. However, it remains a challenging and underexplored problem due to the chaotic nature of the weather system. Recent data-driven studies have shown promising results, but their performance is limited by the inadequate incorporation of climate states and a model tendency to degrade, progressively losing fine-scale details and yielding over-smoothed forecasts. To overcome these limitations, we propose TianQuan-S2S, a global S2S forecasting model that integrates initial weather states with climatological means via incorporating climatology into patch embedding and enhancing variability capture through an uncertainty-augmented Transformer. Extensive experiments on the Earth Reanalysis 5 (ERA5) reanalysis dataset demonstrate that our model yields a significant improvement in both deterministic and ensemble forecasting over the climatology mean, traditional numerical methods, and data-driven models. Ablation studies empirically show the effectiveness of our model designs. Remarkably, our model outperforms skillful numerical ECMWF-S2S and advanced dat
Sea surface temperature (SST) is a fundamental physical parameter characterising the thermal state of sea surface. Due to the intricate thermal interactions between land, sea, and atmosphere, the spatial gradients of SST in coastal waters often appear at finer spatial scales than those in open ocean waters. The Thermal Infrared Sensor (TIRS) onboard Landsat-8, with its 100-meter spatial resolution, offers a unique opportunity to uncover fine-scale coastal SST patterns that would otherwise be overlooked by coarser-resolution thermal sensors. In this study, we first analysed the spatiotemporal patterns of SST in South Australia's temperate coastal waters from 2014 to 2023 by developing an operational approach for SST retrieval from the Landsat-8 TIRS sensor. A buoy was deployed off the coast of Port Lincoln, South Australia, to validate the quality of SST retrievals. Then the daily baseline climatology of SST with 100 m resolution was constructed, which allowed for the detection and analysis of anomalous SST events. Our results suggest the following: (1) the satellite-derived SST data aligned well with the in-situ measured SST values; (2) the semi-enclosed, shallow regions of Upper S
A reliable tropical cyclone (TC) climatology is the key to assessing historical and future changes in TC activities. While global TC records have been systematically maintained since the early 1940s, substantial uncertainties remain for the pre-satellite era during which TC observations relied mostly on scattered aircraft reconnaissance and sporadic ship reports. This study presents a deep learning (DL) approach to reconstruct historical TC activity in the western North Pacific (WNP) basin, with a main focus on the pre-satellite era. Using data feature enrichment tailored for tropical cyclogenesis (TCG), we demonstrate that DL can effectively capture the main characteristics and changes in TCG climatology during the post-satellite era. With additional cross-validations, the reconstruction of TCG climatology is then extended to a pre-satellite period (1940-1960) during which TC base-track datasets are most uncertain. Our DL reconstruction reveals a significant missing of TCG in the current best-track data between September and November during the pre-satellite era. Such a TCG undercount in the best track data occurs mainly around 10-15$^\circ$N in the central WNP, while coastal regi
The Mars ambient space environment evolves with the varying solar activity. Understanding the Martian space environment, particularly the topside ionosphere across different phases of Solar Cycles (SC) 24 \& 25 remains a key research gap in planetary ionospheric science. In this study, we utilized the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission data (150-500 km) from Martian years 32-38 (2015-2024) during solar quiet-time. This study investigated the behavior of topside ionosphere (e-, CO2+, O2+, NO+, OH+, O+, N+ \& C+) across different phases of SC over the northern hemisphere. A significant variation in ionosphere is observed over low-latitude (0-30°N) with higher densities compared to mid-latitude (31-60°N) across SC. Additionally, we found that the Martian northern ionospheric densities were highest during solar maximum phase on both dayside and nightside compared to low active phases. The dayside densities were approximately 1-2 orders higher compared to those on the nightside. The electron and molecular ions densities increased by factors of 1-5 and 1-13, respectively. While O+ ion density was enhanced by nearly 2-2.5 times, along with an upliftment of
This paper investigated the potential of a multivariate Transformer model to forecast the temporal trajectory of the Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) for short (1 month) and long horizon (more than 1 month) periods at the regional level in Europe and North Africa. The input data covers the period from 2002 to 2022 and includes remote sensing and weather data for modelling FAPAR predictions. The model was evaluated using a leave one year out cross-validation and compared with the climatological benchmark. Results show that the transformer model outperforms the benchmark model for one month forecasting horizon, after which the climatological benchmark is better. The RMSE values of the transformer model ranged from 0.02 to 0.04 FAPAR units for the first 2 months of predictions. Overall, the tested Transformer model is a valid method for FAPAR forecasting, especially when combined with weather data and used for short-term predictions.
Because of the rudimentary reporting methods and general lack of documentation, the creation of a severe weather database within the Philippines has been difficult yet relevant target for climatology purposes and historical interest. Previous online severe weather documentation i.e. of tornadoes, waterspouts, and hail events, has also often been few, inconsistent, inactive, or is now completely decommissioned. Several countries or continents support severe weather information through either government-sponsored or independent organizations. For this work, Project SWAP stands as a collaborative exercise, with clear data attribution and open avenues for augmentation, and the creation of a common data model to store the phenomenon's information will assist in maintaining and updating the aforementioned online archive in the Philippines. This paper presents the methods necessary for creating the SWAP database, provide broader climatological analysis of spatio-temporal patterns in severe weather occurrence within the Philippine context, and outline potential use cases for the data. We also highlight the project's current limitations as is to any other existing and far larger database, a