Solar Impact on Climate: Modeling the Coupling Between the Middle and the Lower Atmosphere

Solar variability influences the earth's atmosphere on different time scales. In particular, the impact of the 11-year solar cycle is of interest as it provides the major contribution to natural climate variability. Observations show clear 11-year variations in meteorological variables such as temperature or geopotential height from the upper atmosphere down to the troposphere and the earth's surface. In this paper the mechanisms will be discussed which are assumed to be responsible for the downward transfer of the solar signal within the atmosphere. These involve radiative, dynamical and chemical processes which have been studied in detail in model simulations and will be presented here

A critical evaluation of decadal solar cycle imprints in the MiKlip historical ensemble simulations

Studies concerning solar–terrestrial connections over the last decades claim to have found evidence that the quasi-decadal solar cycle can have an influence on the dynamics in the middle atmosphere in the Northern Hemisphere (NH) during the winter season. It has been argued that feedbacks between the intensity of the UV part of the solar spectrum and low-latitude stratospheric ozone may produce anomalies in meridional temperature gradients which have the potential to alter the zonal-mean flow in middle to high latitudes. Interactions between the zonal wind and planetary waves can lead to a downward propagation of the anomalies, produced in the middle atmosphere, down to the troposphere. More recently, it has been proposed that top-down-initiated decadal solar signals might modulate surface climate and synchronize the North Atlantic Oscillation. A realistic representation of the solar cycle in climate models was suggested to significantly enhance decadal prediction skill. These conclusions have been debated controversial since then due to the lack of realistic decadal prediction model setups and more extensive analysis. In this paper we aim for an objective and improved evaluation of possible solar imprints from the middle atmosphere to the surface and with that from head to toe. Thus, we analyze model output from historical ensemble simulations conducted with the state-of-the-art Max Planck Institute for Meteorology Earth System Model in high-resolution configuration (MPI-ESM-HR). The target of these simulations was to isolate the most crucial model physics to foster basic research on decadal climate prediction and to develop an operational ensemble decadal prediction system within the “Mittelfristige Klimaprognose” (MiKlip) framework. Based on correlations and multiple linear regression analysis we show that the MPI-ESM-HR simulates a realistic, statistically significant and robust shortwave heating rate and temperature response at the tropical stratopause, in good agreement with existing studies. However, the dynamical response to this initial radiative signal in the NH during the boreal winter season is weak. We find a slight strengthening of the polar vortex in midwinter during solar maximum conditions in the ensemble mean, which is consistent with the so-called “top-down” mechanism. The individual ensemble members, however, show a large spread in the dynamical response with opposite signs in response to the solar cycle, which might be a result of the large overall internal variability compensating for rather small solar imprints. We also analyze the possible surface responses to the 11-year solar cycle and review the proposed synchronization between the solar forcing and the North Atlantic Oscillation. We find that the simulated westerly wind anomalies in the lower troposphere, as well as the anomalies in the mean sea level pressure, are most likely independent from the timing of the solar signal in the middle atmosphere and the alleged top-down influences. The pattern rather reflects the decadal internal variability in the troposphere, mimicking positive and negative phases of the Arctic and North Atlantic oscillations throughout the year sporadically, which is then assigned to the solar predictor time series without any plausible physical connection and sound solar contribution. Finally, by applying lead–lag correlations, we find that the proposed synchronization between the solar cycle and the decadal component of the North Atlantic Oscillation might rather be a statistical artifact, affected for example by the internal decadal variability in the ocean, than a plausible physical connection between the UV solar forcing and quasi-decadal variations in the troposphere

Elevated stratopause events in the current and a future climate: A chemistry-climate model study

The characteristics and driving mechanisms of Elevated Stratopause Events (ESEs) are examined in simulations of the ECHAM/MESSy Atmospheric Chemistry (EMAC) chemistry-climate model under present and projected climate conditions. ESEs develop after sudden stratospheric warmings (SSWs) in boreal winter. While the stratopause descends during SSWs, it is reformed at higher altitudes after the SSWs, leading to ESEs in years with a particularly high new stratopause. EMAC reproduces well the frequency and main characteristics of observed ESEs. ESEs occur in 24% of the winters, mostly after major SSWs. They develop in stable polar vortices due to a persistent tropospheric wave forcing leading to a prolonged zonal wind reversal in the lower stratosphere. By wave filtering, this enables a faster re-establishment of the mesospheric westerly jet, polar downwelling and a higher stratopause. We find the presence of a westward-propagating wavenumber-1 planetary wave in the mesosphere following the onset, consistent with in-situ generation by large-scale instability. By the end of the 21st century, the number of ESEs is projected to increase, mainly due to a sinking of the original stratopause after strong tropospheric wave forcing and planetary wave dissipation at lower levels. Future ESEs develop preferably in more intense and cold polar vortices, and tend to be shorter. While in the current climate, planetary wavenumber-2 contributes to the forcing of ESEs, future wave forcing is dominated by wavenumber-1 activity as a result of climate change. Hence, a persistent wave forcing seems to be more relevant for the development of an ESE than the wavenumber decomposition of the forcing

