BESM-OA2.5 is the result of coupling the Center for Weather Forecast and Climate Studies (CPTEC/INPE) Brazilian Atmospheric Model (BAM; ) and the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model version 4p1 via the Flexible Modular System (FMS) (also from GFDL). The dynamical core and physical parameterizations of the atmospheric component of BESM-OA2.5 are the same as those discussed in . BAM is a hydrostatic model, with its dynamical core based on the spectral transform method, which employs global spherical harmonic basis functions. The Eulerian advection scheme option is used in this study, with a two-time-level semi-Lagrangian scheme for the transport of moisture and microphysics prognostic variables, which are carried out completely within the model grid space. Simplified fast physical parameterizations are used here to increase the computational efficiency of long integrations, thus resulting in a decreased computational demand compared with that required by the operational Numerical Weather Prediction (NWP) model. A summary of the main differences in the physical parameterizations between BAM (as used in this paper) and the BAM NWP operational model is provided in Table . The dynamical equations in BAM are discretized following a spectral transform with horizontal resolution truncated at triangular wavenumber 62 (an equivalent grid size of approximately 1.875∘) and 28 layers unevenly spaced in the vertical sigma coordinate with the top level at approximately 2.73 hPa. The oceanic component uses a tripolar grid at a horizontal resolution of 1∘ in longitude, and in the latitudinal direction the grid spacing is 1/4∘ between 10∘ S and 10∘ N, decreasing uniformly to 1∘ at 45∘ and to 2∘ at 90∘ in both hemispheres. The ocean grid has 50 vertical levels with a 10 m resolution in the upper 220 m, decreasing gradually to approximately 370 m at deeper levels.

showed that BESM-OA2.5 can capture the general mean climate state; however, substantial biases appeared in the simulation associated with a double Intertropical Convergence Zone (ITCZ) over the Pacific and Atlantic oceans and regional biases in the precipitation over the Amazon and Indian regions. BESM-OA2.5 can also reproduce the most important large-scale interannual and decadal climate variability patterns. The Atlantic Meridional Mode (AMM) is well simulated by the model in terms of its spatial pattern and temporal variability, whereas this mode is poorly represented in most CMIP5 models . The maximum strength of the Atlantic Meridional Overturning Circulation (AMOC) represented by BESM-OA2.5 is 14 Sv, which is lower than the value determined within the RAPID project (17 Sv; ) but in the range of the observed root mean square variability, and this value is comparable to the ensemble AMOC simulated by the CMIP5 models. Moreover, the spatial structures of both the North Atlantic Oscillation (NAO) and the Pacific Decadal Oscillation (PDO) variability are well captured .

The recently developed BESM-OA2.5 is an advancement of BESM-OA2.3, which was presented by . The main differences between BESM-OA2.5 and the previous version (BESM-OA2.3) lie in the atmospheric model and how some surface layer variables are estimated. The total energy balance at the TOA is better represented in BESM-OA2.5 than in BESM-OA2.3, with a reduced global mean bias of approximately -4 W m-2 compared with -20 W m-2 for the latter. It should be noted that BESM-OA2.5 has a new set of parameterizations, mainly related to an improved representation of microphysical processes. For instance, the precipitation in the previous model was parameterized only in terms of the large-scale condensation. Moreover, BESM-OA2.5 underwent improvements in the representation of wind, humidity, and temperature in the surface layer with the use of the similarity functions method presented by . Based on Monin–Obukhov theory, the wind (u10m), humidity (q2m), and temperature (θ2m) are estimated from the values of the first atmospheric model level and the surface, as described in Eqs. (24), (25), and (26) of . Furthermore, the similarity functions ψm and ψh depend on the stability regimes . For BESM-OA2.5, those regimes are associated with stable (ζ/L>0) and unstable (ζ/L≤0) conditions . These diagnostic variables are important for BESM because they are used in the ocean–atmosphere coupling strategy.

Both versions reproduce the main climate variability, particularly over the Atlantic, such as the AMOC and the AMM, but simulate a weak El Niño–Southern Oscillation (ENSO) interannual variability over the equatorial Pacific . Concerning the general mean present-day climate state, BESM-OA2.5 shows improvements in reproducing the Intertropical Convergence Zone (ITCZ), and it reduces both the global precipitation root mean square error (RMSE) and the sea surface temperature (SST) RMSE compared with those modeled by BESM-OA2.3.

