The evolution of biomass-burning aerosol size distributions due to coagulation: dependence on fire and meteorological details and parameterization

Biomass-burning aerosols have a significant effect on global and regional aerosol climate forcings. To model the magnitude of these effects accurately requires knowledge of the size distribution of the emitted and evolving aerosol particles. Current biomass-burning inventories do not include size distributions, and global and regional models generally assume a fixed size distribution from all biomass-burning emissions. However, biomass-burning size distributions evolve in the plume due to coagulation and net organic aerosol (OA) evaporation or formation, and the plume processes occur on spacial scales smaller than global/regional-model grid boxes. The extent of this size-distribution evolution is dependent on a variety of factors relating to the emission source and atmospheric conditions. Therefore, accurately accounting for biomass-burning aerosol size in global models requires an effective aerosol size distribution that accounts for this sub-grid evolution and can be derived from available emission-inventory and meteorological parameters. In this paper, we perform a detailed investigation of the effects of coagulation on the aerosol size distribution in biomass-burning plumes. We compare the effect of coagulation to that of OA evaporation and formation. We develop coagulation-only parameterizations for effective biomass-burning size distributions using the SAM-TOMAS large-eddy simulation plume model. For the most-sophisticated parameterization, we use the Gaussian Emulation Machine for Sensitivity Analysis (GEM-SA) to build a parameterization of the aged size distribution based on the SAM-TOMAS output and seven inputs: emission median dry diameter, emission distribution modal width, mass emissions flux, fire area, mean boundary-layer wind speed, plume mixing depth, and time/distance since emission. This parameterization was tested against an independent set of SAM-TOMAS simulations and yields R2 values of 0.83 and 0.89 for Dpm and modal width, respectively. The size distribution is particularly sensitive to the mass emissions flux, fire area, wind speed, and time, and we provide simplified fits of the aged size distribution to just these input variables. The simplified fits were tested against 11 aged biomass-burning size distributions observed at the Mt. Bachelor Observatory in August 2015. The simple fits captured over half of the variability in observed Dpm and modal width even though the freshly emitted Dpm and modal widths were unknown. These fits may be used in global and regional aerosol models. Finally, we show that coagulation generally leads to greater changes in the particle size distribution than OA evaporation/formation does, using estimates of OA production/loss from the literature

Statistical Emulation for Environmental Sustainability Analysis

The potential effects of climate change on the environment and society are many. In order to effectively quantify the uncertainty associated with these effects, highly complex simulation models are run with detailed representations of ecosystem processes. These models are computationally expensive and can involve computer runs of several days for their outputs. Computationally cheaper models can be obtained from large ensembles of simulations using a statistical emulation. The purpose of this thesis is to construct cheaper computational models (emulators) from simulation outputs of Lund-Potsdam-Jena-managed Land (LPJmL) which is a dynamic global vegetation and crop model. This research work is part of a project called ERMITAGE. The project links together several key component models into a common framework to better understand how the management and interaction of land, water and the earth’s climate system could be improved. The thesis focuses specifically on emulation of major outputs from the LPJmL model; carbon fluxes (NPP, carbon loss due to heterotrophic respiration and fire carbon) and potential crop yields (cereal, rice, maize and oil crops). Future decadal changes in carbon fluxes and crop yields are modelled as linear functions of climate change and other relevant variables. The emulators are constructed using a combination of statistical techniques of stepwise least squares regression, principal component analysis, weighted least squares regression, censored regression and Gaussian process regression. Further modelling involves sensitivity analyses to identify the relative contribution of each input variable to the total output variance. This used the Sobol global sensitivity method. The data cover the period 2001-2100 and comprise climate scenarios of several GCMs and RCPs. Under cross validation the percentage of variance explained ranges from 52-96% for carbon fluxes, 60-88% for the rainfed crops and 62-93% for the irrigated crops, averaged over climate scenarios

Emulating global climate change impacts on crop yields

The potential effects of climate change on the environment and society are many. In order to effectively quantify the uncertainty associated with these effects, highly complex simulation models are run with detailed representations of ecosystem processes. These models are computationally expensive and can involve a computer run of several days. Computationally cheaper models can be obtained from large ensembles of simulations using statistical emulation. The purpose of this paper is to construct a cheaper computational model (emulator) from simulations of the Lund- Potsdam-Jena managed Land (LPJmL), which is a dynamic global vegetation and crop model. This paper focuses on statistical emulation of potential crop yields from LPJmL and an emulator is constructed using a combination of ordinary least squares, principal component analysis and weighted least squares methods. For five climate models, under cross-validation the percentage of variance explained ranges from 60- 88% for the rainfed crops and 62-93% for the irrigated crops. The emulator can be used to predict potential crop yield change under any future climate scenarios and management options

