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 CO₂. 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 CH₄/year for the preindustrial and current period increase to 170–200 MtCH₄/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/m² and 0.05 K, respectively).     

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 F _ASA by the end of the 20th century relative to the preindustrial state is ?0.34 W/m², lying within the uncertainty range of the corresponding present-day estimates. The absolute value of F _ASA 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.     

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 S _p in the IAP RAS CM changes only slightly until the late 20th century, reaching about 21 million km², and then decreases to 11–12 million km² in 2036–2065 and 4–8 million km² 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 km² 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 km² to 31–37 million km². 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.     

Effect of including land-use driven radiative forcing of the surface albedo of land on climate response in the 16th–21st centuries

A change in ecosystem types, such as through natural-vegetation-agriculture conversion, alters the surface albedo and triggers attendant shortwave radiative forcing (RF). This paper describes numerical experiments performed using the climate model (CM) of the Institute of Atmospheric Physics (IAP), Russian Academy of Sciences, for the 16th–21st centuries; this model simulated the response to a change in the contents of greenhouse gases (tropospheric and stratospheric), sulfate aerosols, solar constant, as well as the response to change in surface albedo of land due to natural-vegetation-agriculture conversion. These forcing estimates relied on actual data until the late 20th century. In the 21st century, the agricultural area was specified according to scenarios of the Land Use Harmonization project and other anthropogenic impacts were specified using SRES scenarios. The change in the surface vegetation during conversion from natural vegetation to agriculture triggers a cooling RF in most regions except for those of natural semiarid vegetation. The global and annual average RF derived from the IAP RAS CM in late 20th century is ?0.11 W m^?2. Including the land-use driven RF in IAP RAS CM appreciably reconciled the model calculations to observations in this historical period. For instance, in addition to the net climate warming, IAP RAS CM predicted an annually average cooling and reduction in precipitation in the subtropics of Eurasia and North America and in Amazonia and central Africa, as well as a local maximum in annually average and summertime warming in East China. The land-use driven RF alters the sign in the dependence that the amplitude of the annual cycle of the near-surface atmospheric temperature has on the annually averaged temperature. One reason for the decrease in precipitation as a result of a change in albedo due to land use may be the suppression of the convective activity in the atmosphere in the warm period (throughout the year in the tropics) and the corresponding decrease in convective precipitation. In the 21st century, the effect that the land-use driven RF has on the climate response for scenarios of anthropogenic impact is generally small.     

Validating and assessing the sensitivity of the climate model with an ocean general circulation model developed at the Institute of Atmospheric Physics,Russian Academy of Sciences

A new version of the Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS), climate model (CM) has been developed using an ocean general circulation model instead of the statistical-dynamical ocean model applied in the previous version. The spatial resolution of the new ocean model is 3° in latitude and 5° in longitude, with 25 unevenly spaced vertical levels. In the previous version of the oceanic model, as in the atmospheric model, the horizontal resolution was 4.5° in latitude and 6° in longitude, with four vertical levels (the upper quasi-homogeneous layer, seasonal thermocline, abyssal ocean, and bottom friction layer). There is no correction for the heat and momentum fluxes between the atmosphere and ocean in the new version of the IAP RAS CM. Numerical experiments with the IAP RAS CM have been performed under current initial and boundary conditions, as well as with an increasing concentration of atmospheric carbon dioxide. The main simulated atmospheric and oceanic fields agree quite well with observational data. The new version’s equilibrium temperature sensitivity to atmospheric CO₂ doubling was found to be 2.9 K. This value lies in the mid-range of estimates (2–4.5 K) obtained from simulations with state-of-the-art models of different complexities.     

Estimating the efficiency of mitigating and preventing global warming with scenarios of controlled aerosol emissions into the stratosphere

