A. V. Eliseev

CLIMBER-2: A climate system model of intermediate complexity. Part I: Model description and performance for present climate

Abstract A 2.5-dimensional climate system model of intermediate complexity CLIMBER-2 and its performance for present climate conditions are presented. The model consists of modules describing atmosphere, ocean, sea ice, land surface processes, terrestrial vegetation cover, and global carbon cycle. The modules interact through the fluxes of momentum, energy, water and carbon. The model has a coarse spatial resolution, nevertheless capturing the major features of the Earths geography. The model describes temporal variability of the system on seasonal and longer time scales. Due to the fact that the model does not employ flux adjustments and has a fast turnaround time, it can be used to study climates significantly different from the present one and to perform long-term (multimillennia) simulations. The comparison of the model results with present climate data show that the model successfully describes the seasonal variability of a large set of characteristics of the climate system, including radiative balance, temperature, precipitation, ocean circulation and cryosphere.

Long-Term Climate Change Commitment and Reversibility: An EMIC Intercomparison

AbstractThis paper summarizes the results of an intercomparison project with Earth System Models of Intermediate Complexity (EMICs) undertaken in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). The focus is on long-term climate projections designed to 1) quantify the climate change commitment of different radiative forcing trajectories and 2) explore the extent to which climate change is reversible on human time scales. All commitment simulations follow the four representative concentration pathways (RCPs) and their extensions to year 2300. Most EMICs simulate substantial surface air temperature and thermosteric sea level rise commitment following stabilization of the atmospheric composition at year-2300 levels. The meridional overturning circulation (MOC) is weakened temporarily and recovers to near-preindustrial values in most models for RCPs 2.6–6.0. The MOC weakening is more persistent for RCP8.5. Elimination of anthropogenic CO2 emissions after 2300 resu...

Interaction of the Methane Cycle and Processes in Wetland Ecosystems in a Climate Model of Intermediate Complexity

The climate model of the Institute of Atmospheric Physics of the Russian Academy of Sciences (IAP RAS CM) has been supplemented with a module of soil thermal physics and the methane cycle, which takes into account the response of methane emissions from wetland ecosystems to climate changes. Methane emissions are allowed only from unfrozen top layers of the soil, with an additional constraint in the depth of the simulated layer. All wetland ecosystems are assumed to be water-saturated. The molar amount of the methane oxidized in the atmosphere is added to the simulated atmospheric concentration of CO2. A control preindustrial experiment and a series of numerical experiments for the 17th–21st centuries were conducted with the model forced by greenhouse gases and tropospheric sulfate aerosols. It is shown that the IAP RAS CM generally reproduces preindustrial and current characteristics of both seasonal thawing/freezing of the soil and the methane cycle. During global warming in the 21st century, the permafrost area is reduced by four million square kilometers. By the end of the 21st century, methane emissions from wetland ecosystems amount to 130–140 Mt CH4/year for the preindustrial and current period increase to 170–200 MtCH4/year. In the aggressive anthropogenic forcing scenario A2, the atmospheric methane concentration grows steadily to ≈3900 ppb. In more moderate scenarios A1B and B1, the methane concentration increases until the mid-21st century, reaching ≈2100–2400 ppb, and then decreases. Methane oxidation in air results in a slight additional growth of the atmospheric concentration of carbon dioxide. Allowance for the interaction between processes in wetland ecosystems and the methane cycle in the IAP RAS CM leads to an additional atmospheric methane increase of 10–20% depending on the anthropogenic forcing scenario and the time. The causes of this additional increase are the temperature dependence of integral methane production and the longer duration of a warm period in the soil. However, the resulting enhancement of the instantaneous greenhouse radiative forcing of atmospheric methane and an increase in the mean surface air temperature are small (globally < 0.1 W/m2 and 0.05 K, respectively).

Decadal climate variability in a coupled atmosphere-ocean climate model of moderate complexity

