scholarly journals Advanced error diagnostics of the CMAQ and Chimere modelling systems within the AQMEII3 model evaluation framework

Author(s):  
Efisio Solazzo ◽  
Christian Hogrefe ◽  
Augustin Colette ◽  
Marta Garcia-Vivanco ◽  
Stefano Galmarini

Abstract. The work here complements the overview analysis of the modelling systems participating in the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) by focusing on the performance for hourly surface ozone by two modelling systems, Chimere for Europe and CMAQ for North America. The evaluation strategy outlined in the course of the three phases of the AQMEII activity, aimed to build up a diagnostic methodology for model evaluation, is pursued here and novel diagnostic methods are proposed. In addition to evaluating the base case simulation in which all model components are configured in their standard mode, the analysis also makes use of sensitivity simulations in which the models have been applied by altering and/or zeroing lateral boundary conditions, emissions of anthropogenic precursors, and ozone dry deposition. To help understand of the causes of model deficiencies, the error components (bias, variance, and covariance) of the base case and of the sensitivity runs are analysed in conjunction with time-scale considerations and error modelling using the available error fields of temperature, wind speed, and NOx concentration. The results reveal the effectiveness and diagnostic power of the methods devised (which remains the main scope of this study), allowing the detection of the time scale and the fields that the two models are most sensitive to. The representation of planetary boundary layers (PBL) dynamics is pivotal to both models. In particular: i) The fluctuations slower than −1.5 days account for 70–85 % of the total ozone quadratic error; ii) A recursive, systematic error with daily periodicity is detected, responsible for 10–20 % of the quadratic total error; iii) Errors in representing the timing of the daily transition between stability regimes in the PBL are responsible for a covariance error as large as 9 ppb (as much as the standard deviation of the network-average ozone observations in summer in both Europe and North America); iv) The CMAQ ozone error has a weak/negligible dependence on the errors in NO2 and wind speed, while the error in NO2 significantly impacts the ozone error produced by Chimere; v) On a continent wide monitoring network-average, a zeroing out of anthropogenic emissions produces an error increase of 45 % (25 %) during summer and of 56 % (null) during winter for Chimere (CMAQ), while a zeroing out of lateral boundary conditions results in an ozone error increase of 30 % during summer and of 180 % during winter (CMAQ).

2017 ◽  
Vol 17 (17) ◽  
pp. 10435-10465 ◽  
Author(s):  
Efisio Solazzo ◽  
Christian Hogrefe ◽  
Augustin Colette ◽  
Marta Garcia-Vivanco ◽  
Stefano Galmarini

Abstract. The work here complements the overview analysis of the modelling systems participating in the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) by focusing on the performance for hourly surface ozone by two modelling systems, Chimere for Europe and CMAQ for North America. The evaluation strategy outlined in the course of the three phases of the AQMEII activity, aimed to build up a diagnostic methodology for model evaluation, is pursued here and novel diagnostic methods are proposed. In addition to evaluating the base case simulation in which all model components are configured in their standard mode, the analysis also makes use of sensitivity simulations in which the models have been applied by altering and/or zeroing lateral boundary conditions, emissions of anthropogenic precursors, and ozone dry deposition. To help understand of the causes of model deficiencies, the error components (bias, variance, and covariance) of the base case and of the sensitivity runs are analysed in conjunction with timescale considerations and error modelling using the available error fields of temperature, wind speed, and NOx concentration. The results reveal the effectiveness and diagnostic power of the methods devised (which remains the main scope of this study), allowing the detection of the timescale and the fields that the two models are most sensitive to. The representation of planetary boundary layer (PBL) dynamics is pivotal to both models. In particular, (i) the fluctuations slower than ∼ 1.5 days account for 70–85 % of the mean square error of the full (undecomposed) ozone time series; (ii) a recursive, systematic error with daily periodicity is detected, responsible for 10–20 % of the quadratic total error; (iii) errors in representing the timing of the daily transition between stability regimes in the PBL are responsible for a covariance error as large as 9 ppb (as much as the standard deviation of the network-average ozone observations in summer in both Europe and North America); (iv) the CMAQ ozone error has a weak/negligible dependence on the errors in NO2, while the error in NO2 significantly impacts the ozone error produced by Chimere; (v) the response of the models to variations of anthropogenic emissions and boundary conditions show a pronounced spatial heterogeneity, while the seasonal variability of the response is found to be less marked. Only during the winter season does the zeroing of boundary values for North America produce a spatially uniform deterioration of the model accuracy across the majority of the continent.


