Climate feedbacks determined using radiative kernels in a multi-thousand member ensemble of AOGCMs

The use of radiative kernels to diagnose climate feedbacks is a recent development that may be applied to existing climate change simulations. We apply the radiative kernel technique to transient simulations from a multi-thousand member perturbed physics ensemble of coupled atmosphere-ocean general circulation models, comparing distributions of model feedbacks with those taken from the CMIP-3 multi GCM ensemble. Although the range of clear sky longwave feedbacks in the perturbed physics ensemble is similar to that seen in the multi-GCM ensemble, the kernel technique underestimates the net clear-sky feedbacks (or the radiative forcing) in some perturbed models with significantly altered humidity distributions. In addition, the compensating relationship between global mean atmospheric lapse rate feedback and water vapor feedback is found to hold in the perturbed physics ensemble, but large differences in relative humidity distributions in the ensemble prevent the compensation from holding at a regional scale. Both ensembles show a similar range of response of global mean net cloud feedback, but the mean of the perturbed physics ensemble is shifted towards more positive values such that none of the perturbed models exhibit a net negative cloud feedback. The perturbed physics ensemble contains fewer models with strong negative shortwave cloud feedbacks and has stronger compensating positive longwave feedbacks. A principal component analysis used to identify dominant modes of feedback variation reveals that the perturbed physics ensemble produces very different modes of climate response to the multi-model ensemble, suggesting that one may not be used as an analog for the other in estimates of uncertainty in future response. Whereas in the multi-model ensemble, the first order variation in cloud feedbacks shows compensation between longwave and shortwave components, in the perturbed physics ensemble the shortwave feedbacks are uncompensated, possibly explaining the larger range of climate sensitivities observed in the perturbed simulations. Regression analysis suggests that the parameters governing cloud formation, convection strength and ice fall speed are the most significant in altering climate feedbacks. Perturbations of oceanic and sulfur cycle parameters have relatively little effect on the atmospheric feedbacks diagnosed by the kernel technique.     

Using climate impacts indicators to evaluate climate model ensembles: temperature suitability of premium winegrape cultivation in the United States

We explore the potential to improve understanding of the climate system by directly targeting climate model analyses at specific indicators of climate change impact. Using the temperature suitability of premium winegrape cultivation as a climate impacts indicator, we quantify the inter- and intra-ensemble spread in three climate model ensembles: a physically uniform multi-member ensemble consisting of the RegCM3 high-resolution climate model nested within the NCAR CCSM3 global climate model; the multi-model NARCCAP ensemble consisting of single realizations of multiple high-resolution climate models nested within multiple global climate models; and the multi-model CMIP3 ensemble consisting of realizations of multiple global climate models. We find that the temperature suitability for premium winegrape cultivation is substantially reduced throughout the high-value growing areas of California and the Columbia Valley region (eastern Oregon and Washington) in all three ensembles in response to changes in temperature projected for the mid-twenty first century period. The reductions in temperature suitability are driven primarily by projected increases in mean growing season temperature and occurrence of growing season severe hot days. The intra-ensemble spread in the simulated climate change impact is smaller in the single-model ensemble than in the multi-model ensembles, suggesting that the uncertainty arising from internal climate system variability is smaller than the uncertainty arising from climate model formulation. In addition, the intra-ensemble spread is similar in the NARCCAP nested climate model ensemble and the CMIP3 global climate model ensemble, suggesting that the uncertainty arising from the model formulation of fine-scale climate processes is not smaller than the uncertainty arising from the formulation of large-scale climate processes. Correction of climate model biases substantially reduces both the inter- and intra-ensemble spread in projected climate change impact, particularly for the multi-model ensembles, suggesting that—at least for some systems—the projected impacts of climate change could be more robust than the projected climate change. Extension of this impacts-based analysis to a larger suite of impacts indicators will deepen our understanding of future climate change uncertainty by focusing on the climate phenomena that most directly influence natural and human systems.     

Towards quantifying uncertainty in transient climate change

