The first application of a numerically exact, higher-order sensitivity analysis approach for atmospheric modelling: implementation of the hyperdual-step method in the Community Multiscale Air Quality Model (CMAQ) version 5.3.2
-
Published:2024-01-24
Issue:2
Volume:17
Page:567-585
-
ISSN:1991-9603
-
Container-title:Geoscientific Model Development
-
language:en
-
Short-container-title:Geosci. Model Dev.
Author:
Liu JiachenORCID, Chen Eric, Capps Shannon L.ORCID
Abstract
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.
Funder
Directorate for Engineering
Publisher
Copernicus GmbH
Reference47 articles.
1. Berman, B., Capps, S. L., Sauvageau, I., Gao, E., Eastham, S. D., and Russell, R. P.: ISORROPIA-MCX: Enabling Sensitivity Analysis With Multicomplex Variables in the Aerosol Thermodynamic Model, ISORROPIA, Earth Space Sci., 10, e2022EA002729, https://doi.org/10.1029/2022EA002729, 2023. 2. Boole, G.: A Treatise on the Calculus of Finite Difference, 2nd Edn., Dover, ISBN 048649523X, 1960. 3. Byun, D. and Schere, K. L.: Review of the Governing Equations, Computational Algorithms, and Other Components of the Models-3 Community Multiscale Air Quality (CMAQ) Modeling System, Appl. Mech. Rev., 59, 51–77, https://doi.org/10.1115/1.2128636, 2006. 4. Capaldo, K. P., Pilinis, C., and Pandis, S. N.: A computationally efficient hybrid approach for dynamic gas/aerosol transfer in air quality models, Atmos. Environ., 34, 3617–3627, https://doi.org/10.1016/S1352-2310(00)00092-3, 2000. 5. Che, W., Zheng, J., Wang, S., Zhong, L., and Lau, A.: Assessment of motor vehicle emission control policies using Model-3/CMAQ model for the Pearl River Delta region, China, Atmos. Environ., 45, 1740–1751, https://doi.org/10.1016/j.atmosenv.2010.12.050, 2011.
|
|