Future changes in the stratosphere-to-troposphere ozone mass flux and the contribution from climate change and ozone recovery

Using a state-of-the-art chemistry–climate model we investigate the future change in stratosphere–troposphere exchange (STE) of ozone, the drivers of this change, as well as the future distribution of stratospheric ozone in the troposphere. Supplementary to previous work, our focus is on changes on the monthly scale. The global mean annual influx of stratospheric ozone into the troposphere is projected to increase by 53 % between the years 2000 and 2100 under the RCP8.5 greenhouse gas scenario. The change in ozone mass flux (OMF) into the troposphere is positive throughout the year with maximal increase in the summer months of the respective hemispheres. In the Northern Hemisphere (NH) this summer maximum STE increase is a result of increasing greenhouse gas (GHG) concentrations, whilst in the Southern Hemisphere(SH) it is due to equal contributions from decreasing levels of ozone depleting substances (ODS) and increasing GHG concentrations. In the SH the GHG effect is dominating in the winter months. A large ODS-related ozone increase in the SH stratosphere leads to a change in the seasonal breathing term which results in a future decrease of the OMF into the troposphere in the SH in September and October. The resulting distributions of stratospheric ozone in the troposphere differ for the GHG and ODS changes due to the following: (a) ozone input occurs at different regions for GHG- (midlatitudes) and ODS-changes (high latitudes); and (b) stratospheric ozone is more efficiently mixed towards lower tropospheric levels in the case of ODS changes, whereas tropospheric ozone loss rates grow when GHG concentrations rise. The comparison between the moderate RCP6.0 and the extreme RCP8.5 emission scenarios reveals that the annual global OMF trend is smaller in the moderate scenario, but the resulting change in the contribution of ozone with stratospheric origin (O3s) to ozone in the troposphere is of comparable magnitude in both scenarios. This is due to the larger tropospheric ozone precursor emissions and hence ozone production in the RCP8.5 scenario

Solar impact on the lower mesospheric subtropical jet: A comparative study with general circulation model simulations

The seasonal and interannual variation in the lower mesospheric subtropical jet (LMSJ) and their dependence on the 11-year solar cycle are studied by comparing observational data with simulations by two general circulation models. In the model simulations, a strengthening of the LMSJs is found in both hemispheres during the winter under the solar maximum condition, similar to the observation. However the model responses are substantially smaller except for one case in the southern hemisphere. It is also found that the stronger LMSJ due to an enhanced solar forcing appears during the period which follows an increasing period of interannual variation. Analysis of the observed seasonal march of the LMSJ in each year shows two different regimes of behavior. For a successful simulation, the model should realistically reproduce the observed interannual variability as well as the climatological mean

Quantifying uncertainties of climate signals in chemistry climate models related to the 11-year solar cycle – Part 1: Annual mean response in heating rates, temperature, and ozone

Variations in the solar spectral irradiance (SSI) with the 11-year sunspot cycle have been shown to have a significant impact on temperatures and the mixing ratios of atmospheric constituents in the stratosphere and mesosphere. Uncertainties in modelling the effects of SSI variations arise from uncertainties in the empirical models reconstructing the prescribed SSI data set as well as from uncertainties in the chemistry–climate model (CCM) formulation. In this study CCM simulations with the ECHAM/MESSy Atmospheric Chemistry (EMAC) model and the Community Earth System Model 1 (CESM1)–Whole Atmosphere Chemistry Climate Model (WACCM) have been performed to quantify the uncertainties of the solar responses in chemistry and dynamics that are due to the usage of five different SSI data sets or the two CCMs. We apply a two-way analysis of variance (ANOVA) to separate the influence of the SSI data sets and the CCMs on the variability of the solar response in shortwave heating rates, temperature, and ozone. The solar response is derived from climatological differences of time slice simulations prescribing SSI for the solar maximum in 1989 and near the solar minimum in 1994. The SSI values for the solar maximum of each SSI data set are created by adding the SSI differences between November 1994 and November 1989 to a common SSI reference spectrum for near-solar-minimum conditions based on ATLAS-3 (Atmospheric Laboratory of Applications and Science-3). The ANOVA identifies the SSI data set with the strongest influence on the variability of the solar response in shortwave heating rates in the upper mesosphere and in the upper stratosphere–lower mesosphere. The strongest influence on the variability of the solar response in ozone and temperature is identified in the upper stratosphere–lower mesosphere. However, in the region of the largest ozone mixing ratio, in the stratosphere from 50 to 10 hPa, the SSI data sets do not contribute much to the variability of the solar response when the Spectral And Total Irradiance REconstructions-T (SATIRE-T) SSI data set is omitted. The largest influence of the CCMs on variability of the solar responses can be identified in the upper mesosphere. The solar response in the lower stratosphere also depends on the CCM used, especially in the tropics and northern hemispheric subtropics and mid-latitudes, where the model dynamics modulate the solar responses. Apart from the upper mesosphere, there are also regions where the largest fraction of the variability of the solar response is explained by randomness, especially for the solar response in temperature