Global simulations that were 1 year long and 6-hourly outputs were performed with BAM configured with surface layer schemes based on and , here called BAM-Arya (the original scheme) and BAM-Jimenez (the new scheme), respectively. The normalized RMSE was computed with respect to the reanalysis NCEP-DOE (National Centers for Environmental Prediction – Department of Energy) version 2 . The normalized RMSE of the wind at 10 m and the temperature and humidity at 2 m for the two surface layer schemes was investigated. Consistent improvements of BAM-Jimenez relative to BAM-Arya were noted in all three variables over the oceanic regions. The normalized RMSE analysis over the continents yielded less consistent results, with an improved BAM-Jimenez representation of both winds and temperature but with an inferior representation of the humidity field (figures not shown).

For this study, climate simulations were performed using BESM-OA2.5 (hereinafter BESM) for the piControl (preindustrial control scenario, run for 300 years with atmospheric CO2 concentration invariant at 274 ppmv) and abrupt4xCO2 (run for 150 years after the abrupt quadrupling of atmospheric CO2 at year 151 of the piControl simulation) scenarios; therefore, both experiments were run in parallel for 150 years. These two scenarios are commonly employed in CMIP5 studies for climate sensitivity assessment . Climate change is evaluated as the difference between the abrupt4xCO2 and piControl experiments. In addition, BESM's results were compared with a selection of 25 CMIP5 models listed in Table . All models, including BESM, were interpolated at 2.5∘ × 2.5∘ longitude–latitude horizontal resolution. All CMIP5 model data are available from the Earth System Grid Federation (ESGF).

Here we estimate the climate feedback with two different methods using either a regression according to or radiative kernels . The Gregory method is more straightforward computationally; however, it returns only a global mean value. Moreover, the ECS can be estimated with this method. On the other hand, it is possible to obtain the seasonal feedback for every latitude–longitude point with the radiative kernel method; furthermore, the feedback can be decomposed into different processes.

Monthly mean climatologies were computed for the last 30 years of the piControl and abrupt4xCO2 runs, and the projected climate response to CO2 increase was evaluated from the difference between these abrupt4xCO2 and piControl monthly mean climatologies. The statistical significance of this difference was calculated based on the Student's t test with a significance level of 90 %. Furthermore, to evaluate how similar two spatial patterns are, we used the spatial inner product calculated as ∑(Ai⋅Bi)/(|A|⋅|B|), where A and B are the 2-D variables and i is the spatial index related to their latitude–longitude coordinates.

Figure shows the linear regressions of ΔR , ΔLW (clear sky), and ΔSW (clear sky) against ΔT‾as from BESM. The linear regressions based on all-sky radiative flux are used to estimate ECS, G, and λ, while the regressions based on clear-sky data are used to obtain ΔCRE (as mentioned in the previous section). BESM features G=8.62 W m-2, λ=-1.45 W m-2 K-1, ΔCRE=-0.13 W m-2 K-1, and ECS=2.96 K.

The parameters G, λ, ΔCRE, and ECS as computed for all models are shown in Table . The climate sensitivities of the 26 CMIP5 coupled models (including BESM-OA2.5) were compiled in an extension of the previous work by , who evaluated 15 CMIP5 coupled models. In the present work, we added the following models: ACCESS1-0, ACCESS1-3, bcc-csm1-1, BESM-OA2.5, BNU-ESM, CCSM4, FGOALS-g2, FGOALS-s2, GISS-E2-H, GISS-E2-R, and inmcm4. In the ECS ranges from 2.07 to 4.74 K for the 15 models analyzed there, which is largely confirmed by our analysis. The small differences can possibly be attributed to the interpolation of the data. G and λ vary from 5.01 to 8.95 W m-2 and from -1.66 to -0.60 W m-2 K-1, respectively. The inter-model spread in G among the models is due to differences in the radiative codes used and the rapid adjustment processes in the troposphere and at the surface . The spread in the ECS is more robustly influenced by λ than by G (Fig. ), as was also suggested by . The correlation coefficient between the ECS and λ is -0.82, which is significant at a 1 % significance level (Fig. b). On the other hand, the correlation between the ECS and G is -0.01, which is not statistically significant (Fig. a). Thus, inter-model variation in the balance of feedbacks explains the dispersion in the ECS better than the initial radiative imbalance triggered by the CO2 increase (related to G). Although BESM yielded one of the highest G values among all of the CMIP5 models, it showed a warming response to CO2 doubling that is well within the range of 3.30±0.76 K as presented by the ensemble models.

