Inter-model comparison of global hydroxyl radical (OH) distributions and their impact on atmospheric methane over the 2000–2016 period

The modeling study presented here aims to estimate how uncertainties in global hydroxyl radical (OH) distributions, variability, and trends may contribute to resolving discrepancies between simulated and observed methane (CH4) changes since 2000. A multi-model ensemble of 14 OH fields was analyzed and aggregated into 64 scenarios to force the offline atmospheric chemistry transport model LMDz (Laboratoire de Meteorologie Dynamique) with a standard CH4 emission scenario over the period 2000–2016. The multi-model simulated global volume-weighted tropospheric mean OH concentration ([OH]) averaged over 2000–2010 ranges between 8:7*10^5 and 12:8*10^5 molec cm-3. The inter-model differences in tropospheric OH burden and vertical distributions are mainly determined by the differences in the nitrogen oxide (NO) distributions, while the spatial discrepancies between OH fields are mostly due to differences in natural emissions and volatile organic compound (VOC) chemistry. From 2000 to 2010, most simulated OH fields show an increase of 0.1–0:3*10^5 molec cm-3 in the tropospheric mean [OH], with year-to-year variations much smaller than during the historical period 1960–2000. Once ingested into the LMDz model, these OH changes translated into a 5 to 15 ppbv reduction in the CH4 mixing ratio in 2010, which represents 7%–20% of the model-simulated CH4 increase due to surface emissions. Between 2010 and 2016, the ensemble of simulations showed that OH changes could lead to a CH4 mixing ratio uncertainty of > 30 ppbv. Over the full 2000–2016 time period, using a common stateof- the-art but nonoptimized emission scenario, the impact of [OH] changes tested here can explain up to 54% of the gap between model simulations and observations. This result emphasizes the importance of better representing OH abundance and variations in CH4 forward simulations and emission optimizations performed by atmospheric inversions

The effect of atmospheric nudging on the stratospheric residual circulation in chemistry–climate models

We perform the first multi-model intercomparison of the impact of nudged meteorology on the stratospheric residual circulation using hindcast simulations from the Chemistry–Climate Model Initiative (CCMI). We examine simulations over the period 1980–2009 from seven models in which the meteorological fields are nudged towards a reanalysis dataset and compare these with their equivalent free-running simulations and the reanalyses themselves. We show that for the current implementations, nudging meteorology does not constrain the mean strength of the stratospheric residual circulation and that the inter-model spread is similar, or even larger, than in the free-running simulations. The nudged models generally show slightly stronger upwelling in the tropical lower stratosphere compared to the free-running versions and exhibit marked differences compared to the directly estimated residual circulation from the reanalysis dataset they are nudged towards. Downward control calculations applied to the nudged simulations reveal substantial differences between the climatological lower-stratospheric tropical upward mass flux (TUMF) computed from the modelled wave forcing and that calculated directly from the residual circulation. This explicitly shows that nudging decouples the wave forcing and the residual circulation so that the divergence of the angular momentum flux due to the mean motion is not balanced by eddy motions, as would typically be expected in the time mean. Overall, nudging meteorological fields leads to increased inter-model spread for most of the measures of the mean climatological stratospheric residual circulation assessed in this study. In contrast, the nudged simulations show a high degree of consistency in the inter-annual variability in the TUMF in the lower stratosphere, which is primarily related to the contribution to variability from the resolved wave forcing. The more consistent inter-annual variability in TUMF in the nudged models also compares more closely with the variability found in the reanalyses, particularly in boreal winter. We apply a multiple linear regression (MLR) model to separate the drivers of inter-annual and long-term variations in the simulated TUMF; this explains up to ∼75 % of the variance in TUMF in the nudged simulations. The MLR model reveals a statistically significant positive trend in TUMF for most models over the period 1980–2009. The TUMF trend magnitude is generally larger in the nudged models compared to their free-running counterparts, but the intermodel range of trends doubles from around a factor of 2 to a factor of 4 due to nudging. Furthermore, the nudged models generally do not match the TUMF trends in the reanalysis they are nudged towards for trends over different periods in the interval 1980–2009. Hence, we conclude that nudging does not strongly constrain long-term trends simulated by the chemistry–climate model (CCM) in the residual circulation. Our findings show that while nudged simulations may, by construction, produce accurate temperatures and realistic representations of fast horizontal transport, this is not typically the case for the slower zonal mean vertical transport in the stratosphere. Consequently, caution is required when using nudged simulations to interpret the behaviour of stratospheric tracers that are affected by the residual circulation