Elevated stratopause events in the current and a future climate: a chemistry-climate model study

The characteristics and driving mechanisms of Elevated Stratopause Events (ESEs) are examined in simulations of the ECHAM/MESSy Atmospheric Chemistry (EMAC) chemistry-climate model under present and projected climate conditions. ESEs develop after sudden stratospheric warmings (SSWs) in boreal winter. While the stratopause descends during SSWs, it is reformed at higher altitudes after the SSWs, leading to ESEs in years with a particularly high new stratopause. EMAC reproduces well the frequency and main characteristics of observed ESEs. ESEs occur in 24% of the winters, mostly after major SSWs. They develop in stable polar vortices due to a persistent tropospheric wave forcing leading to a prolonged zonal wind reversal in the lower stratosphere. By wave filtering, this enables a faster re-establishment of the mesospheric westerly jet, polar downwelling and a higher stratopause. We find the presence of a westward-propagating wavenumber-1 planetary wave in the mesosphere following the onset, consistent with in-situ generation by large-scale instability. By the end of the 21st century, the number of ESEs is projected to increase, mainly due to a sinking of the original stratopause after strong tropospheric wave forcing and planetary wave dissipation at lower levels. Future ESEs develop preferably in more intense and cold polar vortices, and tend to be shorter. While in the current climate, planetary wavenumber-2 contributes to the forcing of ESEs, future wave forcing is dominated by wavenumber-1 activity as a result of climate change. Hence, a persistent wave forcing seems to be more relevant for the development of an ESE than the wavenumber decomposition of the forcing

Chemical effects in 11-year solar cycle simulations with the Freie Universität Berlin Climate Middle Atmosphere Model with online chemistry (FUB-CMAM-CHEM)

The impact of 11-year solar cycle variations on stratospheric ozone (O3) is studied with the Freie Universität Berlin Climate Middle Atmosphere Model with interactive chemistry (FUB-CMAM-CHEM). To consider the effect of variations in charged particle precipitation we included an idealized NO x source in the upper mesosphere representing relativistic electron precipitation (REP). Our results suggest that the NO x source by particles and its transport from the mesosphere to the stratosphere in the polar vortex are important for the solar signal in stratospheric O3. We find a positive dipole O3 signal in the annual mean, peaking at 40–45 km at high latitudes and a negative O3 signal in the tropical lower stratosphere. This is similar to observations, but enhanced due to the idealized NO x source and at a lower altitude compared to the observed minimum. Our results imply that this negative O3 signal arises partly via chemical effects

The role of the stratosphere in Iberian Peninsula rainfall: a preliminary approach in February

This paper attempts to establish a connection between stratospheric anomalies in the North Pole and rainfall on the Iberian Peninsula through the occurrence of major midwinter warmings (MMWs) and cold events (CEs), taking February as a preliminary approach. We define the MMWs as the warmings which break down the polar vortex, whereas the CEs are the episodes in which the polar vortex remains cold and undisturbed. Both anomalies lead to a wind anomaly around the north polar stratosphere, which is connected with a shortly lagged tropospheric anomaly through a stratosphere-troposphere coupling in winter. A T-mode principal component analysis (PCA) was used as an objective pattern classification method for identifying the main daily surface-level pressure (SLP) patterns for February for the 1961-1990 reference period. Subsequently, those February months with an MMW or a CE influence in the troposphere are identified in the whole study period (1958-2000) by means of the Arctic Oscillation Index (AOI). Thus, performing the same analysis for the selected February months, new principal patterns for detecting changes in surface circulation structure and morphology are obtained. The results show a significant decrease in the westerlies and a southward shift of the storm tracks in Western Europe some weeks after an MMW occurrence, leading to an increase in precipitation in western Iberia and a slight decrease on the eastern Mediterranean fringe. The results are quite the opposite under a CE influence: the westerlies are strengthened and shifted northwards due to the displacement of the Atlantic anticyclone towards Central Europe; dry conditions are established throughout Iberia, except for the Mediterranean fringe, where precipitation shows a considerable increase due to the greater frequency of the northeasterly winds. Finally, an 11-year sunspot cycle-quasibiennial oscillation (QBO) modulation might be demonstrated in Iberian rainfall in February through the occurrence of these stratospheric anomalies

Solar forcing for CMIP6 (v3.2)

Abstract. This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models