Ensembles of coupled atmosphere–ocean global circulation model simulations are required to make probabilistic predictions of future climate change. “Perturbed physics” ensembles provide a new approach in which modelling uncertainties are sampled systematically by perturbing uncertain parameters. The aim is to provide a basis for probabilistic predictions in which the impact of prior assumptions and observational constraints can be clearly distinguished. Here we report on the first perturbed physics coupled atmosphere–ocean model ensemble in which poorly constrained atmosphere, land and sea-ice component parameters are varied in the third version of the Hadley Centre model (the variation of ocean parameters will be the subject of future study). Flux adjustments are employed, both to reduce regional sea surface temperature (SST) and salinity biases and also to admit the use of combinations of model parameter values which give non-zero values for the global radiation balance. This improves the extent to which the ensemble provides a credible basis for the quantification of uncertainties in climate change, especially at a regional level. However, this particular implementation of flux-adjustments leads to a weakening of the Atlantic overturning circulation, resulting in the development of biases in SST and sea ice in the North Atlantic and Arctic Oceans. Nevertheless, model versions are produced which are of similar quality to the unperturbed and un-flux-adjusted version. The ensemble is used to simulate pre-industrial conditions and a simple scenario of a 1% per year compounded increase in CO₂. The range of transient climate response (the 20 year averaged global warming at the time of CO₂ doubling) is 1.5–2.6°C, similar to that found in multi-model studies. Measures of global and large scale climate change from the coupled models show simple relationships with associated measures computed from atmosphere-mixed-layer-ocean climate change experiments, suggesting that recent advances in computing the probability density function of climate change under equilibrium conditions using the perturbed physics approach may be extended to the transient case.     

The impact of perturbations to ocean-model parameters on climate and climate change in a coupled model

Assessments of the impacts of uncertainties in parameters on mean climate and climate change in complex climate models have, to date, largely focussed on perturbations to parameters in the atmosphere component of the model. Here we expand on a previously published study which found the global impacts of perturbed ocean parameters on the rate of transient climate change to be small compared to perturbed atmosphere parameters. By separating the climate-change-induced ocean vertical heat transport in each perturbed member into components associated with the resolved flow and each parameterisation scheme, we show that variations in global mean heat uptake in different perturbed versions are an order of magnitude smaller than the average heat uptake. The lack of impact of the perturbations is attributed to (1) the relatively small impact of the perturbation on the direct vertical heat transport associated with the perturbed process and (2) a compensation between those direct changes and indirect changes in heat transport from other processes. Interactions between processes and changes appear to combine in complex ways to limit ensemble spread and uncertainty in the rate of warming. We also investigate regional impacts of the perturbations that may be important for climate change predictions. We find variations across the ensemble that are significant when measured against natural variability. In terms of the experimental set-up used here (models without flux adjustments) we conclude that perturbed physics ensembles with ocean parameter perturbations are an important component of any probabilistic estimate of future climate change, despite the low spread in global mean quantities. Hence, careful consideration should be given to assessing uncertainty in ocean processes in future probabilistic assessments of regional climate change.     

Reliability of multi-model and structurally different single-model ensembles

The performance of several state-of-the-art climate model ensembles, including two multi-model ensembles (MMEs) and four structurally different (perturbed parameter) single model ensembles (SMEs), are investigated for the first time using the rank histogram approach. In this method, the reliability of a model ensemble is evaluated from the point of view of whether the observations can be regarded as being sampled from the ensemble. Our analysis reveals that, in the MMEs, the climate variables we investigated are broadly reliable on the global scale, with a tendency towards overdispersion. On the other hand, in the SMEs, the reliability differs depending on the ensemble and variable field considered. In general, the mean state and historical trend of surface air temperature, and mean state of precipitation are reliable in the SMEs. However, variables such as sea level pressure or top-of-atmosphere clear-sky shortwave radiation do not cover a sufficiently wide range in some. It is not possible to assess whether this is a fundamental feature of SMEs generated with particular model, or a consequence of the algorithm used to select and perturb the values of the parameters. As under-dispersion is a potentially more serious issue when using ensembles to make projections, we recommend the application of rank histograms to assess reliability when designing and running perturbed physics SMEs.     

A multi-physics ensemble of regional climate change projections over the Iberian Peninsula