In this study we determined characteristic temporal modes of atmospheric variability at the decadal and interdecadal timescales. This was done on the basis of 1000 year long integrations of a global coupled atmosphere-ocean climate model of moderate complexity including the troposphere, stratosphere, and mesosphere. The applied model resolves explicitely the basic features of the large-scale long-term atmospheric and oceanic variables. The synoptic-scale processes are described in terms of autocorrelation and crosscorrelation functions. The paper includes an extended description and validation of the model as well as the results of analyses of two 1000 year long model integrations. One model run has been performed with the fully coupled model of the atmosphere-ocean system. The performed time-frequency analyses of atmospheric fields reveal strong decadal and interdecadal modes with periods of about 9, 18, and 30 years. To quantify the influence of the ocean on atmospheric variations an additional run with seasonally varying prescribed sea surface temperatures has been carried out, which is characterized by strong decadal modes with periods of about 9 years. The comparison of both runs suggests that decadal variability can be understood as an inherent atmospheric mode due to the nonlinear dynamics of the large-scale atmospheric circulation patterns whereas interdecadal climate variability has to be regarded as coupled atmosphere-ocean modes.

Future changes in extratropical storm tracks and baroclinicity under climate change

The weather in Eurasia, Australia, and North and South America is largely controlled by the strength and position of extratropical storm tracks. Future climate change will likely affect these storm tracks and the associated transport of energy, momentum, and water vapour. Many recent studies have analyzed how storm tracks will change under climate change, and how these changes are related to atmospheric dynamics. However, there are still discrepancies between different studies on how storm tracks will change under future climate scenarios. Here, we show that under global warming the CMIP5 ensemble of coupled climate models projects only little relative changes in vertically averaged mid-latitude mean storm track activity during the northern winter, but agree in projecting a substantial decrease during summer. Seasonal changes in the Southern Hemisphere show the opposite behaviour, with an intensification in winter and no change during summer. These distinct seasonal changes in northern summer and southern winter storm tracks lead to an amplified seasonal cycle in a future climate. Similar changes are seen in the mid-latitude mean Eady growth rate maximum, a measure that combines changes in vertical shear and static stability based on baroclinic instability theory. Regression analysis between changes in the storm tracks and changes in the maximum Eady growth rate reveal that most models agree in a positive association between the two quantities over mid-latitude regions.

Amplitude-Phase Characteristics of the Annual Cycle of Surface Air Temperature in the Northern Hemisphere

The amplitude phase characteristics (APC) or surface air temperature (SAT) annual cycle (AC) in the Northern Hemisphere are analyzed. From meteorological observations for the 20th century and meteorological reanalyses for its second half, it is found that over land negative correlation of SAT AC amplitude with annual mean SAT dominates. Nevertheless, some exceptions exist. The positive correlation between these two variables is found over the two desert regions: in northern Africa and in Central America. Areas of positive correlations are also found for the northern Pacific and for the tropical Indian and Pacific Oceans. Southward of the characteristic annual mean snow ice boundary (SIB) position, the shape of the SAT AC becomes more sinusoidal under climate warming. In contrast, northward of it, this shape becomes less sinusoidal. The latter is also found for the above-mentioned two desert regions. In the Far East (southward of about 50°N), the SAT AC shifts as a whole: here its spring and autumn phases occur earlier if the annual mean SAT increases. From energy balance climate considerations, those trends for SAT AC APC in the middle and high latitudes are associated with the influence of the albedo SAT feedback due to the SIB movement. In the Far East the trends are attributed to the interannual cloudiness variability, and in the desert regions, to the influence of a further desertification and/or scattering aerosol loading into the atmosphere. In the north Pacific, the exhibited trends could only be explained as a result of the influence of the greenhouse-gases loading on atmospheric opacity. The trends for SAT AC APC related to the SIB movement are simulated reasonably well by the climate model of intermediate complexity (IAP RAS CM) in the experiment with greenhouse gases atmospheric loading. In contrast, the tendencies resulting from the cloudiness variability are not reproduced by this model. The model also partly simulates the tendencies related to the desertification processes.

