Random forest meteorological normalisation models for Swiss PM10 trend analysis

Meteorological normalisation is a technique which accounts for changes in meteorology over time in an air quality time series. Controlling for such changes helps support robust trend analysis because there is more certainty that the observed trends are due to changes in emissions or chemistry, not changes in meteorology. Predictive random forest models (RF; a decision tree machine learning technique) were grown for 31 air quality monitoring sites in Switzerland using surface meteorological, synoptic scale, boundary layer height, and time variables to explain daily PM10 concentrations. The RF models were used to calculate meteorologically normalised trends which were formally tested and evaluated using the Theil–Sen estimator. Between 1997 and 2016, significantly decreasing normalised PM10 trends ranged between −0.09 and −1.16 μg m−3 year−1 with urban traffic sites experiencing the greatest mean decrease in PM10 concentrations at −0.77 μg m−3 year−1. Similar magnitudes have been reported for normalised PM10 trends for earlier time periods in Switzerland which indicates PM10 concentrations are continuing to decrease at similar rates as in the past. The ability for RF models to be interpreted was leveraged using partial dependence plots to explain the observed trends and relevant physical and chemical processes influencing PM10 concentrations. Notably, two regimes were suggested by the models which cause elevated PM10 concentrations in Switzerland: one related to poor dispersion conditions and a second resulting from high rates of secondary PM generation in deep, photochemically active boundary layers. The RF meteorological normalisation process was found to be robust, user friendly and simple to implement, and readily interpretable which suggests the technique could be useful in many air quality exploratory data analysis situations

The radiative effect of ion-induced inorganic nucleation in the free troposphere

To determine the effect of cosmic rays on the Earth's climate via ion-induced nucleation, a parametrisation of inorganic nucleation was formulated based on experiments at the CERN CLOUD experiment. The parametrisation was implemented in the GLOMAP aerosol microphysics model and used to estimate the radiative effect of the change in ionisation experienced over an 11-year solar cycle

Aerosol climate feedback due to decadal increases in Southern Hemisphere wind speeds

Observations indicate that the westerly jet in the Southern Hemisphere troposphere is accelerating. Using a global aerosol model we estimate that the increase in wind speed of 0.45 + /- 0.2 m s(-1) decade(-1) at 50-65 degrees S since the early 1980s caused a higher sea spray flux, resulting in an increase of cloud condensation nucleus concentrations of more than 85% in some regions, and of 22% on average between 50 and 65 degrees S. These fractional increases are similar in magnitude to the decreases over many northern hemisphere land areas due to changes in air pollution over the same period. The change in cloud drop concentrations causes an increase in cloud reflectivity and a summertime radiative forcing between at 50 and 65 degrees S comparable in magnitude but acting against that from greenhouse gas forcing over the same time period, and thus represents a substantial negative climate feedback. However, recovery of Antarctic ozone depletion in the next two decades will likely cause a fall in wind speeds, a decrease in cloud drop concentration and a correspondingly weaker cloud feedback

Statistical constraints on climate model parameters using a scalable cloud-based inference framework

Atmospheric aerosols influence the Earth’s climate, primarily by affecting cloud formation and scattering visible radiation. However, aerosol-related physical processes in climate simulations are highly uncertain. Constraining these processes could help improve model-based climate predictions. We propose a scalable statistical framework for constraining the parameters of expensive climate models by comparing model outputs with observations. Using the C3.AI Suite, a cloud computing platform, we use a perturbed parameter ensemble of the UKESM1 climate model to efficiently train a surrogate model. A method for estimating a data-driven model discrepancy term is described. The strict bounds method is applied to quantify parametric uncertainty in a principled way. We demonstrate the scalability of this framework with 2 weeks’ worth of simulated aerosol optical depth data over the South Atlantic and Central African region, written from the model every 3 hr and matched in time to twice-daily MODIS satellite observations. When constraining the model using real satellite observations, we establish constraints on combinations of two model parameters using much higher time-resolution outputs from the climate model than previous studies. This result suggests that within the limits imposed by an imperfect climate model, potentially very powerful constraints may be achieved when our framework is scaled to the analysis of more observations and for longer time periods

