P. Räisänen

The Continual Intercomparison of Radiation Codes: Results from Phase I

[1] We present results from Phase I of the Continual Intercomparison of Radiation Codes (CIRC), intended as an evolving and regularly updated reference source for evaluation of radiative transfer (RT) codes used in global climate models and other atmospheric applications. CIRC differs from previous intercomparisons in that it relies on an observationally validated catalog of cases. The seven CIRC Phase I baseline cases, five cloud free and two with overcast liquid clouds, are built around observations by the Atmospheric Radiation Measurements program that satisfy the goals of Phase I, namely, to examine RT model performance in realistic, yet not overly complex, atmospheric conditions. Besides the seven baseline cases, additional idealized “subcases” are also employed to facilitate interpretation of model errors. In addition to quantifying individual model performance with respect to reference line-by-line calculations, we also highlight RT code behavior for conditions of doubled CO2, issues arising from spectral specification of surface albedo, and the impact of cloud scattering in the thermal infrared. Our analysis suggests that improvements in the calculation of diffuse shortwave flux, shortwave absorption, and shortwave CO2 forcing as well as in the treatment of spectral surface albedo should be considered for many RT codes. On the other hand, longwave calculations are generally in agreement with the reference results. By expanding the range of conditions under which participating codes are tested, future CIRC phases will hopefully allow even more rigorous examination of RT codes.

The Monte Carlo independent column approximation's conditional random noise : Impact on simulated climate

Abstract The Monte Carlo Independent Column Approximation (McICA) method for computing domain-average radiative fluxes is unbiased with respect to the full ICA, but its flux estimates contain conditional random noise. Results for five experiments are used to assess the impact of McICA-related noise on simulations of global climate made by the NCAR Community Atmosphere Model (CAM). The experiment with the least noise (an order of magnitude below that of basic McICA) is taken as the reference. Two additional experiments help demonstrate how the impact of noise depends on the time interval between calls to the radiation code. Each experiment is an ensemble of seven 15-month simulations. Experiments with very high noise levels feature significant reductions to cloudiness in the lowermost model layer over tropical oceans as well as changes in highly related quantities. This bias appears immediately, stabilizes after a couple of model days, and appears to stem from nonlinear interactions between clouds and radi...

Seasonal Changes in Solar Radiation and Relative Humidity in Europe in Response to Global Warming

AbstractFuture seasonal changes in surface incident solar radiation and relative humidity (RH) over Europe and adjacent ocean areas were assessed based on phase 3 of the Coupled Model Intercomparison Project (CMIP3) model ensemble. Under the A1B scenario, by 2070–99, summertime solar radiation is projected to increase by 5%–10% in central and southern Europe. In winter, radiation decreases in most of northern and eastern Europe by 5%–15%. RH drops in summer in the southern European inland by 8%–12%, whereas in winter a small increase of 2%–3% is projected for northeastern Europe. In spring, the change is an intermediate between those in the extreme seasons, while in autumn the patterns resemble summer. Over the northern Atlantic Ocean, RH increases in all seasons by 1%–2%. The intermodel agreement on the sign of all these shifts is good, and the patterns recur in the responses to the A2 and B1 scenarios. Substantial changes are already simulated to occur before the midcentury, for example, in summer RH de...

Climate impacts of changing aerosol emissions since 1996

Increases in Asian aerosol emissions have been suggested as one possible reason for the hiatus in global temperature increase during the past 15 years. We study the effect of sulphur and black carbon (BC) emission changes between 1996 and 2010 on the global energy balance. We find that the increased Asian emissions have had very little regional or global effects, while the emission reductions in Europe and the U.S. have caused a positive radiative forcing. In our simulations, the global-mean aerosol direct radiative effect changes by 0.06 W/m2 during 1996 to 2010, while the effective radiative forcing (ERF) is 0.42 W/m2. The rather large ERF arises mainly from changes in cloudiness, especially in Europe. In Asia, the BC warming due to sunlight absorption has largely offset the cooling caused by sulphate aerosols. Asian BC concentrations have increased by a nearly constant fraction at all altitudes, and thus, they warm the atmosphere also in cloudy conditions.

Sensitivity of the shortwave radiative effect of dust on particle shape: Comparison of spheres and spheroids