Evidence for a continuous decline in lower stratospheric ozone offsetting ozone layer recovery

Ozone forms in the Earth's atmosphere from the photodissociation of molecular oxygen, primarily in the tropical stratosphere. It is then transported to the extratropics by the Brewer–Dobson circulation (BDC), forming a protective "ozone layer" around the globe. Human emissions of halogen-containing ozone-depleting substances (hODSs) led to a decline in stratospheric ozone until they were banned by the Montreal Protocol, and since 1998 ozone in the upper stratosphere is rising again, likely the recovery from halogen-induced losses. Total column measurements of ozone between the Earth's surface and the top of the atmosphere indicate that the ozone layer has stopped declining across the globe, but no clear increase has been observed at latitudes between 60° S and 60° N outside the polar regions (60–90°). Here we report evidence from multiple satellite measurements that ozone in the lower stratosphere between 60° S and 60° N has indeed continued to decline since 1998. We find that, even though upper stratospheric ozone is recovering, the continuing downward trend in the lower stratosphere prevails, resulting in a downward trend in stratospheric column ozone between 60° S and 60° N. We find that total column ozone between 60° S and 60° N appears not to have decreased only because of increases in tropospheric column ozone that compensate for the stratospheric decreases. The reasons for the continued reduction of lower stratospheric ozone are not clear; models do not reproduce these trends, and thus the causes now urgently need to be established

Ozone sensitivity to varying greenhouse gases and ozone-depleting substances in CCMI-1 simulations

Ozone fields simulated for the first phase of the Chemistry-Climate Model Initiative (CCMI-1) will be used as forcing data in the 6th Coupled Model Intercomparison Project. Here we assess, using reference and sensitivity simulations produced for CCMI-1, the suitability of CCMI-1 model results for this process, investigating the degree of consistency amongst models regarding their responses to variations in individual forcings. We consider the influences of methane, nitrous oxide, a combination of chlorinated or brominated ozone-depleting substances, and a combination of carbon dioxide and other greenhouse gases. We find varying degrees of consistency in the models' responses in ozone to these individual forcings, including some considerable disagreement. In particular, the response of total-column ozone to these forcings is less consistent across the multi-model ensemble than profile comparisons. We analyse how stratospheric age of air, a commonly used diagnostic of stratospheric transport, responds to the forcings. For this diagnostic we find some salient differences in model behaviour, which may explain some of the findings for ozone. The findings imply that the ozone fields derived from CCMI-1 are subject to considerable uncertainties regarding the impacts of these anthropogenic forcings. We offer some thoughts on how to best approach the problem of generating a consensus ozone database from a multi-model ensemble such as CCMI-1

Deriving Global OH Abundance and Atmospheric Lifetimes for Long-Lived Gases: A Search for CH 3 CCl 3 Alternatives