Ensemble numerical experiments with the climate model of intermediate complexity developed at the A.M. Obukhov Institute of Atmospheric Physics of the Russian Academy of Sciences (IAP RAS CM) are conducted to estimate the efficiency of controlled climate forcing (geoengineering) due to stratospheric sulfate aerosol (SSA) emissions in order to compensate for global warming under the SRES A1B anthropogenic emission scenario. Full (or even excessive) compensation for the expected anthropogenic warming in the model is possible with sufficiently intense geoengineering. For ensemble members with values of the governing parameters corresponding to those obtained for the Mt. Pinatubo eruption, global warming is reduced by no more than 0.46 K in the second half of the 21st century, with a residual rise in the global surface temperature T _g comparative to 1961–1990 of 1.0–1.2 K by 2050 and 1.9–2.2 K by 2100. The largest reduction in global warming (with the other parameters of the numerical experiment being equal) is found not for a meridional distribution of SSA concentration peaked at low latitudes (despite the largest (in magnitude) global compensation instantaneous radiative forcing), but for a uniform horizontal aerosol distribution and for a distribution with the SSA concentration maximum in the middle and subpolar latitudes of the Northern Hemisphere. The efficiency of geoengineering in terms of T _g in the second half of the 21st century between the most efficient and the least efficient meridional distributions of stratospheric aerosols differs by as much as one-third, depending on the values of other governing parameters. For meridional distributions of SSA concentration, which produce the largest deceleration of global warming, such a deceleration is regionally most pronounced over high- and subpolarlatitude land areas and in the Arctic. In particular, this is expressed in the smallest reduction in the sea-ice extent and permafrost area under climate warming in the model. The compensation forcing also decelerates a general increase in global annual precipitation P _g during warming. The relative deceleration in precipitation increase is most pronounced in land regions outside the tropics, where a significant deficit in precipitation is currently observed. After the theoretical completion of geoengineering in the first or second decade, its temperature effect vanishes with an abrupt acceleration of global and regional surface warming. For individual members of the ensemble experiment, the global temperature change in this period is five times as large as that in the experiment without geoengineering and ten times as large regionally (in northeastern Siberia).     

The influence of lightning activity and anthropogenic factors on large-scale characteristics of natural fires

A module for simulating of natural fires (NFs) in the climate model of the A.M. Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS CM), is extended with respect to the influence of lightning activity and population density on the ignition frequency and fire suppression. The IAP RAS CM is used to perform numerical experiments in accordance with the conditions of the project that intercompares climate models, CMIP5 (Coupled Models Intercomparison Project, phase 5). The frequency of lightning flashes was assigned in accordance with the LIS/OTD satellite data. In the calculations performed, anthropogenic ignitions play an important role in NF occurrences, except for regions at subpolar latitudes and, to a lesser degree, tropical and subtropical regions. Taking into account the dependence of fire frequency on lightning activity and population density intensifies the influence of characteristics of natural fires on the climate changes in tropics and subtropics as compared to the version of the IAP RAS CM that does not take the influence of ignition sources on the large-scale characteristics of NFs into consideration.     

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.     

Climate and carbon cycle variations in the 20th and 21st centuries in a model of intermediate complexity

The climate model of intermediate complexity developed at the Oboukhov Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS CM), has been supplemented by a zero-dimensional carbon cycle model. With the carbon dioxide emissions prescribed for the second half of the 19th century and for the 20th century, the model satisfactorily reproduces characteristics of the carbon cycle over this period. However, with continued anthropogenic CO₂ emissions (SRES scenarios A1B, A2, B1, and B2), the climate-carbon cycle feedback in the model leads to an additional atmospheric CO₂ increase (in comparison with the case where the influence of climate changes on the carbon exchange between the atmosphere and the underlying surface is disregarded). This additional increase is varied in the range 67–90 ppmv depending on the scenario and is mainly due to the dynamics of soil carbon storage. The climate-carbon cycle feedback parameter varies nonmonotonically with time. Positions of its extremes separate characteristic periods of the change in the intensity of anthropogenic emissions and of climate variations. By the end of the 21st century, depending on the emission scenario, the carbon dioxide concentration is expected to increase to 615–875 ppmv and the global temperature will rise by 2.4–3.4 K relative to the preindustrial value. In the 20th–21st centuries, a general growth of the buildup of carbon dioxide in the atmosphere and ocean and its reduction in terrestrial ecosystems can be expected. In general, by the end of the 21st century, the more aggressive emission scenarios are characterized by a smaller climate-carbon cycle feedback parameter, a lower sensitivity of climate to a single increase in the atmospheric concentration of carbon dioxide, a larger fraction of anthropogenic emissions stored in the atmosphere and the ocean, and a smaller fraction of emissions in terrestrial ecosystems.     

Sensitivity of amplitude-phase characteristics of the surface air temperature annual cycle to variations in annual mean temperature