Simulation of Characteristics of Thermal and Hydrologic Soil Regimes in Equilibrium Numerical Experiments with a Climate Model of Intermediate Complexity

The IAP RAS CM (Institute of Atmospheric Physics, Russian Academy of Sciences, climate model) has been extended to include a comprehensive scheme of thermal and hydrologic soil processes. In equilibrium numerical experiments with specified preindustrial and current concentrations of atmospheric carbon dioxide, the coupled model successfully reproduces thermal characteristics of soil, including the temperature of its surface, and seasonal thawing and freezing characteristics. On the whole, the model also reproduces soil hydrology, including the winter snow water equivalent and river runoff from large watersheds. Evapotranspiration from the soil surface and soil moisture are simulated somewhat worse. The equilibrium response of the model to a doubling of atmospheric carbon dioxide shows a considerable warming of the soil surface, a reduction in the extent of permanently frozen soils, and the general growth of evaporation from continents. River runoff increases at high latitudes and decreases in the subtropics. The results are in qualitative agreement with observational data for the 20th century and with climate model simulations for the 21st century.

Changes in climatic characteristics of Northern Hemisphere extratropical land in the 21st century: Assessments with the IAP RAS climate model

Assessments of future changes in the climate of Northern Hemisphere extratropical land regions have been made with the IAP RAS climate model (CM) of intermediate complexity (which includes a detailed scheme of thermo- and hydrophysical soil processes) under prescribed greenhouse and sulfate anthropogenic forcing from observational data for the 19th and 20th centuries and from the SRES B1, A1B, and A2 scenarios for the 21st century. The annual mean warming of the extratropical land surface has been found to reach 2–5 K (3–10 K) by the middle (end) of the 21st century relative to 1961–1990, depending on the anthropogenic forcing scenario, with larger values in North America than in Europe. Winter warming is greater than summer warming. This is expressed in a decrease of 1–4 K (or more) in the amplitude of the annual harmonic of soil-surface temperature in the middle and high latitudes of Eurasia and North America. The total area extent of perennially frozen ground Sp in the IAP RAS CM changes only slightly until the late 20th century, reaching about 21 million km2, and then decreases to 11–12 million km2 in 2036–2065 and 4–8 million km2 in 2071–2100. In the late 21st century, near-surface permafrost is expected to remain only in Tibet and in central and eastern Siberia. In these regions, depths of seasonal thaw exceed 1 m (2 m) under the SRES B1 (A1B or A2) scenario. The total land area with seasonal thaw or cooling is expected to decrease from the current value of 54–55 million km2 to 38–42 in the late 21st century. The area of Northern Hemisphere snow cover in February is also reduced from the current value of 45–49 million km2 to 31–37 million km2. For the basins of major rivers in the extratropical latitudes of the Northern Hemisphere, runoff is expected to increase in central and eastern Siberia. In European Russia and in southern Europe, runoff is projected to decrease. In western Siberia (the Ob watershed), runoff would increase under the SRES A1B and A2 scenarios until the 2050s–2070s, then it would decrease to values close to present-day ones; under the anthropogenic forcing scenario SRES B1, the increase in runoff will continue up to the late 21st century. Total runoff from Eurasian rivers into the Arctic Ocean in the IAP RAS CM in the 21st century will increase by 8–9% depending on the scenario. Runoff from the North American rivers into the Arctic Ocean has not changed much throughout numerical experiments with the IAP RAS CM.

Influence of direct sulfate-aerosol radiative forcing on the results of numerical experiments with a climate model of intermediate complexity