This study illustrates the sensitivity of regional climate change projections to the model physics. A single-model (MM5) multi-physics ensemble of regional climate simulations over the Iberian Peninsula for present (1970–1999) and future (2070–2099 under the A2 scenario) periods is assessed. The ensemble comprises eight members resulting from the combination of two options of parameterization schemes for the planetary boundary layer, cumulus and microphysics. All the considered combinations were previously evaluated by comparing hindcasted simulations to observations, none of them providing clearly outlying climates. Thus, the differences among the various ensemble members (spread) in the future projections could be considered as a matter of uncertainty in the change signals (as similarly assumed in multi-model studies). The results highlight the great dependence of the spread on the synoptic conditions driving the regional model. In particular, the spread generally amplifies under the future scenario leading to a large spread accompanying the mean change signals, as large as the magnitude of the mean projected changes and analogous to the spread obtained in multi-model ensembles. Moreover, the sign of the projected change varies depending on the choice of the model physics in many cases. This, together with the fact that the key mechanisms identified for the simulation of the climatology of a given period (either present or future) and those introducing the largest spread in the projected changes differ significantly, make further claims for efforts to better understand and model the parameterized subgrid processes.     

Quantifying uncertainties in projections of extremes��a perturbed land surface parameter experiment

Uncertainties in the climate response to a doubling of atmospheric CO₂ concentrations are quantified in a perturbed land surface parameter experiment. The ensemble of 108 members is constructed by systematically perturbing five poorly constrained land surface parameters of global climate model individually and in all possible combinations. The land surface parameters induce small uncertainties at global scale, substantial uncertainties at regional and seasonal scale and very large uncertainties in the tails of the distribution, the climate extremes. Climate sensitivity varies across the ensemble mainly due to the perturbation of the snow albedo parameterization, which controls the snow albedo feedback strength. The uncertainty range in the global response is small relative to perturbed physics experiments focusing on atmospheric parameters. However, land surface parameters are revealed to control the response not only of the mean but also of the variability of temperature. Major uncertainties are identified in the response of climate extremes to a doubling of CO₂. During winter the response both of temperature mean and daily variability relates to fractional snow cover. Cold extremes over high latitudes warm disproportionately in ensemble members with strong snow albedo feedback and large snow cover reduction. Reduced snow cover leads to more winter warming and stronger variability decrease. As a result uncertainties in mean and variability response line up, with some members showing weak and others very strong warming of the cold tail of the distribution, depending on the snow albedo parametrization. The uncertainty across the ensemble regionally exceeds the CMIP3 multi-model range. Regarding summer hot extremes, the uncertainties are larger than for mean summer warming but smaller than in multi-model experiments. The summer precipitation response to a doubling of CO₂ is not robust over many regions. Land surface parameter perturbations and natural variability alter the sign of the response even over subtropical regions.     

Reliability and importance of structural diversity of climate model ensembles

We investigate the performance of the newest generation multi-model ensemble (MME) from the Coupled Model Intercomparison Project (CMIP5). We compare the ensemble to the previous generation models (CMIP3) as well as several single model ensembles (SMEs), which are constructed by varying components of single models. These SMEs range from ensembles where parameter uncertainties are sampled (perturbed physics ensembles) through to an ensemble where a number of the physical schemes are switched (multi-physics ensemble). We focus on assessing reliability against present-day climatology with rank histograms, but also investigate the effective degrees of freedom (EDoF) of the fields of variables which makes the statistical test of reliability more rigorous, and consider the distances between the observation and ensemble members. We find that the features of the CMIP5 rank histograms, of general reliability on broad scales, are consistent with those of CMIP3, suggesting a similar level of performance for present-day climatology. The spread of MMEs tends towards being “over-dispersed” rather than “under-dispersed”. In general, the SMEs examined tend towards insufficient dispersion and the rank histogram analysis identifies them as being statistically distinguishable from many of the observations. The EDoFs of the MMEs are generally greater than those of SMEs, suggesting that structural changes lead to a characteristically richer range of model behaviours than is obtained with parametric/physical-scheme-switching ensembles. For distance measures, the observations and models ensemble members are similarly spaced from each other for MMEs, whereas for the SMEs, the observations are generally well outside the ensemble. We suggest that multi-model ensembles should represent an important component of uncertainty analysis.     

On the generation of climate model ensembles

Climate model ensembles are used to estimate uncertainty in future projections, typically by interpreting the ensemble distribution for a particular variable probabilistically. There are, however, different ways to produce climate model ensembles that yield different results, and therefore different probabilities for a future change in a variable. Perhaps equally importantly, there are different approaches to interpreting the ensemble distribution that lead to different conclusions. Here we use a reduced-resolution climate system model to compare three common ways to generate ensembles: initial conditions perturbation, physical parameter perturbation, and structural changes. Despite these three approaches conceptually representing very different categories of uncertainty within a modelling system, when comparing simulations to observations of surface air temperature they can be very difficult to separate. Using the twentieth century CMIP5 ensemble for comparison, we show that initial conditions ensembles, in theory representing internal variability, significantly underestimate observed variance. Structural ensembles, perhaps less surprisingly, exhibit over-dispersion in simulated variance. We argue that future climate model ensembles may need to include parameter or structural perturbation members in addition to perturbed initial conditions members to ensure that they sample uncertainty due to internal variability more completely. We note that where ensembles are over- or under-dispersive, such as for the CMIP5 ensemble, estimates of uncertainty need to be treated with care.     