2012 ◽  
Vol 140 (11) ◽  
pp. 3525-3551 ◽  
Author(s):  
Michael S. Buban ◽  
Conrad L. Ziegler ◽  
Edward R. Mansell ◽  
Yvette P. Richardson

Abstract A dryline and misocyclones have been simulated in a cloud-resolving model by applying specified initial and time-dependent lateral boundary conditions obtained from analyses of the 22 May 2002 International H2O Project (IHOP_2002) dataset. The initial and lateral boundary conditions were obtained from a combination of the time–spaced Lagrangian analyses for temperature and moisture with horizontal velocities from multiple-Doppler wind syntheses. The simulated dryline, horizontal dry-convective rolls (HCRs) and open cells (OCCs), misocyclones, and cumulus clouds are similar to the corresponding observed features. The misocyclones move northward at nearly the mean boundary layer (BL) wind speed, rotate dryline gradients owing to their circulations, and move the local dryline eastward via their passage. Cumuli develop along a secondary dryline, along HCR and OCC segments between the primary and secondary drylines, along HCR and OCC segments that have moved over the dryline, and within the dryline updraft. After the initial shearing instability develops, misocyclogenesis proceeds from tilting and stretching of vorticity by the persistent secondary dryline circulation. The resulting misocyclone evolution is discussed.


2015 ◽  
Vol 15 (4) ◽  
pp. 5327-5358 ◽  
Author(s):  
Z. Jiang ◽  
D. B. A. Jones ◽  
J. Worden ◽  
H. M. Worden ◽  
D. K. Henze ◽  
...  

Abstract. Chemical transport models (CTMs) driven with high-resolution meteorological fields can better resolve small-scale processes, such as frontal lifting or deep convection, and thus improve the simulation and emission estimates of tropospheric trace gases. In this work, we explore the use of the GEOS-Chem four-dimensional variational (4-D-Var) data assimilation system with the nested high-resolution version of the model (0.5° × 0.67°) to quantify North American CO emissions during the period of June 2004 – May 2005. With optimized lateral boundary conditions, regional inversion analyses can reduce the sensitivity of the CO source estimates to errors in long-range transport and in the distributions of the hydroxyl radical (OH), the main sink for CO. To further limit the potential impact of discrepancies in chemical aging of air in the free troposphere, associated with errors in OH, we use surface level multispectral MOPITT CO retrievals, which have greater sensitivity to CO near the surface and reduced sensitivity in the free troposphere, compared to previous versions of the retrievals. We estimate that the annual total anthropogenic CO emission from the contiguous US 48 states was 97 Tg CO, a 14% increase from the 85 Tg CO in the a priori. This increase is mainly due to enhanced emissions around the Great Lakes region and along the west coast, relative to the a priori. Sensitivity analyses using different OH fields and lateral boundary conditions suggest a possible error, associated with local North America OH distribution, in these emission estimates of 20% during summer 2004, when the CO lifetime is short. This 20% OH-related error is 50% smaller than the OH-related error previously estimated for North American CO emissions using a global inversion analysis. We believe that reducing this OH-related error further will require integrating additional observations to provide a strong constraint on the CO distribution across the domain. Despite these limitations, our results show the potential advantages of combining high-resolution regional inversion analyses with global analyses to better quantify regional CO source estimates.


2018 ◽  
Author(s):  
Marina Astitha ◽  
Ioannis Kioutsoukis ◽  
Ghezae Araya Fisseha ◽  
Roberto Bianconi ◽  
Johannes Bieser ◽  
...  