The ΔCRE from BESM is -0.13 W m-2 K-1, while the CMIP5 models yield ΔCRE values ranging from -0.50 to 0.70 W m-2 K-1. Unlike the ΔCREa, this term does not consider the masking effects of clouds as estimated via the radiative kernel method (Eq. ). Therefore, the ΔCRE cannot be interpreted to reflect a change in the cloud properties alone.

Figure shows the global mean feedbacks for lapse rate, water vapor, lapse rate plus water vapor, albedo, and cloud (SW, LW, and total) for each CMIP5 model. Both radiative kernels are used to test whether the results are sensitive to the choice of radiative kernel and whether the inter-model deviation is greater than the distribution of the radiatively active constituents (temperature, water vapor, albedo, and cloud) of the base model. It is worth clarifying that positive (negative) values of feedbacks contribute to the amplification (damping) of global warming. The strongest positive feedback (Fig. ) is due to water vapor (mean value: 1.39 W m-2 K-1), followed by clouds (mean value: 0.96 W m-2 K-1) and then surface albedo (mean value: 0.32 W m-2 K-1). The global mean Planck feedback is negative, with an average value of -3.60 W m-2 K-1 (not shown in Fig. ), followed by a lapse rate feedback of -0.77 W m-2 K-1. BESM yields values near the ensemble mean for the albedo and cloud feedbacks, i.e., 0.27 and 0.95 W m-2 K-1, respectively. For the lapse rate feedback BESM yields a value of -0.71 W m-2 K-1, slightly underestimating the ensemble mean value in magnitude. In turn, BESM is among the models with the lowest global water vapor feedback average, with a value of around 1.24 W m-2 K-1.

Figure shows the latitudinal profiles for the global mean feedback values in Fig. , allowing us to identify the regions that induce deviations of BESM results from the CMIP ensemble. In Fig. a–b, there is a nearly constant Planck feedback of approximately -3.4 W m-2 K-1 from 90 ∘S to 60∘N, with a notable increase in the ensemble standard deviation in the sub-Antarctic latitudes (around 60 ∘S), which is in accordance with . The exception is in the Arctic region, where the mean value reaches -10 W m-2 K-1 with a similarly increased standard deviation value. In the sub-Antarctic and Arctic latitudes, BESM yields one of the most negative values for the Planck feedback, revealing that BESM has a stronger vertically homogeneous warming (corresponding to large surface warming) among the CMIP5 models. Furthermore, for those same regions, BESM showed more positive lapse rate feedback than the ensemble. Therefore, BESM yields a warming relatively larger at the surface and relatively weaker at the upper troposphere, resulting in a stronger vertical temperature gradient in comparison to the other models.

In the tropics, where there is an intense moist convection, atmospheric warming almost follows a moist adiabat (temperature increase is larger in the upper troposphere compared to that at the surface), implying a negative lapse rate feedback (Fig. c–d) . In accordance with this upper tropospheric warming in the tropics, an increase in the specific humidity occurs (Manabe and Wetherald, 1975), which causes a reinforcement of the greenhouse effect, reflected by a positive water vapor feedback as shown in Fig. d–e. Because of this close link between the lapse rate and water vapor feedbacks, it is common to consider their effects as a sum, as displayed in Fig. . BESM shows a lapse rate feedback near the ensemble mean for the tropics. The greatest BESM deviations are observed near the Antarctic and over the Arctic (Fig. c–d), where this feedback became positive for all models. For the water vapor feedback, greater dispersion of the models was observed in the tropics, with BESM systematically yielding values below the ensemble mean for the same latitude band (Fig. e–f). These lower values extend throughout the Northern Hemisphere, consistent with the low global mean water vapor feedback value relative to the ensemble (shown in Fig. ).