The ERA40 and NCEP/NCAR data over 1958–1998 were used to estimate the sensitivity of amplitude-phase characteristics (APCs) of the annual cycle (AC) of the surface air temperature (SAT) T _s. The results were compared with outputs of the ECHAM4/OPYC3, HadCM3, and INM RAS general circulation models and the IAP RAS climate model of intermediate complexity, which were run with variations in greenhouse gases and sulfate aerosol specified over 1860–2100. The analysis was performed in terms of the linear regression coefficients b of SAT AC APCs on the local annual mean temperature and in terms of the sensitivity characteristic D = br ², which takes into account not only the linear regression coefficient but also its statistical significance (via the correlation coefficient r). The reanalysis data were used to reveal the features of the tendencies of change in the SAT AC APCs in various regions, including areas near the snow-ice boundary, storm-track ocean regions, large desert areas, and the tropical Pacific. These results agree with earlier observations. The model computations are in fairly good agreement with the reanalysis data in regions of statistically significant variations in SAT AC APCs. The differences between individual models and the reanalysis data can be explained, in particular, in terms of the features of the sea-ice schemes used in the models. Over the land in the middle and high latitudes of the Northern Hemisphere, the absolute values of D for the fall phase time and the interval of exceeding exhibit a positive intermodel correlation with the absolute value of D for the annual-harmonic amplitude. Over the ocean, the models reproducing larger (in modulus) sensitivity parameters of the SAT annual-harmonic amplitude are generally characterized by larger (in modulus) negative sensitivity values of the semiannual-harmonic amplitude T _{s, 2}, especially at latitudes characteristic of the sea-ice boundary. In contrast to the averaged fields of AC APCs and their interannual standard deviations, the sensitivity parameters of the SAT AC APCs on a regional scale vary noticeably for various types of anthropogenic forcing.     

Comparison of climatic efficiency of the mechanisms of land-surface albedo changes caused by land use

Changes in ecosystem types, including situations when natural vegetation is replaced by agricultural lands, leads to surface albedo changes and the development of the corresponding short-wave radiative forcing (RF). This work analyzes ensemble numerical experiments with the climatic model (CM) of the Institute of Atmospheric Physics at the Russian Academy of Sciences (IAP RAS) for the 16th–21st centuries. The responses to changes in the contents of greenhouse gases and sulfate aerosols (tropospheric and stratospheric), in the solar constant, and in the land-surface albedo when natural vegetation is replaced by agricultural lands were modeled during these experiments. Different members of these ensemble experiments were obtained by varying the model parameters affecting the RF on the climate during land use: the albedo of agricultural lands was varied within the interval from 0.15 to 0.25 and the parameter controlling the efficiency of snow masking by tree vegetation was varied in the range from the absence of this effect to its maximally possible efficiency. It has been established that changes in surface albedo when natural vegetation is replaced by agricultural lands have the largest influence on the globally averaged annual mean RF at the top of the atmosphere, whereas the influence of snow masking on the RF is substantially less. This phenomenon is caused by the fact that snow masking by tree vegetation can take place only in winter in regions of temperate and high latitudes, when insolation is relatively low. A comparison of the spatial structure of the annual mean response of the surface temperature with the HadCRUT3v and GISS observational data makes it possible to narrow the admissible range of model parameter values. In particular, it can be inferred that the key parameter values which control the influence that land use has on the surface albedo in the IAP RAS CM are close to optimal. In addition, variations in these parameters do not lead to a significant influence of land use on climate change in the 21st century if the Land Use Harmonization (LUH) scenarios of changes in the area of agricultural lands are used: the uncertainty of the model response associated with the uncertainty of values of such controlling parameters in the 21st century does not exceed 0.1 K.     

Estimation of the uncertainty of future changes in atmospheric carbon dioxide concentration and its radiative forcing

An ensemble experiment with the IAP RAS CM was performed to estimate future changes in the atmospheric concentration of carbon dioxide, its radiative forcing, and characteristics of the climate-carbon cycle feedback. Different ensemble members were obtained by varying the governing parameters of the terrestrial carbon cycle of the model. For 1860–2100, anthropogenic CO₂ emissions due to fossil-fuel burning and land use were prescribed from observational estimates for the 19th and 20th centuries. For the 21st century, emissions were taken from the SRES A2 scenario. The ensemble of numerical experiments was analyzed via Bayesian statistics, which made the uncertainty range of estimates much narrower. To distinguish between realistic and unrealistic ensemble members, the observational characteristics of the carbon cycle for the 20th century were used as a criterion. For the given emission scenario, the carbon dioxide concentration expected by the end of the 21st century falls into the range 818 ± 46 ppm (an average plus or minus standard deviation). The corresponding global instantaneous radiative forcing at the top of the atmosphere (relative to the preindustrial state) lies in the uncertainty range 6.8 ± 0.4 W m^?2. The uncertainty range of the strength of the climate-carbon cycle feedback by the end of the 21st century reaches 59 ± 98 ppm in terms of the atmospheric carbon dioxide concentration and 0.4 ± 0.7 W m^?2 in terms of the radiative forcing.     

