Dependence of fast changes in global and local precipitation on the geographical location of absorbing aerosol

Anthropogenic aerosol interacts strongly with incoming solar radiation, perturbing Earth’s energy budget and precipitation on both local and global scales. Understanding these changes in precipitation has proven particularly difficult for the case of absorbing aerosol, which absorbs a significant amount of incoming solar radiation and hence acts as a source of localized diabatic heating to the atmosphere. In this work, we use an ensemble of atmosphere-only climate model simulations forced by identical absorbing aerosol perturbations in different geographical locations across the globe to develop a basic physical understanding of how this localized heating impacts the atmosphere and how these changes impact on precipitation both globally and locally. In agreement with previous studies we find that absorbing aerosol causes a decrease in global-mean precipitation, but we also show that even for identical aerosol optical depth perturbations, the global-mean precipitation change varies by over an order of magnitude depending on the location of the aerosol burden. Our experiments also demonstrate that the local precipitation response to absorbing aerosol is opposite in sign between the tropics and the extratropics, as found by previous work. We then show that this contrasting response can be understood in terms of different mechanisms by which the large-scale circulation responds to heating in the extratropics and in the tropics. We provide a simple theory to explain variations in the local precipitation response to absorbing aerosol in the tropics. Our work highlights that the spatial pattern of absorbing aerosol and its interactions with circulation are a key determinant of its overall climate impact and must be taken into account when developing our understanding of aerosol–climate interactions

Exploring Randomly Wired Neural Networks for Climate Model Emulation

Exploring the climate impacts of various anthropogenic emissions scenarios is key to making informed decisions for climate change mitigation and adaptation. State-of-the-art Earth system models can provide detailed insight into these impacts, but have a large associated computational cost on a per-scenario basis. This large computational burden has driven recent interest in developing cheap machine learning models for the task of climate model emulation. In this manuscript, we explore the efficacy of randomly wired neural networks for this task. We describe how they can be constructed and compare them to their standard feedforward counterparts using the ClimateBench dataset. Specifically, we replace the serially connected dense layers in multilayer perceptrons, convolutional neural networks, and convolutional long short-term memory networks with randomly wired dense layers and assess the impact on model performance for models with 1 million and 10 million parameters. We find average performance improvements of 4.2% across model complexities and prediction tasks, with substantial performance improvements of up to 16.4% in some cases. Furthermore, we find no significant difference in prediction speed between networks with standard feedforward dense layers and those with randomly wired layers. These findings indicate that randomly wired neural networks may be suitable direct replacements for traditional dense layers in many standard models

Detecting anthropogenic cloud perturbations with deep learning

One of the most pressing questions in climate science is that of the effect of anthropogenic aerosol on the Earth's energy balance. Aerosols provide the `seeds' on which cloud droplets form, and changes in the amount of aerosol available to a cloud can change its brightness and other physical properties such as optical thickness and spatial extent. Clouds play a critical role in moderating global temperatures and small perturbations can lead to significant amounts of cooling or warming. Uncertainty in this effect is so large it is not currently known if it is negligible, or provides a large enough cooling to largely negate present-day warming by CO2. This work uses deep convolutional neural networks to look for two particular perturbations in clouds due to anthropogenic aerosol and assess their properties and prevalence, providing valuable insights into their climatic effects.Comment: Awarded Best Paper and Spotlight Oral at Climate Change: How Can AI Help? (Workshop) at International Conference on Machine Learning (ICML), Long Beach, California, 201

Cumulo: A Dataset for Learning Cloud Classes

One of the greatest sources of uncertainty in future climate projections comes from limitations in modelling clouds and in understanding how different cloud types interact with the climate system. A key first step in reducing this uncertainty is to accurately classify cloud types at high spatial and temporal resolution. In this paper, we introduce Cumulo, a benchmark dataset for training and evaluating global cloud classification models. It consists of one year of 1km resolution MODIS hyperspectral imagery merged with pixel-width 'tracks' of CloudSat cloud labels. Bringing these complementary datasets together is a crucial first step, enabling the Machine-Learning community to develop innovative new techniques which could greatly benefit the Climate community. To showcase Cumulo, we provide baseline performance analysis using an invertible flow generative model (IResNet), which further allows us to discover new sub-classes for a given cloud class by exploring the latent space. To compare methods, we introduce a set of evaluation criteria, to identify models that are not only accurate, but also physically-realistic. CUMULO can be download from https://www.dropbox.com/sh/6gca7f0mb3b0ikz/AADq2lk4u7k961Qa31FwIDEpa?dl=0