The albedo feedback profiles from BESM and the CMIP5 models are compared in Fig. g–h. Nonzero results mostly occur over the polar regions, where there is a reduction in sea ice and snow cover . The positive albedo feedback signals yielded by all of the models in such regions imply that the reduction in the albedo corresponds to an increase in the radiation budget at the TOA due to reduced upward shortwave radiation. As expected, the polar regions present a large dispersion among the models, which is related to how fast the sea ice melts in the different climate models. The regions over the Arctic and the ocean near the Antarctic show the largest surface warming, and this positive albedo, together with the positive lapse rate feedbacks, is the main factor responsible for a phenomenon known as polar amplification . BESM yielded an albedo feedback greater than the ensemble standard deviation over the Southern Ocean at around 60∘ S. This same latitude is where BESM shows negative Planck and positive lapse rate feedbacks outside the models limits, as previously discussed.

Finally, regarding cloud feedback, most of the inter-model spread arises from the SW component (Figs.  and i–j). This dispersion is also reflected in the standard deviation and in the limit between the minimum and maximum of the zonally averaged cloud feedback shown in Fig. i–j. The SW cloud feedback ranges from -0.28 to 1.40 W m-2 K-1, while the LW cloud effect ranges from 0.10 to 0.96 W m-2 K-1. The combined SW and LW cloud effects result in positive cloud feedback ranging from 0.35 to 1.69 W m-2 K-1. This result is similar to that found by for CMIP3 (IPCC AR4, ) models, who also reported a near-zero to positive cloud feedback. BESM presents positive values of around 0.5 W m-2 K-1 for both SW and LW cloud feedback, resulting in a total cloud feedback of about 1.0 W m-2 K-1 (as shown in Fig. ). The highest positive values are in regions with strong albedo feedback (Fig. i–j).

Although BESM yielded global area-averaged feedbacks near the model ensemble mean values, differences are found mainly at high latitudes. In fact, for cloud feedback, BESM is an outlier due to a strong shortwave component response over both the Arctic and the Southern Ocean near the Antarctic. This effect is evident via the decomposition of the cloud feedback into the SW and LW components for the ensemble and BESM values, as shown in Fig. . To assess the analytical causes of this strong BESM shortwave cloud feedback departure from the ensemble, we separately computed the contributors to the SW cloud feedback, i.e., the SW CRE (as described by ) and the feedback cloud masks, as in Eq. (). For the shortwave component, the feedback cloud masks are obtained for the albedo and SW water vapor feedbacks by calculating the all-sky minus clear-sky radiation flux for each feedback. We find that the higher BESM cloud feedback values (Fig. g–h) are mainly consequences of the sum of the SW CRE (Fig. a–b) and the effect of cloud masking for the albedo feedback (Fig. c–d) in the sub-Antarctic and Arctic regions. In turn, the cloud masking for the SW water vapor (Fig. e–f) does not contribute to the higher positive BESM values. As shown in Fig. , it is possible to attribute BESM's status as an outlier over the Arctic region to the SW CRE, while in the Southern Ocean (around 60∘ S), the major contribution comes from the albedo feedback cloud mask.

A deepened physical analysis to understand the BESM cloud feedback behavior in high latitudes is obtained by examining the zonal mean of the change in the cloud vertical profile for BESM, as shown in Fig. . Over the Arctic and near the Antarctic, BESM showed an increase in the cloud fraction above 850 hPa and a decrease below that level, indicating an upward shift of low-level clouds (Fig. a). However, the increase in cloud cover above 850 hPa is stronger than the reduction below. Because of this increase in the total cloud fraction, a negative SW CRE (not adjusted) appears in those regions (Fig. a–b), consistent with an increase in sun shading (Fig. b). Moreover, the SW cooling is smaller than the heating provided by LW radiation due to the upward shift of the low-level clouds, as evident in the net effect (Fig. d). Furthermore, the positive albedo and lapse rate feedbacks (Fig. c–f) are consistent with this vertical cloud shifting. In this manner, a loss of SW energy at the surface associated with an increase in the total cloud fraction explains the negative SW CRE, and the gain of LW energy is responsible for the sea ice melting. Consequently, this gain of LW energy is indirectly linked to the albedo feedback cloud mask for BESM, since the mask (Δa/ΔT‾as(Ka-Kacs)) is proportional to the albedo change (Δa). As discussed before, both the SW CRE and albedo feedback cloud mask contribute to the large positive cloud feedback over the Arctic and sub-Antarctic areas observed in BESM.