Climate of the last millennium: a sensitivity study

Seventy-one sensitivity experiments have been performed using a two-dimensional sector-averaged global climate model to assess the potential impact of six different factors on the last millennium climate and in particular on the surface air temperature evolution. Both natural (i.e, solar and volcanism) and anthropogenically-induced (i.e. deforestation, additional greenhouse gases, and tropospheric aerosol burden) climate forcings have been considered.
Comparisons of climate reconstructions with model results indicate that all the investigated forcings are needed to simulate the surface air temperature evolution. Due to uncertainties in historical climate forcings and temperature reconstructions, the relative importance of a particular forcing in the explanation of the recorded temperature variance is largely function of the forcing time series used. Nevertheless, our results indicate that whatever the historical solar and volcanic reconstructions may be, these externally driven natural climate forcings are unable to give climate responses comparable in magnitude and time to the late–20th-century temperature warming while for earlier periods combination of solar and volcanic forcings can explain the Little Ice Age and the Medieval Warm Period. Only the greenhouse gas forcing allows the model to simulate an accelerated warming rate during the last three decades. The best guess simulation (largest similarity with the reconstruction) for the period starting 1850 AD requires however to include anthropogenic sulphate forcing as well as the impact of deforestation to constrain the magnitude of the greenhouse gas twentieth century warming to better fit the observation. On the contrary, prior to 1850 AD mid-latitude land clearance tends to reinforce the Little Ice age in our simulations.     

20世纪太平洋海温变化中人为因子与自然因子贡献的模拟研究

基于中国科学院大气物理研究所大气科学和地球流体力学国家重点实验室(LASG/IAP)发展的气候系统模式FGOALS_gl对20世纪太平洋海温变化的模拟,讨论了自然因子和人为因子对20世纪太平洋海温变化的相对贡献。观测资料表明,20世纪太平洋平均的SST变化主要分为3个时段:20世纪上半叶的增暖,40—70年代的微弱变冷,70年代之后的迅速增暖。20世纪太平洋SST变化的主导模态是全海盆尺度的振荡上升模态,其次为PDO振荡型,在70年代末PDO存在明显的年代际转型。通过全强迫试验、自然强迫试验、控制试验对上述现象进行归因分析,结果表明,人为因子和内部变率都对第一次增暖有贡献,而人类活动(主要是温室气体的增加)是70年代之后太平洋SST迅速增暖的主要原因。分区域来看,在两个增暖时段中,影响黑潮延伸体区SST变化的主要是自然因子和内部变率,影响其它海域SST变化的则主要是人为因子。全强迫试验可以较好的模拟出前两个模态的空间分布及时间序列。在没有人为因子的影响下,PDO成为太平洋海温变化的主导模态,其年代际转变发生在60年代中期,意味着人为因子是全海盆振荡增暖的主导原因,并且它使得年代际转型滞后了10a。因此,自然因子是导致SST年代际转型中的主导因子,人为因子有"调谐"作用。     

Changes of cooling near mesopause under global warming from observations and model simulations

The results of joint analysis of temperature variations near mesopause from long-term measurements at the Zvenigorod Scientific Station of the Obukhov Institute of Atmospheric Physics RAS in 1960–2015 and variations of surface air temperature characterizing global climate change. Together with variations of temperature at the mesopause T _ms from measurements of the hydroxyl emissions we analyzed the temperature variations near mesopause T _m reduced to the same level of solar activity. The observed strong decrease in temperature near mesopause during last decades, particularly in winter, with its tendency to slow down since the 1980’s is was detected against the background of general increase in the surface air temperature of the Northern Hemisphere T _NHs and the Earth as a whole. It was revealed a sharp drop in winter temperature near mesopause in 1970s. and its synchronicity with the shift in climatic features at the surface associated with changes in formation of El Nino events and their impact on the global climate. The general significant negative correlation of temperature variations near mesopause and T _NHs detected from 56-year observational data was not accompanied by any significant coherence between the most long-period temperature variations from the cross-wavelet analysis. To assess the possible manifestation of this coherence the results of numerical simulations with a global climate model were used. According to model simulations for the 20–21 centuries taking into account anthropogenic forcings for significant coherence between long-term variations T _m and T _NHs the prolonged observations are required for temperature near mesopause–about a century or more.     

