skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: ISORROPIA‐MCX: Enabling Sensitivity Analysis With Multicomplex Variables in the Aerosol Thermodynamic Model, ISORROPIA
Abstract Sensitivity analysis with atmospheric chemical transport models may be used to quantify influences of specific emissions on pollutant concentrations. This information facilitates efficient environmental decision‐making regarding emissions control strategies for pollutants that affect human health and public welfare. The multicomplex step method (MCX) is a sensitivity analysis approach that enables calculation of first‐ and higher‐order sensitivities of a nonlinear algorithm with analytical accuracy. Compared to the well‐known finite difference method, the MCX method is also straight‐forward to compute yet does not suffer from precision errors due to subtracting numbers with common leading digits and eliminates the requirement of tuning the step size. The aerosol inorganic equilibrium thermodynamic model, ISORROPIA, which treats ammonium, chloride, nitrate, sodium, sulfate, calcium, potassium, and magnesium, was augmented to leverage the multicomplex step method (ISORROPIA‐MCX) to analyze the influence that the total amount of a pollutant has on concentrations partitioned into different phases. This enables simultaneous calculation of the first‐order, second‐order, and cross‐sensitivity terms in the Taylor Series expansion when evaluating the impact of changes in input parameters on an output variable, increasing the accuracy of the estimated effect when the functions are nonlinear. ISORROPIA encodes highly nonlinear processes which showcases the computational advantages of the multicomplex step method as well as the limitations of the approach for fractured solution surfaces. With ISORROPIA‐MCX, the influence of total concentrations of aerosol precursors on aerosol acidity are evaluated with cross‐sensitivity terms for the first time.  more » « less
Award ID(s):
1944669
PAR ID:
10426714
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Earth and Space Science
Volume:
10
Issue:
6
ISSN:
2333-5084
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Yilong Wang (Ed.)
    Abstract. Sensitivity analysis in chemical transport models quantifies the response of output variables to changes in input parameters. This information is valuable for researchers engaged in data assimilation and model development. Additionally, environmental decision-makers depend upon these expected responses of concentrations to emissions when designing and justifying air pollution control strategies. Existing sensitivity analysis methods include the finite-difference method, the direct decoupled method (DDM), the complex variable method, and the adjoint method. These methods are either prone to significant numerical errors when applied to nonlinear models with complex components (e.g. finite difference and complex step methods) or difficult to maintain when the original model is updated (e.g. direct decoupled and adjoint methods). Here, we present the implementation of the hyperdual-step method in the Community Multiscale Air Quality Model (CMAQ) version 5.3.2 as CMAQ-hyd. CMAQ-hyd can be applied to compute numerically exact first- and second-order sensitivities of species concentrations with respect to emissions or concentrations. Compared to CMAQ-DDM and CMAQ-adjoint, CMAQ-hyd is more straightforward to update and maintain, while it remains free of subtractive cancellation and truncation errors, just as those augmented models do. To evaluate the accuracy of the implementation, the sensitivities computed by CMAQ-hyd are compared with those calculated with other traditional methods or a hybrid of the traditional and advanced methods. We demonstrate the capability of CMAQ-hyd with the newly implemented gas-phase chemistry and biogenic aerosol formation mechanism in CMAQ. We also explore the cross-sensitivity of monoterpene nitrate aerosol formation to its anthropogenic and biogenic precursors to show the additional sensitivity information computed by CMAQ-hyd. Compared with the traditional finite difference method, CMAQ-hyd consumes fewer computational resources when the same sensitivity coefficients are calculated. This novel method implemented in CMAQ is also computationally competitive with other existing methods and could be further optimized to reduce memory and computational time overheads. 
    more » « less
  2. null (Ed.)
    Particle acidity (aerosol pH) is an important driver of atmospheric chemical processes and the resulting effects on human and environmental health. Understanding the factors that control aerosol pH is critical when enacting control strategies targeting specific outcomes. This study characterizes aerosol pH at a land-water transition site near Baltimore, MD during summer 2018 as part of the second Ozone Water-Land Environmental Transition Study (OWLETS-2) field campaign. Inorganic fine mode aerosol composition, gas-phase NH3 measurements, and all relevant meteorological parameters were used to characterize the effects of temperature, aerosol liquid water (ALW), and composition on predictions of aerosol pH. Temperature, the factor linked to the control of NH3 partitioning, was found to have the most significant effect on aerosol pH during OWLETS-2. Overall, pH varied with temperature at a rate of −0.047 K−1 across all observations, though the sensitivity was −0.085 K−1 for temperatures > 293 K. ALW had a minor effect on pH, except at the lowest ALW levels (< 1 µg m−3) which caused a significant increase in aerosol acidity (decrease in pH). Aerosol pH was generally insensitive to composition (SO42− , SO42−:NH4+ , Tot-NH3 = NH3 + NH4+), consistent with recent studies in other locations. In a companion paper, the sources of episodic NH3 events (95th percentile concentrations, NH3 > 7.96 µg m−3) during the study are analyzed; aerosol pH was higher by only ~0.1–0.2 pH units during these events compared to the study mean. A case study was analyzed to characterize the response of aerosol pH to nonvolatile cations (NVCs) during a period strongly influenced by primary Chesapeake Bay emissions. Depending on the method used, aerosol pH was estimated to be either weakly (~0.1 pH unit change based on NH3 partitioning calculation) or strongly (~1.4 pH unit change based on ISORROPIA thermodynamic model predictions) affected by NVCs. The case study suggests a strong pH gradient with size during the event and underscores the need to evaluate assumptions of aerosol mixing state applied to pH calculations. Unique features of this study, including the urban land-water transition site and the strong influence of NH3 emissions from both agricultural and industrial sources, add to the understanding of aerosol pH and its controlling factors in diverse environments. 
    more » « less
  3. Abstract. Particle acidity (aerosol pH) is an important driver of atmospheric chemical processes and the resulting effects on human and environmentalhealth. Understanding the factors that control aerosol pH is critical when enacting control strategies targeting specific outcomes. This studycharacterizes aerosol pH at a land–water transition site near Baltimore, MD, during summer 2018 as part of the second Ozone Water-Land EnvironmentalTransition Study (OWLETS-2) field campaign. Inorganic fine-mode aerosol composition, gas-phase NH3 measurements, and all relevantmeteorological parameters were used to characterize the effects of temperature, aerosol liquid water (ALW), and composition on predictions ofaerosol pH. Temperature, the factor linked to the control of NH3 partitioning, was found to have the most significant effect on aerosol pHduring OWLETS-2. Overall, pH varied with temperature at a rate of −0.047 K−1 across all observations, though the sensitivity was−0.085 K−1 for temperatures > 293 K. ALW had a minor effect on pH, except at the lowest ALW levels(< 1 µg m−3), which caused a significant increase in aerosol acidity (decrease in pH). Aerosol pH was generally insensitive tocomposition (SO42-, SO42-:NH4+, total NH3 (Tot-NH3) = NH3 + NH4+), consistentwith recent studies in other locations. In a companion paper, the sources of episodic NH3 events (95th percentile concentrations,NH3 > 7.96 µg m−3) during the study are analyzed; aerosol pH was higher by only ∼ 0.1–0.2 pH unitsduring these events compared to the study mean. A case study was analyzed to characterize the response of aerosol pH to nonvolatile cations (NVCs)during a period strongly influenced by primary Chesapeake Bay emissions. Depending on the method used, aerosol pH was estimated to be either weakly(∼ 0.1 pH unit change based on NH3 partitioning calculation) or strongly (∼ 1.4 pH unit change based onISORROPIA thermodynamic model predictions) affected by NVCs. The case study suggests a strong pH gradient with size during the event and underscores the need to evaluate assumptions of aerosol mixing state applied to pH calculations. Unique features of this study, including the urban land–water transition site and the strong influence of NH3 emissions from both agricultural and industrial sources, add to the understanding of aerosol pH and its controlling factors in diverseenvironments. 
    more » « less
  4. Abstract During summer, ammonia emissions in Southeast Asia influence air pollution and cloud formation. Convective transport by the South Asian monsoon carries these pollutant air masses into the upper troposphere and lower stratosphere (UTLS), where they accumulate under anticyclonic flow conditions. This air mass accumulation is thought to contribute to particle formation and the development of the Asian Tropopause Aerosol Layer (ATAL). Despite the known influence of ammonia and particulate ammonium on air pollution, a comprehensive understanding of the ATAL is lacking. In this modelling study, the influence of ammonia on particle formation is assessed with emphasis on the ATAL. We use the EMAC chemistry-climate model, incorporating new particle formation parameterisations derived from experiments at the CERN CLOUD chamber. Our diurnal cycle analysis confirms that new particle formation mainly occurs during daylight, with a 10-fold enhancement in rate. This increase is prominent in the South Asian monsoon UTLS, where deep convection introduces high ammonia levels from the boundary layer, compared to a baseline scenario without ammonia. Our model simulations reveal that this ammonia-driven particle formation and growth contributes to an increase of up to 80% in cloud condensation nuclei (CCN) concentrations at cloud-forming heights in the South Asian monsoon region. We find that ammonia profoundly influences the aerosol mass and composition in the ATAL through particle growth, as indicated by an order of magnitude increase in nitrate levels linked to ammonia emissions. However, the effect of ammonia-driven new particle formation on aerosol mass in the ATAL is relatively small. Ammonia emissions enhance the regional aerosol optical depth (AOD) for shortwave solar radiation by up to 70%. We conclude that ammonia has a pronounced effect on the ATAL development, composition, the regional AOD, and CCN concentrations. 
    more » « less
  5. The formation of inorganic fine particulate matter (i.e., iPM2.5) is controlled by the thermodynamic equilibrium partitioning of NH3-NH4+. To develop effective control strategies of PM2.5, we aim to understand the impacts of changes in different precursor gases on iPM2.5 concentrations and partitioning of NH3-NH4+. To understand partitioning of NH3-NH4+ in the southeastern U.S., responses of iPM2.5 to precursor gases in four seasons were investigated using field measurements of iPM2.5, precursor gases, and meteorological conditions. The ISORROPIA II model was used to examine the effects of changes in total ammonia (gas + aerosol), total sulfuric acid (aerosol), and total nitric acid (gas + aerosol) on iPM2.5 concentrations and partitioning of NH3-NH4+. The results indicate that reduction in total H2SO4 is more effective than reduction in total HNO3 and total NH3 to reduce iPM2.5 especially under NH3-rich condition. The reduction in total H2SO4 may change partitioning of NH3-NH4+ towards gas-phase and may also lead to an increase in NO3− under NH3-rich conditions, which does not necessarily lead to full neutralization of acidic gases (pH < 7). Thus, future reduction in iPM2.5 may necessitate the coordinated reduction in both H2SO4 and HNO3 in the southeastern U.S. It is also found that the response of iPM2.5 to the change in total H2SO4 is more sensitive in summer than winter due to the dominance of SO42− salts in iPM2.5 and the high temperature in summer. The NH3 emissions from Animal Feeding Operations (AFOs) at an agricultural rural site (YRK) had great impacts on partitioning of NH3-NH4+. The Multiple Linear Regression (MLR) model revealed a strong positive correlation between cation-NH4+ and anions-SO42− and NO3−. This research provides an insight into iPM2.5 formation mechanism for the advancement of PM2.5 control and regulation in the southeastern U.S. 
    more » « less