Figure shows the annual mean surface temperature differences between the abrupt4xCO2 and piControl scenarios for the ensemble of 25 CMIP5 models and BESM. As clearly shown in Fig. , despite the generalized increase in the surface temperature over most of the globe, BESM shows a generally lower temperature increase, principally over the continental areas. The CMIP5 ensemble yields a mean continental temperature increase of 6.78 K, while BESM yields a value of 5.57 K. Nevertheless, the spatial patterns of the temperature increases are similar, as measured by the spatial inner product (as described in the previous section) between Fig. a and b, which results in a value of 0.96 (values near 1 indicate that both variables have a similar spatial pattern, whereas values near 0 mean that there is hardly any pattern correlation between variables). One point of interest within the scientific community is the relatively low temperature increase over the subpolar North Atlantic, also referred to as the warming hole . In the CMIP5 ensemble mean, the North Atlantic does not show a temperature decrease, although it is the region with the smallest temperature increase globally, while BESM shows an area of temperature decrease in this region. Such a decrease is also present in six other analyzed models (CSIRO-Mk3-6-0, FGOALS-s2, GFDL-ESM2G, GFDL-ESM2M, GISS-E2-R, and inmcm4). These results are consistent with those of , who showed that both observations and CMIP5 models present maximum cooling in the center of the subpolar gyre. The authors argue that there is evidence that both the subpolar gyre and the AMOC adjust in concert with different time lags.

The regions with the largest temperature increase in the abrupt4xCO2 scenario are the polar regions, particularly over the Arctic. The equatorial Pacific shows a relative maximum in warming when the abrupt4xCO2 scenario is compared with the piControl both in the CMIP5 ensemble and BESM. Such changes in the Pacific mean state are consistent with the IPCC AR5, which reports that the Pacific Ocean becomes warmer near the Equator compared to the subtropics in the CMIP5 projections . The scatter plot of the global average under abrupt4xCO2 conditions versus the piControl conditions presented in Fig.  provides additional information that helps to understand the dispersion around the mean value among the different models. Even though the outputs of most of the models lie in either quadrant 1 or 3 (top right and bottom left, respectively), it is not possible to claim any robust linear relationship. This result indicates that models with warmer or cooler mean climates in the piControl runs apparently do not show a corresponding warmer or cooler climate in the abrupt4xCO2 experiments. BESM yields a temperature near that of the ensemble in both the piControl and abrupt4xCO2 runs; consequently, it also showed a temperature increase near the ensemble mean, consistent with its Planck feedback (Figs.  and a–b).

Figure shows the precipitation changes between the abrupt4xCO2 and piControl scenarios for the multi-model ensemble and the BESM. The results are, in general, similar to those of , with wet regions becoming wetter (near-equatorial and subpolar regions) and dry regions becoming drier (centered around 30∘ in both hemispheres). The precipitation pattern in the CMIP5 ensemble shows increased precipitation over the equatorial Pacific, which could be related to the equatorial Pacific warming pattern shown in the temperature change (Fig. ). Furthermore, the CMIP5 ensemble shows a decrease in precipitation in northern South America. The BESM precipitation pattern is similar to the spatial patterns in the CMIP5 ensemble but with some notable exceptions. For example, the decreased precipitation over the South Pacific shown in the CMIP5 ensemble plot is extended into the Indonesian region in BESM. It is also worth noting that in the BESM simulation, the South Pacific convergence zone (SPCZ) shifts southward in the abrupt4xCO2 scenario compared with its position in the piControl scenario. In both the multi-model ensemble and BESM, the precipitation change pattern over South America is similar to that which occurs during El Niño years , with increased precipitation over southeastern South America and decreased precipitation over northern–northeastern South America. The scatter plot in Fig.  emphasizes a linear relationship between the experiments, indicating that models with higher (lower) global average precipitation in the piControl scenario show higher (lower) precipitation in the abrupt4xCO2 scenario. As is obvious from Fig. , the BESM performance perfectly matches the ensemble mean behavior in global mean.