The sensitivity of direct shortwave radiative effects of dust (DRE) to assumed particle shape is investigated. Radiative transfer simulations are conducted using optical properties of either spheres, mass-equivalent spheroids (mass-conserving case), or (mass-equivalent) spheroids whose number concentration is modified so that they have the same midvisible optical thickness (tau(545 nm)) as spheres (tau-conserving case). The impact of particle shape on DRE is investigated for different dust particle effective radii, optical thickness of the dust cloud, solar zenith angle, and spectral surface albedo (ocean, grass, and desert). It is found that the influence of particle shape on the DRE is strongest over ocean. It also depends very strongly on the shape distribution of spheroids used, to a degree that the results for two distributions of spheroids may deviate more from each other than from those for spheres. Finally, the effects of nonsphericity largely depend on whether the mass- or tau-conserving case is considered. For example, when using a shape distribution of spheroids recommended in a recent study for approximating the single-scattering properties of dust, the DRE at the surface differs at most 5% from that from spherical particles in the mass-conserving case. This stems from compensating nonsphericity effects on optical thickness, asymmetry parameter, and single-scattering albedo. However, in the tau-conserving case, the negative DRE at the surface can be up to 15% weaker for spheroids than spheres.

Tests of Monte Carlo Independent Column Approximation in the ECHAM5 Atmospheric GCM

The Monte Carlo Independent Column Approximation (McICA) method for computing domain-average radiative fluxes allows a flexible treatment of unresolved cloud structure, and it is unbiased with respect to the full ICA, but its flux estimates contain conditional random noise. Here, tests of McICA in the ECHAM5 atmospheric GCM are reported. ECHAM5 provides an interesting test bed for McICA because it carries prognostic variables for the subgrid-scale probability distribution of total water content, which allows us to determine subgrid-scale cloud variability directly from the resolved-scale model variables. Three experiments with differing levels of radiative noise, each consisting of ten 6-yr runs, are performed to estimate the impact of McICA noise on simulated climate. In an experiment that attempted to deliberately maximize McICA noise, a systematic reduction in low cloud fraction occurred. For a more reasonable implementation of McICA, the impact of noise is very small, although statistically discernible. In terms of the impacts of noise, McICA appears to be a viable approach for use in ECHAM5. However, to improve the simulation of cloud radiative effects, realistic representation of both unresolved and resolved cloud structures is needed, which remains a challenging problem. Comparison of ECHAM5 data with a global cloud system–resolving model dataset and with International Satellite Cloud Climatology Project data suggested two problems related to unresolved cloud structures. First, ECHAM5 appears to underestimate subgrid-scale cloud variability. This problem seems partly related to the use of the beta distribution scheme for total water content in ECHAM5: in its current form, the scheme is unable to generate highly inhomogeneous clouds (relative standard deviation of condensate amount 1). Second, it appears that in ECHAM5, overcast cloud layers occur too frequently and partially cloudy layers too rarely. This problem is not unique to the beta distribution scheme; in fact, it is more pronounced when using an alternative, relative humidity–based cloud fraction scheme.

Seasonal soil moisture and drought occurrence in Europe in CMIP5 projections for the 21st century

Projections for near-surface soil moisture content in Europe for the 21st century were derived from simulations performed with 26 CMIP5 global climate models (GCMs). Two Representative Concentration Pathways, RCP4.5 and RCP8.5, were considered. Unlike in previous research in general, projections were calculated separately for all four calendar seasons. To make the moisture contents simulated by the various GCMs commensurate, the moisture data were normalized by the corresponding local maxima found in the output of each individual GCM. A majority of the GCMs proved to perform satisfactorily in simulating the geographical distribution of recent soil moisture in the warm season, the spatial correlation with an satellite-derived estimate varying between 0.4 and 0.8. In southern Europe, long-term mean soil moisture is projected to decline substantially in all seasons. In summer and autumn, pronounced soil drying also afflicts western and central Europe. In northern Europe, drying mainly occurs in spring, in correspondence with an earlier melt of snow and soil frost. The spatial pattern of drying is qualitatively similar for both RCP scenarios, but weaker in magnitude under RCP4.5. In general, those GCMs that simulate the largest decreases in precipitation and increases in temperature and solar radiation tend to produce the most severe soil drying. Concurrently with the reduction of time-mean soil moisture, episodes with an anomalously low soil moisture, occurring once in 10 years in the recent past simulations, become far more common. In southern Europe by the late 21st century under RCP8.5, such events would be experienced about every second year.