Pollution tracker: finding industrial sources of aerosol emission in satellite imagery

The effects of anthropogenic aerosol, solid or liquid particles suspended in the air, are the biggest contributor to uncertainty in current climate perturbations. Heavy industry sites, such as coal power plants and steel manufacturers, emit large amounts of aerosol in a small area. This makes them ideal places to study aerosol interactions with radiation and clouds. However, existing data sets of heavy industry locations are either not public, or suffer from reporting gaps. Here, we develop a deep learning algorithm to detect unreported industry sites in high-resolution satellite data. For the pipeline to be viable at global scale, we employ a two-step approach. The first step uses 10 m resolution data, which is scanned for potential industry sites, before using 1.2 m resolution images to confirm or reject detections. On held out test data, the models perform well, with the lower resolution one reaching up to 94% accuracy. Deployed to a large test region, the first stage model yields many false positive detections. The second stage, higher resolution model shows promising results at filtering these out, while keeping the true positives. In the deployment area, we find five new heavy industry sites which were not in the training data. This demonstrates that the approach can be used to complement data sets of heavy industry sites

Community Intercomparison Suite (CIS) v1.4.0: a tool for intercomparing models and observations

The Community Intercomparison Suite (CIS) is an easy-to-use command-line tool which has been developed to allow the straightforward intercomparison of remote sensing, in situ and model data. While there are a number of tools available for working with climate model data, the large diversity of sources (and formats) of remote sensing and in situ measurements necessitated a novel software solution. Developed by a professional software company, CIS supports a large number of gridded and ungridded data sources "out-of-the-box", including climate model output in NetCDF or the UK Met Office pp file format, CloudSat, CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization), MODIS (MODerate resolution Imaging Spectroradiometer), Cloud and Aerosol CCI (Climate Change Initiative) level 2 satellite data and a number of in situ aircraft and ground station data sets. The open-source architecture also supports user-defined plugins to allow many other sources to be easily added. Many of the key operations required when comparing heterogenous data sets are provided by CIS, including subsetting, aggregating, collocating and plotting the data. Output data are written to CF-compliant NetCDF files to ensure interoperability with other tools and systems. The latest documentation, including a user manual and installation instructions, can be found on our website (http://cistools.net). Here, we describe the need which this tool fulfils, followed by descriptions of its main functionality (as at version 1.4.0) and plugin architecture which make it unique in the field

Tobac 1.2: Towards a flexible framework for tracking and analysis of clouds in diverse datasets

We introduce tobac (Tracking and Object-Based Analysis of Clouds), a newly developed framework for tracking and analysing individual clouds in different types of datasets, such as cloud-resolving model simulations and geostationary satellite retrievals. The software has been designed to be used flexibly with any two-or three-dimensional timevarying input. The application of high-level data formats, such as Iris cubes or xarray arrays, for input and output allows for convenient use of metadata in the tracking analysis and visualisation. Comprehensive analysis routines are provided to derive properties like cloud lifetimes or statistics of cloud properties along with tools to visualise the results in a convenient way. The application of tobac is presented in two examples. We first track and analyse scattered deep convective cells based on maximum vertical velocity and the threedimensional condensate mixing ratio field in cloud-resolving model simulations. We also investigate the performance of the tracking algorithm for different choices of time resolution of the model output. In the second application, we show how the framework can be used to effectively combine information from two different types of datasets by simultaneously tracking convective clouds in model simulations and in geostationary satellite images based on outgoing longwave radiation. The tobac framework provides a flexible new way to include the evolution of the characteristics of individual clouds in a range of important analyses like model intercomparison studies or model assessment based on observational data. © 2019 Author(s)

Increased water vapour lifetime due to global warming

Water vapour in the atmosphere is the source of a major climate feedback mechanism and potential increases in the availability of water vapour could have important consequences for mean and extreme precipitation. Future precipitation changes further depend on how the hydrological cycle responds to drivers of climate change, such as greenhouse gases and aerosols. Currently, neither the total anthropogenic influence on the hydrological cycle, nor those from individual drivers, are constrained sufficiently to make solid projections. We investigate how integrated water vapour (IWV) responds to different drivers of climate change. Results from 11 global climate models have been used, based on simulations where CO2, methane, solar irradiance, black carbon (BC), and sulphate have been perturbed separately. While the global-mean IWV is usually assumed to increase by ~7% per degree K surface temperature change, we find that the feedback response of IWV differs somewhat between drivers. Fast responses, which include the initial radiative effect and rapid adjustments to an external forcing, amplify these differences. The resulting net changes in IWV range from 6.4±0.9%/K for sulphate to 9.8±2%/K for BC. We further calculate the relationship between global changes in IWV and precipitation, which can be characterized by quantifying changes in atmospheric water vapour lifetime. Global climate models simulate a substantial increase in the lifetime, from 8.2±0.5 to 9.9±0.7 days between 1986-2005 and 2081-2100 under a high emission scenario, and we discuss to what extent the water vapour lifetime provides additional information compared to analysis of IWV and precipitation separately. We conclude that water vapour lifetime changes are an important indicator of changes in precipitation patterns and that BC is particularly efficient in prolonging the distance between evaporation and precipitation

How well are aerosol–cloud interactions represented in climate models? – Part 1: Understanding the sulfate aerosol production from the 2014–15 Holuhraun eruption

For over 6 months, the 2014–2015 effusive eruption at Holuhraun, Iceland, injected considerable amounts of sulfur dioxide (SO2) into the lower troposphere with a daily rate of up to one-third of the global emission rate, causing extensive air pollution across Europe. The large injection of SO2, which oxidises to form sulfate aerosol (SO42-), provides a natural experiment offering an ideal opportunity to scrutinise state-of-the-art general circulation models' (GCMs) representation of aerosol–cloud interactions (ACIs). Here we present Part 1 of a two-part model inter-comparison using the Holuhraun eruption as a framework to analyse ACIs. We use SO2 retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) instrument and ground-based measurements of SO2 and SO42- mass concentrations across Europe, in conjunction with a trajectory analysis using the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model, to assess the spatial and chemical evolution of the volcanic plume as simulated by five GCMs and a chemical transport model (CTM). IASI retrievals of plume altitude and SO2 column load reveal that the volcanic perturbation is largely contained within the lower troposphere. Compared to the satellite observations, the models capture the spatial evolution and vertical variability of the plume reasonably well, although the models often overestimate the plume altitude. HYSPLIT trajectories are used to attribute to Holuhraun emissions 111 instances of elevated sulfurous surface mass concentrations recorded at European Monitoring and Evaluation Programme (EMEP) stations during September and October 2014. Comparisons with the simulated concentrations show that the modelled ratio of SO2 to SO42- during these pollution episodes is often underestimated and overestimated for the young and mature plume, respectively. Models with finer vertical resolutions near the surface are found to better capture these elevated sulfurous ground-level concentrations. Using an exponential function to describe the decay of observed surface mass concentration ratios of SO2 to SO42- with plume age, the in-plume oxidation rate constant is estimated as 0.032 ± 0.002 h−1 (1.30 ± 0.08 d e-folding time), with a near-vent ratio of 25 ± 5 (µg m−3 of SO2 / µg m−3 of SO42-). The majority of the corresponding derived modelled oxidation rate constants are lower than the observed estimate. This suggests that the representation of the oxidation pathway/s in the simulated plumes is too slow. Overall, despite their coarse spatial resolutions, the six models show reasonable skill in capturing the spatial and chemical evolution of the Holuhraun plume. This capable representation of the underlying aerosol perturbation is essential to enable the investigation of the eruption's impact on ACIs in the second part of this study.</p

ClimateBench v1.0: A benchmark for data-driven climate projections

Many different emission pathways exist that are compatible with the Paris climate agreement, and many more are possible that miss that target. While some of the most complex Earth System Models have simulated a small selection of Shared Socioeconomic Pathways, it is impractical to use these expensive models to fully explore the space of possibilities. Such explorations therefore mostly rely on one-dimensional impulse response models, or simple pattern scaling approaches to approximate the physical climate response to a given scenario. Here we present ClimateBench - a benchmarking framework based on a suite of CMIP, AerChemMIP and DAMIP simulations performed by a full complexity Earth System Model, and a set of baseline machine learning models that emulate its response to a variety of forcers. These emulators can predict annual mean global distributions of temperature, diurnal temperature range and precipitation (including extreme precipitation) given a wide range of emissions and concentrations of carbon dioxide, methane and aerosols, allowing them to efficiently probe previously unexplored scenarios. We discuss the accuracy and interpretability of these emulators and consider their robustness to physical constraints such as total energy conservation. Future opportunities incorporating such physical constraints directly in the machine learning models and using the emulators for detection and attribution studies are also discussed. This opens a wide range of opportunities to improve prediction, consistency and mathematical tractability. We hope that by laying out the principles of climate model emulation with clear examples and metrics we encourage others to tackle this important and demanding challenge