Ensembles of global climate model variants designed for the quantification and constraint of uncertainty in aerosols and their radiative forcing

Tropospheric aerosol radiative forcing has persisted for many years as one of the major causes of uncertainty in global climate model simulations. To sample the range of plausible aerosol and atmospheric states and perform robust statistical analyses of the radiative forcing, it is important to account for the combined effects of many sources of model uncertainty, which is rarely done due to the high computational cost. This paper describes the designs of two ensembles of the HadGEM-UKCA global climate model and provides the first analyses of the uncertainties in aerosol radiative forcing and their causes. The first ensemble was designed to comprehensively sample uncertainty in the aerosol state, while the other samples additional uncertainties in the physical model related to clouds, humidity and radiation, thereby allowing an analysis of uncertainty in the aerosol effective radiative forcing. Each ensemble consists of around 200 simulations of the pre-industrial and present-day atmospheres. The uncertainty in aerosol radiative forcing in our ensembles is comparable to the range of estimates from multi-model intercomparison projects. The mean aerosol effective radiative forcing is –1.45 W m–2 (credible interval –2.07 to –0.81 W m–2), which encompasses but is more negative than the –1.17 W m–2 in the 2013 Atmospheric Chemistry and Climate Model Intercomparison Project and –0.90 W m–2 in the IPCC 5th Assessment Report. The ensembles can be used to reduce aerosol radiative forcing uncertainty by challenging them with multiple measurements as well as to isolate potential causes of multi-model differences

Emulation of a chemical transport model to assess air quality under future emission scenarios for the southwest of Western Australia

Simulation outputs from chemical transport models (CTMs) are essential to plan effective air quality policies. A key strength of these models is their ability to separate out source-specific components which facilitate the simulation of the potential impact of policy on future air quality. However, configuring and running these models is complex and computationally intensive, making the evaluation of multiple scenarios less accessible to many researchers and policy experts. The aim of this work is to present how Gaussian process emulation can provide a top-down approach to interrogating and interpreting the outputs from CTMs at minimal computational cost. A case study is presented (based on fine particle sources in the southwest of Western Australia) to illustrate how an emulator can be constructed to simultaneously evaluate changes in emissions from on-road transport and electricity sectors. This study demonstrates how emulation provides a flexible way of exploring local impacts of electric vehicles and wider regional effects of emissions from electricity generation. The potential for emulators to be applied to other settings involving air quality research is discussed

Report from the Tri-Agency Cosmological Simulation Task Force

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.Comment: 36 pages, 3 figures. Delivered to NASA, NSF, and DOE in Dec 201

Statistical Emulation of Winter Ambient Fine Particulate Matter Concentrations From Emission Changes in China

Air pollution exposure remains a leading public health problem in China. The use of chemical transport models to quantify the impacts of various emission changes on air quality is limited by their large computational demands. Machine learning models can emulate chemical transport models to provide computationally efficient predictions of outputs based on statistical associations with inputs. We developed novel emulators relating emission changes in five key anthropogenic sectors (residential, industry, land transport, agriculture, and power generation) to winter ambient fine particulate matter (PM2.5) concentrations across China. The emulators were optimized based on Gaussian process regressors with Matern kernels. The emulators predicted 99.9% of the variance in PM2.5 concentrations for a given input configuration of emission changes. PM2.5 concentrations are primarily sensitive to residential (51%–94% of first‐order sensitivity index), industrial (7%–31%), and agricultural emissions (0%–24%). Sensitivities of PM2.5 concentrations to land transport and power generation emissions are all under 5%, except in South West China where land transport emissions contributed 13%. The largest reduction in winter PM2.5 exposure for changes in the five emission sectors is by 68%–81%, down to 15.3–25.9 μg m−3, remaining above the World Health Organization annual guideline of 10 μg m−3. The greatest reductions in PM2.5 exposure are driven by reducing residential and industrial emissions, emphasizing the importance of emission reductions in these key sectors. We show that the annual National Air Quality Target of 35 μg m−3 is unlikely to be achieved during winter without strong emission reductions from the residential and industrial sectors

AFIT School of Engineering Contributions to Air Force Research and Technology Calendar Year 1973

This report contains abstracts of Master of Science Theses, Doctoral dissertations, and selected faculty publications completed during the 1973 calendar year at the School of Engineering, Air Force Institute of Technology, at Wright-Patterson Air Force Base, Ohio