Data assimilation in atmospheric chemistry models: current status and future prospects for coupled chemistry meteorology models

Abstract. Data assimilation is used in atmospheric chemistry models to improve air quality forecasts, construct re-analyses of three-dimensional chemical (including aerosol) concentrations and perform inverse modeling of input variables or model parameters (e.g., emissions). Coupled chemistry meteorology models (CCMM) are atmospheric chemistry models that simulate meteorological processes and chemical transformations jointly. They offer the possibility to assimilate both meteorological and chemical data; however, because CCMM are fairly recent, data assimilation in CCMM has been limited to date. We review here the current status of data assimilation in atmospheric chemistry models with a particular focus on future prospects for data assimilation in CCMM. We first review the methods available for data assimilation in atmospheric models, including variational methods, ensemble Kalman filters, and hybrid methods. Next, we review past applications that have included chemical data assimilation in chemical transport models (CTM) and in CCMM. Observational data sets available for chemical data assimilation are described, including surface data, surface-based remote sensing, airborne data, and satellite data. Several case studies of chemical data assimilation in CCMM are presented to highlight the benefits obtained by assimilating chemical data in CCMM. A case study of data assimilation to constrain emissions is also presented. There are few examples to date of joint meteorological and chemical data assimilation in CCMM and potential difficulties associated with data assimilation in CCMM are discussed. As the number of variables being assimilated increases, it is essential to characterize correctly the errors; in particular, the specification of error cross-correlations may be problematic. In some cases, offline diagnostics are necessary to ensure that data assimilation can truly improve model performance. However, the main challenge is likely to be the paucity of chemical data available for assimilation in CCMM

A Kalman-filter bias correction of ozone deterministic, ensemble-averaged, and probabilistic forecasts

Kalman filtering (KF) is used to postprocess numerical-model output to estimate systematic errors in surface ozone forecasts. It is implemented with a recursive algorithm that updates its estimate of future ozone-concentration bias by using past forecasts and observations. KF performance is tested for three types of ozone forecasts: deterministic, ensemble-averaged, and probabilistic forecasts. Eight photochemical models were run for 56 days during summer 2004 over northeastern USA and southern Canada as part of the International Consortium for Atmospheric Research on Transport and Transformation New England Air Quality (AQ) Study. The raw and KF-corrected predictions are compared with ozone measurements from the Aerometric Information Retrieval Now data set, which includes roughly 360 surface stations. The completeness of the data set allowed a thorough sensitivity test of key KF parameters. It is found that the KF improves forecasts of ozone-concentration magnitude and the ability to predict rare events, both for deterministic and ensemble-averaged forecasts. It also improves the ability to predict the daily maximum ozone concentration, and reduces the time lag between the forecast and observed maxima. For this case study, KF considerably improves the predictive skill of probabilistic forecasts of ozone concentration greater than thresholds of 10 to 50 ppbv, but it degrades it for thresholds of 70 to 90 ppbv. Moreover, KF considerably reduces probabilistic forecast bias. The significance of KF postprocessing and ensemble-averaging is that they are both effective for real-time AQ forecasting. KF reduces systematic errors, whereas ensemble-averaging reduces random errors. When combined they produce the best overall forecast

Application of dynamic linear regression to improve the skill of ensemble-based deterministic ozone forecasts

Forecasts from seven air quality models and surface ozone data collected over the eastern USA and southern Canada during July and August 2004 provide a unique opportunity to assess benefits of ensemble-based ozone forecasting and devise methods to improve ozone forecasts. In this investigation, past forecasts from the ensemble of models and hourly surface ozone measurements at over 350 sites are used to issue deterministic 24-h forecasts using a method based on dynamic linear regression. Forecasts of hourly ozone concentrations as well as maximum daily 8-h and 1-h averaged concentrations are considered. It is shown that the forecasts issued with the application of this method have reduced bias and root mean square error and better overall performance scores than any of the ensemble members and the ensemble average. Performance of the method is similar to another method based on linear regression described previously by Pagowski et al., but unlike the latter, the current method does not require measurements from multiple monitors since it operates on individual time series. Improvement in the forecasts can be easily implemented and requires minimal computational cost

A novel method to improve temperature simulations of general circulation models based on ensemble empirical mode decomposition and its application to multi-model ensembles