Aerosol and physical atmosphere model parameters are both important sources of uncertainty in aerosol ERF

Changes in aerosols cause a change in net top-of-the-atmosphere (ToA) short-wave and long-wave radiative fluxes, rapid adjustments in clouds, water vapour and temperature, and cause an effective radiative forcing (ERF) of the planetary energy budget. The diverse sources of model uncertainty and the computational cost of running climate models make it difficult to isolate the main causes of aerosol ERF uncertainty and to understand how observations can be used to constrain it. We explore the aerosol ERF uncertainty by using fast model emulators to generate a very large set of aerosol-climate model variants that span the model uncertainty due to twenty-seven parameters related to atmospheric and aerosol processes. Sensitivity analyses shows that the uncertainty in the ToA flux is dominated (around 80 %) by uncertainties in the physical atmosphere model, particularly parameters that affect cloud reflectivity. However, uncertainty in the change in ToA flux caused by aerosol emissions over the industrial period (the aerosol ERF) is controlled by a combination of uncertainties in aerosol (around 60 %) and physical atmosphere (around 40 %) parameters. Four of the atmospheric and aerosol parameters that cause uncertainty in short-wave ToA flux (mostly parameters that directly scale cloud reflectivity, cloud water content or cloud droplet concentrations) also account for around 60% of the aerosol ERF uncertainty. The common causes of uncertainty mean that constraining the modelled planetary brightness to tightly match satellite observations changes the lower 95 % credible aerosol ERF value from −2.65 Wm−2 to −2.37 Wm−2. This suggests the strongest forcings (below around −2.4 Wm−2) are inconsistent with observations. These results show that, regardless of the fact that the ToA flux is an order of magnitude larger than the aerosol ERF, the observed flux can constrain the uncertainty in ERF because their values are connected by constrainable process parameters. The key to reducing the aerosol ERF uncertainty further will be to identify observations that can additionally constrain individual parameter ranges and/or combined parameter effects, which can be achieved through sensitivity analysis of perturbed parameter ensembles

A Simple Model of Global Aerosol Indirect Effects

Most estimates of the global mean indirect effect of anthropogenic aerosol on the Earth's energy balance are from simulations by global models of the aerosol lifecycle coupled with global models of clouds and the hydrologic cycle. Extremely simple models have been developed for integrated assessment models, but lack the flexibility to distinguish between primary and secondary sources of aerosol. Here a simple but more physically based model expresses the aerosol indirect effect (AIE) using analytic representations of cloud and aerosol distributions and processes. Although the simple model is able to produce estimates of AIEs that are comparable to those from some global aerosol models using the same global mean aerosol properties, the estimates by the simple model are sensitive to preindustrial cloud condensation nuclei concentration, preindustrial accumulation mode radius, width of the accumulation mode, size of primary particles, cloud thickness, primary and secondary anthropogenic emissions, the fraction of the secondary anthropogenic emissions that accumulates on the coarse mode, the fraction of the secondary mass that forms new particles, and the sensitivity of liquid water path to droplet number concentration. Estimates of present-day AIEs as low as 5 W/sq m and as high as 0.3 W/sq m are obtained for plausible sets of parameter values. Estimates are surprisingly linear in emissions. The estimates depend on parameter values in ways that are consistent with results from detailed global aerosol-climate simulation models, which adds to understanding of the dependence on AIE uncertainty on uncertainty in parameter values

Explaining global surface aerosol number concentrations in terms of primary emissions and particle formation