The sensitivity of the meridional overturning circulation to modelling uncertainty in a perturbed physics ensemble without flux adjustment

Ocean dynamics play a key role in the climate system, by redistributing heat and freshwater. The uncertainty of how these processes are represented in climate models, and how this uncertainty affects future climate projections can be investigated using perturbed physics ensembles of global circulation models (GCMs). Techniques such as flux adjustments should be avoided since they can impact the sensitivity of the ensemble to the imposed forcing. In this study a method for developing an coupled ensemble with a GCM that does not use flux adjustment is presented. The ensemble is constrained by using information from a prior ensemble with a mixed layer ocean coupled to an atmosphere GCM, to reduce drifts in the coupled ensemble. Constraints on parameter perturbations are derived by using observational constraints on surface temperature, and top of the atmosphere radiative fluxes. As an example of such an ensemble developed with this methodology, uncertainty in response of the meridional overturning circulation (MOC) to increased CO₂ concentrations is investigated. The ensemble mean MOC strength is 17.1?Sv and decreases by 2.1?Sv when greenhouse gas concentrations are doubled. No rapid changes or shutdown of the MOC are seen in any of the ensemble members. There is a strong negative relationship between global mean temperature and MOC strength across the ensemble which is not seen in a multimodel ensemble. A positive relationship between climate sensitivity and the decrease of MOC strength is also seen.     

Quantification of uncertainty sources in a probabilistic climate change assessment of future water shortages

As the incorporation of probabilistic climate change information into UK water resource management gathers apace, understanding the relative scales of the uncertainty sources in projections of future water shortage metrics is necessary for the resultant information to be understood and used effectively. Utilising modified UKCP09 weather generator data and a multi-model approach, this paper represents a first attempt at extending an uncertainty assessment of future stream flows under forced climates to consider metrics of water shortage based on the triggering of reservoir control curves. It is found that the perturbed physics ensemble uncertainty, which describes climate model parameter error uncertainty, is the cause of a far greater proportion of both the overall flow and water shortage per year probability uncertainty than that caused by SRES emissions scenario choice in the 2080s. The methodology for producing metrics of future water shortage risk from UKCP09 weather generator information described here acts as the basis of a robustness analysis of the North Staffordshire WRZ to climate change, which provides an alternative approach for making decisions despite large uncertainties, which will follow.     

Climate model dependence and the replicate Earth paradigm

Multi-model ensembles are commonly used in climate prediction to create a set of independent estimates, and so better gauge the likelihood of particular outcomes and better quantify prediction uncertainty. Yet researchers share literature, datasets and model code—to what extent do different simulations constitute independent estimates? What is the relationship between model performance and independence? We show that error correlation provides a natural empirical basis for defining model dependence and derive a weighting strategy that accounts for dependence in experiments where the multi-model mean would otherwise be used. We introduce the “replicate Earth” ensemble interpretation framework, based on theoretically derived statistical relationships between ensembles of perfect models (replicate Earths) and observations. We transform an ensemble of (imperfect) climate projections into an ensemble whose mean and variance have the same statistical relationship to observations as an ensemble of replicate Earths. The approach can be used with multi-model ensembles that have varying numbers of simulations from different models, accounting for model dependence. We use HadCRUT3 data and the CMIP3 models to show that in out of sample tests, the transformed ensemble has an ensemble mean with significantly lower error and much flatter rank frequency histograms than the original ensemble.     

CMIP5多模式对东亚地区地面气温年际变率的回报研究