Abstract. This study evaluates simulated vertical ozone profiles produced in the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) against ozonesonde observations in North America for the year 2010. Four research groups from the United States (U.S.) and Europe have provided ozone vertical profiles to conduct this analysis. Because some of the modeling systems differ in their meteorological drivers, wind speed and temperature are also included in the analysis. In addition to the seasonal ozone profile evaluation for 2010, we also analyze chemically inert tracers designed to track the influence of lateral boundary conditions on simulated ozone profiles within the modeling domain. Finally, cases of stratospheric ozone intrusions during May–June 2010 are investigated by analyzing ozonesonde measurements and the corresponding model simulations at Intercontinental Chemical Transport Experiment Ozonesonde Network Study (IONS) experiment sites in the western United States. The evaluation of the seasonal ozone profiles reveals that at a majority of the stations, ozone mixing ratios are under-estimated in the 1–6 km range. The seasonal change noted in the errors follows the one seen in the variance of ozone mixing ratios, with the majority of the models exhibiting less variability than the observations. The analysis of chemically inert tracers highlights the importance of lateral boundary conditions up to 250 hPa for the lower tropospheric ozone mixing ratios (0–2 km). Finally, for the stratospheric intrusions, the models are generally able to reproduce the location and timing of most intrusions but underestimate the magnitude of the maximum mixing ratios in the 2–6 km range and overestimate ozone up to the first km possibly due to marine air influences that are not accurately described by the models. The choice of meteorological driver appears to be a greater predictor of model skill in this altitude range than the choice of air quality model.


2002 ◽  
Vol 35 ◽  
pp. 503-509 ◽  
Author(s):  
Olivier Gagliardini ◽  
Jacques Meyssonnier

AbstractA local two-dimensional flow model which accounts for the anisotropic behaviour of polar ice and the evolution of its strain-induced anisotropy is briefly reviewed. Due to its complexity, it is not yet possible to use this model to simulate the flow of a whole ice sheet, and its potential applications are presently restricted to limited spatial domains around existing drilling sites. In order to calculate the local flow of ice, boundary conditions must be applied on the lateral edges of the studied domain. Since these limits correspond to fictitious sections of the ice sheet, the type of boundary condition to adopt is not obvious. In the present paper, different kinds of boundary conditions of the Dirichlet type, applied at the lateral boundary of an idealized ice sheet of simplified geometry, are discussed. This will serve as a first step towards the coupling of the local flow model with a global ice-sheet flow model.


2021 ◽  
Author(s):  
Mehmet Ilicak ◽  
Ivan Federico ◽  
Ivano Barletta ◽  
Nadia Pinardi ◽  
Stefania Angela Ciliberti ◽  
...  

<p>Marmara Sea including Bosphorus and Dardanelles Straits (i.e. Turkish Strait Systems, TSS) is the connection between the Black Sea and the Mediterranean. The exchange flow that occurs in the Straits is crucial to set the deep water properties in the Black Sea and the surface water conditions in the Northern Aegean Sea. We have developed a new high-resolution unstructured grid model (U-TSS) for the Marmara Sea including the Bosporus and Dardanelles Straits using the System of HydrodYnamic Finite Element Modules (SHYFEM). Using an unstructured grid in the horizontal better resolves geometry of the Turkish Straits. The new model has a resolution between 500 meter in the deep to 50 meter in the shallow areas, and 93 geopotential coordinate levels in the vertical. We conducted a 4 year hindcast simulation between 2016 and 2019 using lateral boundary conditions from CMEMS (https://marine.copernicus.eu/) analysis, in particular Black Sea Forecasting System (BS-FS) for the northern boundary and Mediterranean Sea Forecasting System (MS-FS) for the southern boundary. Atmospheric boundary conditions fare from the ECMWF dataset.</p><p>Mean averaged surface circulation shows that there is a cyclonic gyre in the middle of the basin due to Bosphorus jet flowing to the south and turning to west after reaching the southern Marmara coast. The U-TSS model has been validated against the seasonal in situ observations obtained from four different cruises between 2017 and 2018. The maximum bias occurs at around halocline depth between 20 to 30 meters.  We also found that root mean square error field is higher in the mixed layer interface. We conclude that capturing shallow mixed layer depth is very in the Marmara Sea due to the interplay of air-sea fluxes and mixing parametrizations uncertainties. Maximum salinity bias and rms in the new U-TSS model are around 3 psu which is a significant improvement with respect to previous studies. This new model will be used as an operational forecasting system and will provide lateral boundary conditions for the BS-FS and MS-FS models in the future.</p>


Sign in / Sign up

Export Citation Format

Share Document