We use observations of total particle number concentration at 36 worldwide sites and a global aerosol model to quantify the primary and secondary sources of particle number. We show that emissions of primary particles can reasonably reproduce the spatial pattern of observed condensation nuclei (CN) (R2=0.51) but fail to explain the observed seasonal cycle at many sites (R2=0.1). The modeled CN concentration in the free troposphere is biased low (normalised mean bias, NMB=−88%) unless a secondary source of particles is included, for example from binary homogeneous nucleation of sulfuric acid and water (NMB=−25%). Simulated CN concentrations in the continental boundary layer (BL) are also biased low (NMB=−74%) unless the number emission of anthropogenic primary particles is increased or an empirical BL particle formation mechanism based on sulfuric acid is used. We find that the seasonal CN cycle observed at continental BL sites is better simulated by including a BL particle formation mechanism (R2=0.3) than by increasing the number emission from primary anthropogenic sources (R2=0.18). Using sensitivity tests we derive optimum rate coefficients for this nucleation mechanism, which agree with values derived from detailed case studies at individual sites

New Approaches to Quantifying Aerosol Influence on the Cloud Radiative Effect

The topic of cloud radiative forcing associated with the atmospheric aerosol has been the focus of intense scrutiny for decades. The enormity of the problem is reflected in the need to understand aspects such as aerosol composition, optical properties, cloud condensation, and ice nucleation potential, along with the global distribution of these properties, controlled by emissions, transport, transformation, and sinks. Equally daunting is that clouds themselves are complex, turbulent, microphysical entities and, by their very nature, ephemeral and hard to predict. Atmospheric general circulation models represent aerosol−cloud interactions at ever-increasing levels of detail, but these models lack the resolution to represent clouds and aerosol−cloud interactions adequately. There is a dearth of observational constraints on aerosol−cloud interactions. We develop a conceptual approach to systematically constrain the aerosol−cloud radiative effect in shallow clouds through a combination of routine process modeling and satellite and surface-based shortwave radiation measurements. We heed the call to merge Darwinian and Newtonian strategies by balancing microphysical detail with scaling and emergent properties of the aerosol−cloud radiation system

Results from the CERN pilot CLOUD experiment

During a 4-week run in October–November 2006, a pilot experiment was performed at the CERN Proton Synchrotron in preparation for the Cosmics Leaving OUtdoor Droplets (CLOUD) experiment, whose aim is to study the possible influence of cosmic rays on clouds. The purpose of the pilot experiment was firstly to carry out exploratory measurements of the effect of ionising particle radiation on aerosol formation from trace H2SO4 vapour and secondly to provide technical input for the CLOUD design. A total of 44 nucleation bursts were produced and recorded, with formation rates of particles above the 3 nm detection threshold of between 0.1 and 100 cm -3 s -1, and growth rates between 2 and 37 nm h -1. The corresponding H2O concentrations were typically around 106 cm -3 or less. The experimentally-measured formation rates and htwosofour concentrations are comparable to those found in the atmosphere, supporting the idea that sulphuric acid is involved in the nucleation of atmospheric aerosols. However, sulphuric acid alone is not able to explain the observed rapid growth rates, which suggests the presence of additional trace vapours in the aerosol chamber, whose identity is unknown. By analysing the charged fraction, a few of the aerosol bursts appear to have a contribution from ion-induced nucleation and ion-ion recombination to form neutral clusters. Some indications were also found for the accelerator beam timing and intensity to influence the aerosol particle formation rate at the highest experimental SO2 concentrations of 6 ppb, although none was found at lower concentrations. Overall, the exploratory measurements provide suggestive evidence for ion-induced nucleation or ion-ion recombination as sources of aerosol particles. However in order to quantify the conditions under which ion processes become significant, improvements are needed in controlling the experimental variables and in the reproducibility of the experiments. Finally, concerning technical aspects, the most important lessons for the CLOUD design include the stringent requirement of internal cleanliness of the aerosol chamber, as well as maintenance of extremely stable temperatures (variations below 0.1 °C