对CMIP5全球气候模式中年代际回报试验的气温资料及其简单集合平均(Multi-model ensemble mean,EMN)和贝叶斯模式平均的结果(Bayesian Model Averaging,BMA)进行经验正交函数(Empirical Orthogonal Function,EOF)分解和Morlet小波分析,检验评估各个模式及其EMN和BMA对东亚地面气温的方差、气温时空分布特征及周期变化的回报能力。结果表明,10个模式、EMN、BMA都能很好地回报出1981—2010年东亚地面气温的方差分布,其中BMA回报效果最好。EOF分析表明,BMA能较好地回报出东亚地面气温第一模态的时空分布。MIROC5能较好地回报出第二模态的趋势变化,但却不能回报出气温的年际变率。绝大多数模式和EMN、BMA虽然能回报出东亚地面气温的变化趋势,但是对气温年际变率的回报仍然是比较困难的。CMCC-CM对气温变化主模态的3~5 a的周期变化特征回报效果最好,和NCEP资料的结果最为接近。     

Projected increased risk of water deficit over major West African river basins under future climates

Despite decades of research, large multi-model uncertainty remains about the Earth’s equilibrium climate sensitivity to carbon dioxide forcing as inferred from state-of-the-art Earth system models (ESMs). Statistical treatments of multi-model uncertainties are often limited to simple ESM averaging approaches. Sometimes models are weighted by how well they reproduce historical climate observations. Here, we propose a novel approach to multi-model combination and uncertainty quantification. Rather than averaging a discrete set of models, our approach samples from a continuous distribution over a reduced space of simple model parameters. We fit the free parameters of a reduced-order climate model to the output of each member of the multi-model ensemble. The reduced-order parameter estimates are then combined using a hierarchical Bayesian statistical model. The result is a multi-model distribution of reduced-model parameters, including climate sensitivity. In effect, the multi-model uncertainty problem within an ensemble of ESMs is converted to a parametric uncertainty problem within a reduced model. The multi-model distribution can then be updated with observational data, combining two independent lines of evidence. We apply this approach to 24 model simulations of global surface temperature and net top-of-atmosphere radiation response to abrupt quadrupling of carbon dioxide, and four historical temperature data sets. Our reduced order model is a 2-layer energy balance model. We present probability distributions of climate sensitivity based on (1) the multi-model ensemble alone and (2) the multi-model ensemble and observations.     

The benefits of quantifying climate model uncertainty in climate change impacts assessment: an example with heat-related mortality change estimates

The majority of climate change impacts assessments account for climate change uncertainty by adopting the scenario-based approach. This typically involves assessing the impacts for a small number of emissions scenarios but neglecting the role of climate model physics uncertainty. Perturbed physics ensemble (PPE) climate simulations offer a unique opportunity to explore this uncertainty. Furthermore, PPEs mean it is now possible to make risk-based impacts estimates because they allow for a range of estimates to be presented to decision-makers, which spans the range of climate model physics uncertainty inherent from a given climate model and emissions scenario, due to uncertainty associated with the understanding of physical processes in the climate model. This is generally not possible with the scenario-based approach. Here, we present the first application of a PPE to estimate the impact of climate change on heat-related mortality. By using the estimated impacts of climate change on heat-related mortality in six cities, we demonstrate the benefits of quantifying climate model physics uncertainty in climate change impacts assessment over the more common scenario-based approach. We also show that the impacts are more sensitive to climate model physics uncertainty than they are to emissions scenario uncertainty, and least sensitive to whether the climate change projections are from a global climate model or a regional climate model. The results demonstrate the importance of presenting model uncertainties in climate change impacts assessments if the impacts are to be placed within a climate risk management framework.     

On the estimation of systematic error in regression-based predictions of climate sensitivity

An extension of a regression-based methodology for constraining climate forecasts using a multi-thousand member ensemble of perturbed climate models is presented, using the multi-model CMIP-3 ensemble to estimate the systematic model uncertainty in the prediction, with the caveat that systematic biases common to all models are not accounted for. It is shown that previous methodologies for estimating the systematic uncertainty in predictions of climate sensitivity are dependent on arbitrary choices relating to ensemble sampling strategy. Using a constrained regression approach, a multivariate predictor may be derived based upon the mean climatic state of each ensemble member, but components of this predictor are excluded if they cannot be validated within the CMIP-3 ensemble. It is found that the application of the CMIP-3 constraint serves to decrease the upper bound of likelihood for climate sensitivity when compared with previous studies, with 10th and 90th percentiles of probability at 1.5 K and 4.3 K respectively.     

Regional and seasonal intercomparison of CMIP3 and CMIP5 climate model ensembles for temperature and precipitation