Improved power-law estimates from multiple samples provided by millennium climate simulations

Using the long annual mean temperature time series provided by millennium Earth System Model simulations and a method of discrete Fourier transform with varying starting point and length of time window together with averaging, we get good fits to power laws between two characteristic oscillatory timescales of the model climate: multidecadal (50–80 years) and El Nino (3–6 years) timescales. For global mean temperature, we fit β ∼ 0.35 in a relation S(f) ∼ f−β in a simulation without external climate forcing and β over 0.7 in a simulation with external forcing included. The power law is found both with and without external forcing despite the forcings, e.g. the volcanic forcing, not showing similar behaviour, indicating a nonlinear temperature response to time-varying forcing. We also fit a power law with β ∼ 8 to the narrow frequency range between El Nino frequencies (up to 1/(3.2 years)) and the Nyquist frequency (1/(2 years)). Also, monthly mean temperature time series are considered and a decent power-law fit for frequencies above 1/year is obtained. Regional variability in best-fit β is explored, and the impact of choosing the frequency range on the result is illustrated. When all resolved frequencies are used, land areas seem to have lower βs than ocean areas on average, but when fits are restricted to frequencies below 1/(6 years), this difference disappears, while regional differences still remain. Results compare well with measurements both for global mean temperature and for the central England temperature record.

Effect of aerosol size distribution changes on AOD, CCN and cloud droplet concentration: Case studies from Erfurt and Melpitz, Germany

[1] For the period of 1990 to 2000, atmospheric particulate mass concentrations have decreased in Central Europe. Simultaneously, the amount of shortwave radiation reaching the ground increased during clear sky conditions. The aerosol indirect effect has not been seen as clearly, as the radiation reaching the ground during overcast conditions has not increased as much as might be expected. Here we show that this may be caused by the condensation kinetics of water during cloud droplet formation. The decrease in the particulate mass led to a clear decrease in the number concentration of cloud condensation nuclei (CCN). However, in urban areas a relatively larger decrease in the number of particles in the upper end of the accumulation mode has led to slower condensation of water. As a result, a higher maximum supersaturation is reached during the cloud droplet formation. This compensates for the effect of decreased CCN concentrations. For example in Erfurt between 1991 and 1996, the aerosol properties changed so that aerosol optical depth decreased by 58% and CCN concentration decreased by 25 to 50%. These led to a 4 to 12% reduction in cloud droplet number concentration (CDNC) and less than a 2 Wm � 2 increase in shortwave radiation during overcast conditions. These results demonstrate that locally the aerosol direct effect can be much larger than the aerosol indirect effect. Furthermore, even though AOD appears to be a valid proxy for CCN, the correlation between AOD and CDNC is not straightforward and thus AOD cannot be used as a proxy for CDNC.

Anthropogenic aerosol forcing n insights from multi-estimates from aerosol-climate models with reduced complexity

The radiative forcing of anthropogenic aerosol remains a key uncertainty in the understanding of climate change. This study quantifies the model spread in aerosol forcing associated with (i) variability internal to the atmosphere and (ii) differences in the model representation of weather. We do so by performing ensembles of atmosphere-only simulations with four state-of-the-art Earth system models, three of which will be used in the sixth coupled model inter-comparison project (CMIP6, Eyring et al., 2016). In those models we reduce the complexity of the anthropogenic aerosol by prescribing the same 5 annually-repeating patterns of the anthropogenic aerosol optical properties and associated effects on the cloud reflectivity. We quantify a comparably small model spread in the long-term averaged ERF compared to the overall possible range in annual ERF estimates associated with model-internal variability. This implies that identifying the true model spread in ERF associated with differences in the representation of meteorological processes and natural aerosol requires averaging over a sufficiently large number of annual estimates. We characterize the model diversity in clouds and use satellite products as benchmarks. 10 Despite major inter-model differences in natural aerosol and clouds, all models show only a small change in the global-mean ERF due to the substantial change in the global anthropogenic aerosol distribution between the mid-1970s and mid-2000s, the ensemble mean ERF being -0.47 Wm−2 for the mid-1970s and -0.51 Wm−2 for the mid-2000s. This result suggests that inter-comparing ERF changes between two periods rather than absolute magnitudes relative to pre-industrial might provide a more stringent test for a model’s ability for representing climate evolutions. 15