An accurate estimate of global hydroxyl radical (OH) abundance is important for projections of air quality, climate, and stratospheric ozone recovery. As the atmospheric mixing ratios of methyl chloroform (CH₃CCl₃) (MCF), the commonly used OH reference gas, approaches zero, it is important to find alternative approaches to infer atmospheric OH abundance and variability. The lack of global bottom‐up emission inventories is the primary obstacle in choosing a MCF alternative. We illustrate that global emissions of long‐lived trace gases can be inferred from their observed mixing ratio differences between the Northern Hemisphere (NH) and Southern Hemisphere (SH), given realistic estimates of their NH‐SH exchange time, the emission partitioning between the two hemispheres, and the NH versus SH OH abundance ratio. Using the observed long‐term trend and emissions derived from the measured hemispheric gradient, the combination of HFC‐32 (CH₂F₂), HFC‐134a (CH₂FCF₃, HFC‐152a (CH₃CHF₂), and HCFC‐22 (CHClF₂), instead of a single gas, will be useful as a MCF alternative to infer global and hemispheric OH abundance and trace gas lifetimes. The primary assumption on which this multispecies approach relies is that the OH lifetimes can be estimated by scaling the thermal reaction rates of a reference gas at 272 K on global and hemispheric scales. Thus, the derived hemispheric and global OH estimates are forced to reconcile the observed trends and gradient for all four compounds simultaneously. However, currently, observations of these gases from the surface networks do not provide more accurate OH abundance estimate than that from MCF

Tyrosine kinase inhibitor therapy-induced changes in humoral immunity in patients with chronic myeloid leukemia

Purpose Tyrosine kinase inhibitors (TKIs) have well-characterized immunomodulatory effects on T and NK cells, but the effects on the humoral immunity are less well known. In this project, we studied TKI-induced changes in B cell-mediated immunity. Methods We collected peripheral blood (PB) and bone marrow (BM) samples from chronic myeloid leukemia (CML) patients before and during first-line imatinib (n = 20), dasatinib (n = 16), nilotinib (n = 8), and bosutinib (n = 12) treatment. Plasma immunoglobulin levels were measured, and different B cell populations in PB and BM were analyzed with flow cytometry. Results Imatinib treatment decreased plasma IgA and IgG levels, while dasatinib reduced IgM levels. At diagnosis, the proportion of patients with IgA, IgG, and IgM levels below the lower limit of normal (LLN) was 0, 11, and 6% of all CML patients, respectively, whereas at 12 months timepoint the proportions were 6% (p = 0.13), 31% (p = 0.042) and 28% (p = 0.0078). Lower initial Ig levels predisposed to the development of hypogammaglobulinemia during TKI therapy. Decreased Ig levels in imatinibtreated patients were associated with higher percentages of immature BM B cells. The patients, who had low Ig levels during the TKI therapy, had significantly more frequent minor infections during the follow-up compared with the patients with normal Ig values (33% vs. 3%, p = 0.0016). No severe infections were reported, except recurrent upper respiratory tract infections in one imatinib-treated patient, who developed severe hypogammaglobulinemia. Conclusions TKI treatment decreases plasma Ig levels, which should be measured in patients with recurrent infections.Peer reviewe

The representation of solar cycle signals in stratospheric ozone –Part 2: Analysis of global models

The impact of changes in incoming solar irradiance on stratospheric ozone abundances should be included in climate simulations to aid in capturing the atmospheric response to solar cycle variability. This study presents the first systematic comparison of the representation of the 11-year solar cycle ozone response (SOR) in chemistry–climate models (CCMs) and in pre-calculated ozone databases specified in climate models that do not include chemistry, with a special focus on comparing the recommended protocols for the Coupled Model Intercomparison Project Phase 5 and Phase 6 (CMIP5 and CMIP6). We analyse the SOR in eight CCMs from the Chemistry–Climate Model Initiative (CCMI-1) and compare these with results from three ozone databases for climate models: the Bodeker Scientific ozone database, the SPARC/Atmospheric Chemistry and Climate (AC&C) ozone database for CMIP5 and the SPARC/CCMI ozone database for CMIP6. The peak amplitude of the annual mean SOR in the tropical upper stratosphere (1–5hPa) decreases by more than a factor of 2, from around 5 to 2%, between the CMIP5 and CMIP6 ozone databases. This substantial decrease can be traced to the CMIP5 ozone database being constructed from a regression model fit to satellite and ozonesonde measurements, while the CMIP6 database is constructed from CCM simulations. The SOR in the CMIP6 ozone database therefore implicitly resembles the SOR in the CCMI-1 models. The structure in latitude of the SOR in the CMIP6 ozone database and CCMI-1 models is considerably smoother than in the CMIP5 database, which shows unrealistic sharp gradients in the SOR across the middle latitudes owing to the paucity of long-term ozone measurements in polar regions. The SORs in the CMIP6 ozone database and the CCMI-1 models show a seasonal dependence with enhanced meridional gradients at mid- to high latitudes in the winter hemisphere. The CMIP5 ozone database does not account for seasonal variations in the SOR, which is unrealistic. Sensitivity experiments with a global atmospheric model without chemistry (ECHAM6.3) are performed to assess the atmospheric impacts of changes in the representation of the SOR and solar spectral irradiance (SSI) forcing between CMIP5 and CMIP6. The larger amplitude of the SOR in the CMIP5 ozone database compared to CMIP6 causes a likely overestimation of the modelled tropical stratospheric temperature response between 11-year solar cycle minimum and maximum by up to 0.55K, or around 80% of the total amplitude. This effect is substantially larger than the change in temperature response due to differences in SSI forcing between CMIP5 and CMIP6. The results emphasize the importance of adequately representing the SOR in global models to capture the impact of the 11-year solar cycle on the atmosphere. Since a number of limitations in the representation of the SOR in the CMIP5 ozone database have been identified, we recommend that CMIP6 models without chemistry use the CMIP6 ozone database and the CMIP6 SSI dataset to better capture the climate impacts of solar variability. The SOR coefficients from the CMIP6 ozone database are published with this paper

The representation of solar cycle signals in stratospheric ozone – Part 2: Analysis of global models

The impact of changes in incoming solar irradiance on stratospheric ozone abundances should be included in climate simulations to aid in capturing the atmospheric response to solar cycle variability. This study presents the first systematic comparison of the representation of the 11-year solar cycle ozone response (SOR) in chemistry–climate models (CCMs) and in pre-calculated ozone databases specified in climate models that do not include chemistry, with a special focus on comparing the recommended protocols for the Coupled Model Intercomparison Project Phase 5 and Phase 6 (CMIP5 and CMIP6). We analyse the SOR in eight CCMs from the Chemistry–Climate Model Initiative (CCMI-1) and compare these with results from three ozone databases for climate models: the Bodeker Scientific ozone database, the SPARC/Atmospheric Chemistry and Climate (AC&C) ozone database for CMIP5 and the SPARC/CCMI ozone database for CMIP6. The peak amplitude of the annual mean SOR in the tropical upper stratosphere (1–5hPa) decreases by more than a factor of 2, from around 5 to 2%, between the CMIP5 and CMIP6 ozone databases. This substantial decrease can be traced to the CMIP5 ozone database being constructed from a regression model fit to satellite and ozonesonde measurements, while the CMIP6 database is constructed from CCM simulations. The SOR in the CMIP6 ozone database therefore implicitly resembles the SOR in the CCMI-1 models. The structure in latitude of the SOR in the CMIP6 ozone database and CCMI-1 models is considerably smoother than in the CMIP5 database, which shows unrealistic sharp gradients in the SOR across the middle latitudes owing to the paucity of long-term ozone measurements in polar regions. The SORs in the CMIP6 ozone database and the CCMI-1 models show a seasonal dependence with enhanced meridional gradients at mid- to high latitudes in the winter hemisphere. The CMIP5 ozone database does not account for seasonal variations in the SOR, which is unrealistic. Sensitivity experiments with a global atmospheric model without chemistry (ECHAM6.3) are performed to assess the atmospheric impacts of changes in the representation of the SOR and solar spectral irradiance (SSI) forcing between CMIP5 and CMIP6. The larger amplitude of the SOR in the CMIP5 ozone database compared to CMIP6 causes a likely overestimation of the modelled tropical stratospheric temperature response between 11-year solar cycle minimum and maximum by up to 0.55K, or around 80% of the total amplitude. This effect is substantially larger than the change in temperature response due to differences in SSI forcing between CMIP5 and CMIP6. The results emphasize the importance of adequately representing the SOR in global models to capture the impact of the 11-year solar cycle on the atmosphere. Since a number of limitations in the representation of the SOR in the CMIP5 ozone database have been identified, we recommend that CMIP6 models without chemistry use the CMIP6 ozone database and the CMIP6 SSI dataset to better capture the climate impacts of solar variability. The SOR coefficients from the CMIP6 ozone database are published with this paper

The influence of mixing on the stratospheric age of air changes in the 21st century

Climate models consistently predict an acceleration of the Brewer–Dobson circulation (BDC) due to climate change in the 21st century. However, the strength of this acceleration varies considerably among individual models, which constitutes a notable source of uncertainty for future climate projections. To shed more light upon the magnitude of this uncertainty and on its causes, we analyse the stratospheric mean age of air (AoA) of 10 climate projection simulations from the Chemistry-Climate Model Initiative phase 1 (CCMI-I), covering the period between 1960 and 2100. In agreement with previous multi-model studies, we find a large model spread in the magnitude of the AoA trend over the simulation period. Differences between future and past AoA are found to be predominantly due to differences in mixing (reduced aging by mixing and recirculation) rather than differences in residual mean transport. We furthermore analyse the mixing efficiency, a measure of the relative strength of mixing for given residual mean transport, which was previously hypothesised to be a model constant. Here, the mixing efficiency is found to vary not only across models, but also over time in all models. Changes in mixing efficiency are shown to be closely related to changes in AoA and quantified to roughly contribute 10 % to the longterm AoA decrease over the 21st century. Additionally, mixing efficiency variations are shown to considerably enhance model spread in AoA changes. To understand these mixing efficiency variations, we also present a consistent dynamical framework based on diffusive closure, which highlights the role of basic state potential vorticity gradients in controlling mixing efficiency and therefore aging by mixing