Figure shows a scatter plot of the ECS versus the precipitation difference between the abrupt4xCO2 and piControl scenarios (ΔPr) for all of the considered models. It is worth noting that all the models show increased global mean precipitation upon a quadrupling of atmospheric CO2 with piControl preindustrial CO2 concentrations (positive values on the y axis in Fig. ). An apparent linear relationship between these differences (abrupt4xCO2 minus piControl) in the global mean precipitation and ECS is also evident in Fig. , in which the warmest models tend to have the largest precipitation changes. The slope of the linear regression reflects a 2.5 % precipitation change per Kelvin, which is close to that found by . This slope is much lower than that predicted by the Clausius–Clapeyron relation, i.e., an approximately 6.5 % change in precipitation per Kelvin. In fact, such precipitation increases are not governed by moisture availability but rather by the surface and tropospheric energy balance, which incorporates the surface radiative heating, surface latent heat flux, and radiative cooling of the troposphere .

MRI-CGCM3, ACCESS1-0, and HadGEM2-ES show greater deviations from the linear fit shown in Fig. . Furthermore, BESM is marginally out of the residual standard error interval, with a 9.5 % increase in precipitation (the error limit is 9.2 %). ACCESS1-0 and HadGEM2-ES use the same atmospheric model , which could explain the lower increase in precipitation in both coupled models.

In addition to temperature and precipitation changes, we are also interested in understanding the changes in the BESM atmospheric circulation (compared to other models) following a quadrupling of the CO2 concentration. The sea level pressure (SLP) response patterns shown in Fig.  depict a poleward shift in the subtropical high-pressure cells in both the CMIP5 ensemble and BESM. Furthermore, when the models are subjected to the increased atmospheric CO2 concentration, a decrease in the SLP over the polar regions is evident. This, connected with the increase in the midlatitudes, indicates a positive trend in Arctic oscillation (AO) and Antarctic oscillation (AAO) episodes, which have already been reported by , , and . It is also interesting to note the statistically significant SLP decrease (increase) over the eastern (western) Pacific, a pattern that might indicate an ENSO-like pattern in scenarios with an increased CO2 concentration. This pattern is consistent with those depicted in Fig.  for the SST changes in a 4xCO2 scenario.

The results from the piControl scenario (the contours in Fig. ) show that the Southern Hemisphere subtropical jet, reflected by the core of the maximum eastward zonal wind, is localized around 35∘ S at 200–150 hPa in both the CMIP5 ensemble and BESM. In Fig. a and b (BESM and the CMIP5 ensemble), we note that the regions with the strongest increases in westerly winds at all levels show a southward jet displacement. This observation is consistent with the poleward displacement of the high SLP center shown in Fig. . Furthermore, as the high-pressure centers experienced a poleward shift, the pressure gradients intensified in the subpolar areas; consequently, the near-surface wind velocity increased, following the geostrophic approximation (u≈-(1/fρ)(∂p/∂y)), where f is the Coriolis parameter and ρ is the air density.

Figure shows the average differences from 5∘ N to 5∘ S (Walker circulation) between the abrupt4xCO2 and piControl scenarios for vertical motion (omega; shaded) and zonal wind and vertical motion (vectors). According to the vertical motion pattern in the piControl (contours), the multi-model ensemble and BESM show subsidence over an extensive area in the Pacific (150∘ E–90∘ W) whose intensity is lower in the abrupt4xCO2 simulation, as shown in Fig.  (blue). This finding is coherent with near-surface temperature patterns (Fig. ), which show an equatorial warming pattern in the mean state (e.g., during El Niño years, a weakening of the Walker circulation occurs). Furthermore, there is enhanced subsidence for the difference between the two scenarios over South America (around 75∘ W), consistent with the decrease in precipitation in tropical South America, in both BESM and the CMIP5 ensemble (Fig. ).