A novel method based on the ensemble empirical mode decomposition (EEMD) method was developed to improve model performance. This method was evaluated by applying it to global surface air temperatures, which were simulated by eight general circulation models from the Coupled Model Intercomparison Project Phase 5 (CMIP5). The temperature simulations of the eight models were separated into their different components by EEMD. The model's performance improved after the first high-frequency component was removed from the original simulations by EEMD for each model, on both the global and continental scale. Moreover, EEMD was more effective in improving the model's performance compared to the wavelet transform method. The multi-model ensembles (MMEs) were calculated based on the EEMD-improved model simulations using the Average Ensemble Mean, Multiple Linear Regression, Singular Value Decomposition and Bayesian Model Averaging methods. The results showed that the MME forecasts performed better when the calculations were based on the EEMD-improved simulations as opposed to the original simulations on both the global and continental scale. Therefore, the results of the MME were further improved by using the EEMD-improved model simulations. This new method provides a simple way to improve model performance and can be easily applied to further improve MME simulations

One-Way Coupling of an Atmospheric and a Hydrologic Model in Colorado

This paper examines the accuracy of high-resolution nested mesoscale model simulations of surface climate. The nesting capabilities of the atmospheric fifth-generation Pennsylvania State University (PSU)–National Center for Atmospheric Research (NCAR) Mesoscale Model (MM5) were used to create high-resolution, 5-yr climate simulations (from 1 October 1994 through 30 September 1999), starting with a coarse nest of 20 km for the western United States. During this 5-yr period, two finer-resolution nests (5 and 1.7 km) were run over the Yampa River basin in northwestern Colorado. Raw and bias-corrected daily precipitation and maximum and minimum temperature time series from the three MM5 nests were used as input to the U.S. Geological Survey’s distributed hydrologic model [the Precipitation Runoff Modeling System (PRMS)] and were compared with PRMS results using measured climate station data. The distributed capabilities of PRMS were provided by partitioning the Yampa River basin into hydrologic response units (HRUs). In addition to the classic polygon method of HRU definition, HRUs for PRMS were defined based on the three MM5 nests. This resulted in 16 datasets being tested using PRMS. The input datasets were derived using measured station data and raw and bias-corrected MM5 20-, 5-, and 1.7-km output distributed to 1) polygon HRUs and 2) 20-, 5-, and 1.7-km-gridded HRUs, respectively. Each dataset was calibrated independently, using a multiobjective, stepwise automated procedure. Final results showed a general increase in the accuracy of simulated runoff with an increase in HRU resolution. In all steps of the calibration procedure, the station-based simulations of runoff showed higher accuracy than the MM5-based simulations, although the accuracy of MM5 simulations was close to station data for the high-resolution nests. Further work is warranted in identifying the causes of the biases in MM5 local climate simulations and developing methods to remove them.This article is from J. Hydrometeor, 7, 569–589. doi: http://dx.doi.org/10.1175/JHM512.1.</p

Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin

Viral erythrocytic necrosis (VEN) affects over 20 species of marine and anadromous fishes in the North Atlantic and North Pacific Oceans. However, the distribution and strain variation of its viral causative agent, erythrocytic necrosis virus (ENV), has not been well characterized within Pacific salmon. Here, metatranscriptomic sequencing of Chinook salmon revealed that ENV infecting salmon was closely related to ENV from Pacific herring, with inferred amino-acid sequences from Chinook salmon being 99% identical to those reported for herring. Sequence analysis also revealed 89 protein-encoding sequences attributed to ENV, greatly expanding the amount of genetic information available for this virus. High-throughput PCR of over 19,000 fish showed that ENV is widely distributed in the NE Pacific Ocean and was detected in 12 of 16 tested species, including in 27% of herring, 38% of anchovy, 17% of pollock, and 13% of sand lance. Despite frequent detection in marine fish, ENV prevalence was significantly lower in fish from freshwater (0.03%), as assessed with a generalized linear mixed effects model (p = 5.5 × 10⁻⁸). Thus, marine fish are likely a reservoir for the virus. High genetic similarity between ENV obtained from salmon and herring also suggests that transmission between these hosts is likely.Science, Faculty ofNon UBCBotany, Department ofEarth, Ocean and Atmospheric Sciences, Department ofMicrobiology and Immunology, Department ofOceans and Fisheries, Institute for theReviewedFacult