Estimates of changes in the rate of methane sink from the atmosphere under climate warming

The temperature dependence of the methane oxidation rate is estimated. The methane lifetime in the atmosphere is shown to decrease by about 3% from 1900 to 2005. The overwhelming fraction of the total methane content is removed from the atmosphere at intratropical latitudes during the daytime. The methane oxidation rate growth due to the temperature increase in the troposphere generates negative feedback in the methane cycle and, accordingly, climatic feedback with the same sign. According to the estimates performed, the halt in methane concentration growth in the atmosphere observed in recent years can be associated with a decrease in the lifetime of methane in the atmosphere. According to the results of numerical experiments with the climatic model of the Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS), the climatic effect of negative feedback of the tropospheric temperature and the methane lifetime in the atmosphere is not large and is comparable with the climatic forcing of the methane emission growth from bog ecosystems.     

Methane cycle in the INM RAS climate model

The atmosphere-ocean general circulation model with the carbon cycle is coupled to a model of methane evolution, in which methane sources in the soil of wetlands and methane evolution in the atmosphere are calculated. A numerical experiment on the simulation of climate and methane-cycle changes in 1860–2100 has been conducted with the model forced by methane emissions prescribed from scenario A1B. The distribution of the sources of methane from soil agrees with the available estimates and amounts to about 240 Mt/year in the 20th century. The methane flux from soil increases to 340 Mt/year by the end of the 21st century. The model adequately reproduces an increase in the atmospheric methane concentration from 800 ppb in 1860 to about 1800 ppb in 2000, but does not produce the observed stabilization of methane concentration in the early 21st century. By 2060, the methane concentration in the model attains 2700 ppb. The increase in atmospheric methane concentration is due mainly to anthropogenic emissions. A similar numerical experiment with fixed sources of methane from soil at the 1860–1900 level suggests that the maximum methane concentration in the model in this case could amount to 2400 ppb. A temperature increase at the end of the 21st century relative to the 19th century is 3.5° for a simulated change in the methane flux from soil and 0.25° less for a fixed methane flux.     

Decadal variability of the shallow Pacific meridional overturning circulation: Relation to tropical sea surface temperatures in observations and climate change models

In this study, we use existing observational datasets to evaluate 20th century climate simulations of the tropical Pacific. The emphasis of our work is decadal variability of the shallow meridional overturning circulation, which links the tropical and subtropical Pacific Ocean. In observations, this circulation is characterized by equatorward geostrophic volume transport convergence in the interior ocean pycnocline across 9°N and 9°S. Historical hydrographic data indicate that there has been a decreasing trend in this convergence over the period 1953–2001 of about 11 Sverdrup (1 Sv = 10⁶ m³ s⁻¹), with maximum decade-to-decade variations of 7–11 Sv. The transport time series is highly anti-correlated with sea surface temperature (SST) anomalies in the central and eastern tropical Pacific, implying that variations in meridional overturning circulation are directly linked to decadal variability and trends in tropical SST. These relationships are explored in 18 model simulations of 20th century climate from 14 state-of-the-art coupled climate models. Significant correlation exists between meridional volume transport convergence and tropical SST in the majority of the models over the last half century. However, the magnitude of transport variability on decadal time scales in the models is underestimated while at the same time modeled SST variations are more sensitive to that transport variability than in the observations. The effects of the meridional overturning circulation on SST trends in most the models is less clear. Most models show no trend in meridional transport convergence and underestimate the trend in eastern tropical Pacific SST. The eddy permitting MIROCH model is the only model that reasonably reproduces the observed trends in transport convergence, tropical Pacific SST, and SST gradient along the equator over the last half century. If the observed trends and those simulated in the MIROCH model are ultimately related to greenhouse gas forcing, these results suggest that the Bjerknes feedback, by affecting pycnocline transport convergences, may enhance warming that arises from anthropogenic forcing in the eastern tropical Pacific.     

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

ensemble 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⁻¹ 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 CO₂ 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⁻¹ (from 0.3 to 3.1 PgC yr⁻¹). 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 CO₂ 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⁻¹ (8.3 ± 0.5, 5.6 ± 0.3 PgC yr⁻¹), 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).