Regional and seasonal temperature and precipitation over land are compared across two generations of global climate model ensembles, specifically, CMIP5 and CMIP3, through historical twentieth century skills and multi-model agreement, and twenty first century projections. A suite of diagnostic and performance metrics, ranging from spatial bias or model-consensus maps and aggregate time series plots, to measures of equivalence between probability density functions and Taylor diagrams, are used for the intercomparisons. Pairwise and multi-model ensemble comparisons were performed for 11 models, which were selected based on data availability and resolutions. Results suggest little change in the central tendency or variability or uncertainty of historical skills or consensus across the two generations of models. However, there are regions and seasons, at different levels of aggregation, where significant changes, performance improvements, and even degradation in skills, are suggested. The insights may provide directions for further improvements in next generations of climate models, and in the meantime, help inform adaptation and policy.     

Identifying uncertainties in Arctic climate change projections

Wide ranging climate changes are expected in the Arctic by the end of the 21st century, but projections of the size of these changes vary widely across current global climate models. This variation represents a large source of uncertainty in our understanding of the evolution of Arctic climate. Here we systematically quantify and assess the model uncertainty in Arctic climate changes in two CO₂ doubling experiments: a multimodel ensemble (CMIP3) and an ensemble constructed using a single model (HadCM3) with multiple parameter perturbations (THC-QUMP). These two ensembles allow us to assess the contribution that both structural and parameter variations across models make to the total uncertainty and to begin to attribute sources of uncertainty in projected changes. We find that parameter uncertainty is an major source of uncertainty in certain aspects of Arctic climate. But also that uncertainties in the mean climate state in the 20th century, most notably in the northward Atlantic ocean heat transport and Arctic sea ice volume, are a significant source of uncertainty for projections of future Arctic change. We suggest that better observational constraints on these quantities will lead to significant improvements in the precision of projections of future Arctic climate change.     

Future changes in tropical cyclone genesis in fully dynamic ocean- and mixed layer ocean-coupled climate models: a low-resolution model study

The global and regional projected changes in tropical cyclone (TC) genesis due to increased CO₂ concentrations has been investigated through a large-scale TC genesis parameter (convective seasonal genesis parameter, ConvGP) in two perturbed physics ensembles. The ensembles are based on the third generation Hadley Centre atmosphere?Cocean general circulation model with the first ensemble using a coupled fully dynamic ocean (HadCM3) and the second coupled to a simplified mixed layer thermodynamic ocean (HadSM3) both consisting of 17 members. In each ensemble, parameters are identically perturbed to provide a wide range of climate sensitivity whilst retaining a credible present-day climate simulation. It is found, by comparing the ConvGP climatology from reanalysis data with the best track genesis, that it is possible to reproduce the observed genesis distribution. Future changes in the spatial ConvGP distribution are explored with respect to each tropical ocean basin. Whilst there is a similarity in the gross pattern of the ensemble-mean projected ConvGP change between HadCM3 and HadSM3, there is a non-trivial difference in the tropical Pacific Ocean, arising from different patterns of tropical Pacific sea surface temperature change. This indicates that ocean representation can be important for regional scale projections. The quantitative contribution of individual constituent parameters (i.e. vorticity parameter, shear parameter and convective potential) to the projected ConvGP change is estimated. It is found that all three large-scale parameters generally contribute constructively, but with different magnitude, in the regions where a large doubled CO₂ response is found.     

Comparison of four ensemble methods combining regional climate simulations over Asia

A number of uncertainties exist in climate simulation because the results of climate models are influenced by factors such as their dynamic framework, physical processes, initial and driving fields, and horizontal and vertical resolution. The uncertainties of the model results may be reduced, and the credibility can be improved by employing multi-model ensembles. In this paper, multi-model ensemble results using 10-year simulations of five regional climate models (RCMs) from December 1988 to November 1998 over Asia are presented and compared. The simulation results are derived from phase II of the Regional Climate Model Inter-comparison Project (RMIP) for Asia. Using the methods of the arithmetic mean, the weighted mean, multivariate linear regression, and singular value decomposition, the ensembles for temperature, precipitation, and sea level pressure are carried out. The results show that the multi-RCM ensembles outperform the single RCMs in many aspects. Among the four ensemble methods used, the multivariate linear regression, based on the minimization of the root mean square errors, significantly improved the ensemble results. With regard to the spatial distribution of the mean climate, the ensemble result for temperature was better than that for precipitation. With an increasing number of models used in the ensembles, the ensemble results were more accurate. Therefore, a multi-model ensemble is an efficient approach to improve the results of regional climate simulations.