The climate model of intermediate complexity developed at the Institute of Atmospheric Physics of the Russian Academy of Sciences (IAP RAS CM) is extended by a block for the direct anthropogenic sulfate-aerosol (SA) radiative forcing. Numerical experiments have been performed with prescribed scenarios of the greenhouse and anthropogenic sulfate radiative forcings from observational estimates for the 19th and 20th centuries and from SRES scenarios A1B, A2, and B1 for the 21st century. The globally averaged direct anthropogenic SA radiative forcing FASA by the end of the 20th century relative to the preindustrial state is −0.34 W/m2, lying within the uncertainty range of the corresponding present-day estimates. The absolute value of FASA is the largest in Europe, North America, and southeastern Asia. A general increase in direct radiative forcing in the numerical experiments that have been performed continues until the mid-21st century. With both the greenhouse and the sulfate loadings included, the global climate warming in the model is 1.5–2.8 K by the end of the 21st century relative to the late 20th century, depending on the scenario, and 2.1–3.4 K relative to the preindustrial period. The sulfate aerosol reduces global warming by 0.1–0.4 K in different periods depending on the scenario. The largest slowdown (>1.5 K) occurs over land at middle and high latitudes in the Northern Hemisphere in the mid-21st century for scenario A2. The IAP RAS CM response to the greenhouse and the aerosol forcing is not additive.

Estimation of changes in characteristics of the climate and carbon cycle in the 21st century accounting for the uncertainty of terrestrial biota parameter values

Abstractensemble simulations with the A.M. Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS) climate model (CM) for the 21st century are analyzed taking into account anthropogenic forcings in accordance with the Special Report on Emission Scenarios (SRES) A2, A1B, and B1, whereas agricultural land areas were assumed to change in accordance with the Land Use Harmonization project scenarios. Different realizations within these ensemble experiments were constructed by varying two governing parameters of the terrestrial carbon cycle. The ensemble simulations were analyzed with the use of Bayesian statistics, which makes it possible to suppress the influence of unrealistic members of these experiments on their results. It is established that, for global values of the main characteristics of the terrestrial carbon cycle, the SRES scenarios used do not differ statistically from each other, so within the framework of the model, the primary productivity of terrestrial vegetation will increase in the 21st century from 74 ± 1 to 102 ± 13 PgC yr−1 and the carbon storage in terrestrial vegetation will increase from 511 ± 8 to 611 ± 8 PgC (here and below, we indicate the mean ± standard deviations). The mutual compensation of changes in the soil carbon stock in different regions will make global changes in the soil carbon storage in the 21st century statistically insignificant. The global CO2 uptake by terrestrial ecosystems will increase in the first half of the 21st century, whereupon it will decrease. The uncertainty interval of this variable in the middle (end) of the 21st century will be from 1.3 to 3.4 PgC yr−1 (from 0.3 to 3.1 PgC yr−1). In most regions, an increase in the net productivity of terrestrial vegetation (especially outside the tropics), the accumulation of carbon in this vegetation, and changes in the amount of soil carbon stock (with the total carbon accumulation in soils of the tropics and subtropics and the regions of both accumulation and loss of soil carbon at higher latitudes) will be robust within the ensemble in the 21st century, as will the CO2 uptake from the atmosphere only by terrestrial ecosystems located at extratropical latitudes of Eurasia, first and foremost by the Siberian taiga. However, substantial differences in anthropogenic emissions between the SRES scenarios in the 21st century lead to statistically significant differences between these scenarios in the carbon dioxide uptake by the ocean, the carbon dioxide content in the atmosphere, and changes in the surface air temperature. In particular, according to the SRES A2 (A1B, B1) scenario, in 2071–2100 the carbon flux from the atmosphere to the ocean will be 10.6 ± 0.6 PgC yr−1 (8.3 ± 0.5, 5.6 ± 0.3 PgC yr−1), and the carbon dioxide concentration in the atmosphere will reach 773 ± 28 ppmv (662 ± 24, 534 ± 16 ppmv) by 2100. The annual mean warming in 2071–2100 relatively to 1961–1990 will be 3.19 ± 0.09 K (2.52 ± 0.08, 1.84 ± 0.06 K).