An efficient climate forecasting method using an intermediate complexity Earth System Model and the ensemble Kalman filter

We present the implementation and results of a model tuning and ensemble forecasting experiment using an ensemble Kalman filter for the simultaneous estimation of 12 parameters in a low resolution coupled atmosphere-ocean Earth System Model by tuning it to realistic data sets consisting of Levitus ocean temperature/salinity climatology, and NCEP/NCAR atmospheric temperature/humidity reanalysis data. The resulting ensemble of tuned model states is validated by comparing various diagnostics, such as mass and heat transports, to observational estimates and other model results. We show that this ensemble has a very reasonable climatology, with the 3-D ocean in particular having comparable realism to much more expensive coupled numerical models, at least in respect of these averaged indicators. A simple global warming experiment is performed to investigate the response and predictability of the climate to a change in radiative forcing, due to 100 years of 1% per annum atmospheric CO₂ increase. The equilibrium surface air temperature rise for this CO₂ increase is 4.2±0.1°C, which is approached on a time scale of 1,000 years. The simple atmosphere in this version of the model is missing several factors which, if included, would substantially increase the uncertainty of this estimate. However, even within this ensemble, there is substantial regional variability due to the possibility of collapse of the North Atlantic thermohaline circulation (THC), which switches off in more than one third of the ensemble members. For these cases, the regional temperature is not only 3–5°C colder than in the warmed worlds where the THC remains switched on, but is also 1–2°C colder than the current climate. Our results, which illustrate how objective probabilistic projections of future climate change can be efficiently generated, indicate a substantial uncertainty in the long-term future of the THC, and therefore the regional climate of western Europe. However, this uncertainty is only apparent in long-term integrations, with the initial transient response being similar across the entire ensemble. Application of this ensemble Kalman filtering technique to more complete climate models would improve the objectivity of probabilistic forecasts and hence should lead to significantly increased understanding of the uncertainty of our future climate.     

Climate model benchmarking with glacial and mid-Holocene climates

Past climates provide a test of models’ ability to predict climate change. We present a comprehensive evaluation of state-of-the-art models against Last Glacial Maximum and mid-Holocene climates, using reconstructions of land and ocean climates and simulations from the Palaeoclimate Modelling and Coupled Modelling Intercomparison Projects. Newer models do not perform better than earlier versions despite higher resolution and complexity. Differences in climate sensitivity only weakly account for differences in model performance. In the glacial, models consistently underestimate land cooling (especially in winter) and overestimate ocean surface cooling (especially in the tropics). In the mid-Holocene, models generally underestimate the precipitation increase in the northern monsoon regions, and overestimate summer warming in central Eurasia. Models generally capture large-scale gradients of climate change but have more limited ability to reproduce spatial patterns. Despite these common biases, some models perform better than others.     

200?years of European temperature variability: insights from and tests of the proxy surrogate reconstruction analog method

Spatially resolved climate reconstructions are commonly derived from long instrumental series and proxy data via linear regression based approaches that use the main modes of the climate system. Such reconstructions have been shown to underestimate climate variability and are based upon the assumption that the main modes of climate variability are stationary back in time. Climate models simulate physically consistent climate fields but cannot be taken to represent the real past climate trajectory because of their necessarily simplified scope and chaotic internal variability. Here, we present sensitivity tests of, and a 200-year temperature reconstruction from, the PSR (Proxy Surrogate Reconstruction) method. This method simultaneously capitalizes on the individual strengths of instrumental/proxy data based reconstructions and model simulations by selecting the model states (analogs) that are most similar with proxy/instrumental data available at specific places and specific moments of time. Sensitivity experiments reveal an optimal PSR configuration and indicate that 6,500 simulation years of existing climate models provide a sufficient pool of possible analogs to skillfully reconstruct monthly European temperature fields during the past 200?years. Reconstruction verification based upon only seven instrumental stations indicates potential for extensions back in time using sparse proxy data. Additionally the PSR method allows evaluation of single time series, in this case the homogeneity of instrumental series, by identifying inconsistencies with the reconstructed climate field. We present an updated European temperature reconstruction including newly homogenized instrumental records performed with the computationally efficient PSR method that proves to capture the total variance of the target.     

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资料的结果最为接近。     

Future projections of the surface heat and water budgets of the Mediterranean Sea in an ensemble of coupled atmosphere–ocean regional climate models

Within the CIRCE project “Climate change and Impact Research: the Mediterranean Environment”, an ensemble of high resolution coupled atmosphere–ocean regional climate models (AORCMs) are used to simulate the Mediterranean climate for the period 1950–2050. For the first time, realistic net surface air-sea fluxes are obtained. The sea surface temperature (SST) variability is consistent with the atmospheric forcing above it and oceanic constraints. The surface fluxes respond to external forcing under a warming climate and show an equivalent trend in all models. This study focuses on the present day and on the evolution of the heat and water budget over the Mediterranean Sea under the SRES-A1B scenario. On the contrary to previous studies, the net total heat budget is negative over the present period in all AORCMs and satisfies the heat closure budget controlled by a net positive heat gain at the strait of Gibraltar in the present climate. Under climate change scenario, some models predict a warming of the Mediterranean Sea from the ocean surface (positive net heat flux) in addition to the positive flux at the strait of Gibraltar for the 2021–2050 period. The shortwave and latent flux are increasing and the longwave and sensible fluxes are decreasing compared to the 1961–1990 period due to a reduction of the cloud cover and an increase in greenhouse gases (GHGs) and SSTs over the 2021–2050 period. The AORCMs provide a good estimates of the water budget with a drying of the region during the twenty-first century. For the ensemble mean, he decrease in precipitation and runoff is about 10 and 15% respectively and the increase in evaporation is much weaker, about 2% compared to the 1961–1990 period which confirm results obtained in recent studies. Despite a clear consistency in the trends and results between the models, this study also underlines important differences in the model set-ups, methodology and choices of some physical parameters inducing some difference in the various air-sea fluxes. An evaluation of the uncertainty sources and possible improvement for future generation of AORCMs highlights the importance of the parameterisation of the ocean albedo, rivers and cloud cover.     

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.     

The role of boundary organizations in climate change adaptation from the perspective of municipal practitioners

In climate change impact research it is crucial to carefully select the meteorological input for impact models. We present a method for model selection that enables the user to shrink the ensemble to a few representative members, conserving the model spread and accounting for model similarity. This is done in three steps: First, using principal component analysis for a multitude of meteorological parameters, to find common patterns of climate change within the multi-model ensemble. Second, detecting model similarities with regard to these multivariate patterns using cluster analysis. And third, sampling models from each cluster, to generate a subset of representative simulations. We present an application based on the ENSEMBLES regional multi-model ensemble with the aim to provide input for a variety of climate impact studies. We find that the two most dominant patterns of climate change relate to temperature and humidity patterns. The ensemble can be reduced from 25 to 5 simulations while still maintaining its essential characteristics. Having such a representative subset of simulations reduces computational costs for climate impact modeling and enhances the quality of the ensemble at the same time, as it prevents double-counting of dependent simulations that would lead to biased statistics.     

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.     

Climate sensitivity estimated from ensemble simulations of glacial climate

The concentration of greenhouse gases (GHGs) in the atmosphere continues to rise, hence estimating the climate system’s sensitivity to changes in GHG concentration is of vital importance. Uncertainty in climate sensitivity is a main source of uncertainty in projections of future climate change. Here we present a new approach for constraining this key uncertainty by combining ensemble simulations of the last glacial maximum (LGM) with paleo-data. For this purpose we used a climate model of intermediate complexity to perform a large set of equilibrium runs for (1) pre-industrial boundary conditions, (2) doubled CO₂ concentrations, and (3) a complete set of glacial forcings (including dust and vegetation changes). Using proxy-data from the LGM at low and high latitudes we constrain the set of realistic model versions and thus climate sensitivity. We show that irrespective of uncertainties in model parameters and feedback strengths, in our model a close link exists between the simulated warming due to a doubling of CO₂, and the cooling obtained for the LGM. Our results agree with recent studies that annual mean data-constraints from present day climate prove to not rule out climate sensitivities above the widely assumed sensitivity range of 1.5–4.5°C (Houghton et al. 2001). Based on our inferred close relationship between past and future temperature evolution, our study suggests that paleo-climatic data can help to reduce uncertainty in future climate projections. Our inferred uncertainty range for climate sensitivity, constrained by paleo-data, is 1.2–4.3°C and thus almost identical to the IPCC estimate. When additionally accounting for potential structural uncertainties inferred from other models the upper limit increases by about 1°C.     

Weak oceanic heat transport as a cause of the instability of glacial climates

The stability of the thermohaline circulation of modern and glacial climates is compared with the help of a two dimensional ocean—atmosphere—sea ice coupled model. It turns out to be more unstable as less freshwater forcing is required to induce a polar halocline catastrophy in glacial climates. The large insulation of the ocean by the extensive sea ice cover changes the temperature boundary condition and the deepwater formation regions moves much further South. The nature of the instability is of oceanic origin, identical to that found in ocean models under mixed boundary conditions. With similar strengths of the oceanic circulation and rates of deep water formation for warm and cold climates, the loss of stability of the cold climate is due to the weak thermal stratification caused by the cooling of surface waters, the deep water temperatures being regulated by the temperature of freezing. Weaker stratification with similar overturning leads to a weakening of the meridional oceanic heat transport which is the major negative feedback stabilizing the oceanic circulation. Within the unstable regime periodic millennial oscillations occur spontaneously. The climate oscillates between a strong convective thermally driven oceanic state and a weak one driven by large salinity gradients. Both states are unstable. The atmosphere of low thermal inertia is carried along by the oceanic overturning while the variation of sea ice is out of phase with the oceanic heat content. During the abrupt warming events that punctuate the course of a millennial oscillation, sea ice variations are shown respectively to damp (amplify) the amplitude of the oceanic (atmospheric) response. This sensitivity of the oceanic circulation to a reduced concentration of greenhouse gases and to freshwater forcing adds support to the hypothesis that the millennial oscillations of the last glacial period, the so called Dansgaard—Oeschger events, may be internal instabilities of the climate system.     

Past and future polar amplification of climate change: climate model intercomparisons and ice-core constraints

Climate model simulations available from the PMIP1, PMIP2 and CMIP (IPCC-AR4) intercomparison projects for past and future climate change simulations are examined in terms of polar temperature changes in comparison to global temperature changes and with respect to pre-industrial reference simulations. For the mid-Holocene (MH, 6,000 years ago), the models are forced by changes in the Earth’s orbital parameters. The MH PMIP1 atmosphere-only simulations conducted with sea surface temperatures fixed to modern conditions show no MH consistent response for the poles, whereas the new PMIP2 coupled atmosphere–ocean climate models systematically simulate a significant MH warming both for Greenland (but smaller than ice-core based estimates) and Antarctica (consistent with the range of ice-core based range). In both PMIP1 and PMIP2, the MH annual mean changes in global temperature are negligible, consistent with the MH orbital forcing. The simulated last glacial maximum (LGM, 21,000 years ago) to pre-industrial change in global mean temperature ranges between 3 and 7°C in PMIP1 and PMIP2 model runs, similar to the range of temperature change expected from a quadrupling of atmospheric CO₂ concentrations in the CMIP simulations. Both LGM and future climate simulations are associated with a polar amplification of climate change. The range of glacial polar amplification in Greenland is strongly dependent on the ice sheet elevation changes prescribed to the climate models. All PMIP2 simulations systematically underestimate the reconstructed glacial–interglacial Greenland temperature change, while some of the simulations do capture the reconstructed glacial–interglacial Antarctic temperature change. Uncertainties in the prescribed central ice cap elevation cannot account for the temperature change underestimation by climate models. The variety of climate model sensitivities enables the exploration of the relative changes in polar temperature with respect to changes in global temperatures. Simulated changes of polar temperatures are strongly related to changes in simulated global temperatures for both future and LGM climates, confirming that ice-core-based reconstructions provide quantitative insights on global climate changes. An erratum to this article can be found at     

Climate projections of future extreme events accounting for modelling uncertainties and historical simulation biases

A methodology is presented for providing projections of absolute future values of extreme weather events that takes into account key uncertainties in predicting future climate. This is achieved by characterising both observed and modelled extremes with a single form of non-stationary extreme value (EV) distribution that depends on global mean temperature and which includes terms that account for model bias. Such a distribution allows the prediction of future “observed” extremes for any period in the twenty-first century. Uncertainty in modelling future climate, arising from a wide range of atmospheric, oceanic, sulphur cycle and carbon cycle processes, is accounted for by using probabilistic distributions of future global temperature and EV parameters. These distributions are generated by Bayesian sampling of emulators with samples weighted by their likelihood with respect to a set of observational constraints. The emulators are trained on a large perturbed parameter ensemble of global simulations of the recent past, and the equilibrium response to doubled CO₂. Emulated global EV parameters are converted to the relevant regional scale through downscaling relationships derived from a smaller perturbed parameter regional climate model ensemble. The simultaneous fitting of the EV model to regional model data and observations allows the characterisation of how observed extremes may change in the future irrespective of biases that may be present in the regional models simulation of the recent past climate. The clearest impact of a parameter perturbation in this ensemble was found to be the depth to which plants can access water. Members with shallow soils tend to be biased hot and dry in summer for the observational period. These biases also appear to have an impact on the potential future response for summer temperatures with some members with shallow soils having increases for extremes that reduce with extreme severity. We apply this methodology for London, using the A1B future emissions scenario to obtain projections of the 50 year return values for the 20 year period centred on 2050. We obtain 10th to 90th percentile ranges of 35.9–42.1 °C for summer daily maximum temperature, 35.5–52.4 mm for summer daily rainfall and 79.2, 97.0 mm for autumn 5 day total rainfall, compared to observed estimates for 1961–1990 of 35.7 °C, 42.1 and 78.4 mm respectively.     

Impact of Greenland and Antarctic ice sheet interactions on climate sensitivity

We use the Earth system model of intermediate complexity LOVECLIM to show the effect of coupling interactive ice sheets on the climate sensitivity of the model on a millennial time scale. We compare the response to a 2×CO₂ warming scenario between fully coupled model versions including interactive Greenland and Antarctic ice sheet models and model versions with fixed ice sheets. For this purpose an ensemble of different parameter sets have been defined for LOVECLIM, covering a wide range of the model??s sensitivity to greenhouse warming, while still simulating the present-day climate and the climate evolution over the last millennium within observational uncertainties. Additional freshwater fluxes from the melting ice sheets have a mitigating effect on the model??s temperature response, leading to generally lower climate sensitivities of the fully coupled model versions. The mitigation is effectuated by changes in heat exchange within the ocean and at the sea?Cair interface, driven by freshening of the surface ocean and amplified by sea?Cice-related feedbacks. The strength of the effect depends on the response of the ice sheets to the warming and on the model??s climate sensitivity itself. The effect is relatively strong in model versions with higher climate sensitivity due to the relatively large polar amplification of LOVECLIM. With the ensemble approach in this study we cover a wide range of possible model responses.     

How dynamical models can learn from the data��an example with a simplified ENSO model

Most dynamical models of the natural system contain a number of empirical parameters which reflect our limited understanding of the simulated system or describe unresolved subgrid-scale processes. While the parameterizations basically introduce some uncertainty to the model results, they also hold the prospect of tuning the model. In general, a deterministic tuning is related to an inversion of the model which is often impossible or requires considerable computing effort for most climate models. Another way to adjust the model parameters to a specific observed process is stochastic fitting where a set of parameters and model output are taken as random variables. Here, we present a dynamical?Cstatistical approach with a simplified model of the El Ni?o?CSouthern Oscillation (ENSO) cycle whose parameters are adjusted to simulated and observed data by means of Bayesian statistics. As ENSO model, we employ the Schop?CSuarez delay oscillator model. Monte Carlo experiments highlight the large sensitivity of the model results to varied model parameters and initial values. The statistical adjustment is done by Bayesian model averaging of the Monte Carlo experiments. Applying the method to simulated data, the posterior ensemble mean is much closer to the reference data than the prior ensemble mean. The learning effect of the model is evident in the leading empirical orthogonal functions and statistically significant in the mean state. When the method is applied to the observed ENSO time series, the ENSO model in its classical setup is not able to account for the temporally varying periodicity of the observed ENSO phenomenon. An improved setup with continuous adjustment periods and extended parameter range is developed in order to allow the model to learn from the data gradually. The improved setup leads to promising results during the twentieth century and even a weak forecast skill over 6?months. Thus, the described method offers a promising tool for data assimilation in dynamical weather and climate models. However, the simplified ENSO model is barely appropriate for operational ENSO forecasts owing to its limited physical complexity.     

A modified Köppen classification applied to model simulations of glacial and interglacial climates

A series of experiments was done using an atmospheric general circulation model to simulate climates from full glacial time at 18 ka (thousands of years before the present) to the present at 3000 year intervals, and at 126 ka, the previous interglacial period. A modified Köppen climate classification was developed to aid in the interpretation of the results of the circulation model experiments. The climate classification scheme permits the characterization of eleven distinct seasonal temperature and precipitation regimes. For the modern climate, the modified classification agrees well with a classification of natural vegetation zones, and provides an easily-assimilated depiction of climate changes resulting from the varying boundary conditions in the past. At 18 ka, the time of glacial maximum, 45% of the land surface had climate classifications different from the present. At 126 ka, a time when northern hemisphere summer radiation was much greater than at present owing to changes in the date of perihelion and tilt of the earth's axis, the corresponding difference was 32%. For all experiments -3 to 18 ka and 126 ka - only 30% of the land surface showed no change in climate classification from the present. Core areas showing no change included the Amazon basin, the northern Sahara and Australia.     

CMIP1 evaluation and intercomparison of coupled climate models

 The climates simulated by 15 coupled atmosphere/ocean climate models participating in the first phase of the Coupled Model Intercomparison Project (CMIP1) are intercompared and evaluated. Results for global means, zonal averages, and geographical distributions of basic climate variables are assembled and compared with observations. The current generation of climate models reproduce the major features of the observed distribution of the basic climate parameters, but there is, nevertheless, a considerable scatter among model results and between simulated and observed values. This is particularly true for oceanic variables. Flux adjusted models generally produce simulated climates which are in better accord with observations than do non-flux adjusted models; however, some non-flux adjusted model results are closer to observations than some flux adjusted model results. Other model differences, such as resolution, do not appear to provide a clear distinction among model results in this generation of models. Many of the systematic differences (those differences common to most models), evident in previous intercomparison studies are exhibited also by the CMIP1 group of models although often with reduced magnitudes. As is characteristic of intercomparison results, different climate variables are simulated with different levels of success by different models and no one model is “best” for all variables. There is some evidence that the “mean model” result, obtained by averaging over the ensemble of models, provides an overall best comparison to observations for climatological mean fields. The model deficiencies identified here do not suggest immediate remedies and the overall success of the models in simulating the behaviour of the complex non-linear climate system apparently depends on the slow improvement in the balance of approximations that characterize a coupled climate model. Of course, the results of this and similar studies provide only an indication, at a particular time, of the current state and the moderate but steady evolution and improvement of coupled climate models. Received: 26 January 2000 / Accepted: 9 June 2000     

Transient climate changes in a perturbed parameter ensemble of emissions-driven earth system model simulations

We describe results from a 57-member ensemble of transient climate change simulations, featuring simultaneous perturbations to 54 parameters in the atmosphere, ocean, sulphur cycle and terrestrial ecosystem components of an earth system model (ESM). These emissions-driven simulations are compared against the CMIP3 multi-model ensemble of physical climate system models, used extensively to inform previous assessments of regional climate change, and also against emissions-driven simulations from ESMs contributed to the CMIP5 archive. Members of our earth system perturbed parameter ensemble (ESPPE) are competitive with CMIP3 and CMIP5 models in their simulations of historical climate. In particular, they perform reasonably well in comparison with HadGEM2-ES, a more sophisticated and expensive earth system model contributed to CMIP5. The ESPPE therefore provides a computationally cost-effective tool to explore interactions between earth system processes. In response to a non-intervention emissions scenario, the ESPPE simulates distributions of future regional temperature change characterised by wide ranges, and warm shifts, compared to those of CMIP3 models. These differences partly reflect the uncertain influence of global carbon cycle feedbacks in the ESPPE. In addition, the regional effects of interactions between different earth system feedbacks, particularly involving physical and ecosystem processes, shift and widen the ESPPE spread in normalised patterns of surface temperature and precipitation change in many regions. Significant differences from CMIP3 also arise from the use of parametric perturbations (rather than a multimodel ensemble) to represent model uncertainties, and this is also the case when ESPPE results are compared against parallel emissions-driven simulations from CMIP5 ESMs. When driven by an aggressive mitigation scenario, the ESPPE and HadGEM2-ES reveal significant but uncertain impacts in limiting temperature increases during the second half of the twenty-first century. Emissions-driven simulations create scope for development of errors in properties that were previously prescribed in coupled ocean–atmosphere models, such as historical CO₂ concentrations and vegetation distributions. In this context, historical intra-ensemble variations in the airborne fraction of CO₂ emissions, and in summer soil moisture in northern hemisphere continental regions, are shown to be potentially useful constraints, subject to uncertainties in the relevant observations. Our results suggest that future climate-related risks can be assessed more comprehensively by updating projection methodologies to support formal combination of emissions-driven perturbed parameter and multi-model earth system model simulations with suitable observational constraints. This would provide scenarios underpinned by a more complete representation of the chain of uncertainties from anthropogenic emissions to future climate outcomes.     

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.     

Uncertainties due to transport-parameter sensitivity in an efficient 3-D ocean-climate model

A simplified climate model is presented which includes a fully 3-D, frictional geostrophic (FG) ocean component but retains an integration efficiency considerably greater than extant climate models with 3-D, primitive-equation ocean representations (20 kyears of integration can be completed in about a day on a PC). The model also includes an Energy and Moisture Balance atmosphere and a dynamic and thermodynamic sea-ice model. Using a semi-random ensemble of 1,000 simulations, we address both the inverse problem of parameter estimation, and the direct problem of quantifying the uncertainty due to mixing and transport parameters. Our results represent a first attempt at tuning a 3-D climate model by a strictly defined procedure, which nevertheless considers the whole of the appropriate parameter space. Model estimates of meridional overturning and Atlantic heat transport are well reproduced, while errors are reduced only moderately by a doubling of resolution. Model parameters are only weakly constrained by data, while strong correlations between mean error and parameter values are mostly found to be an artefact of single-parameter studies, not indicative of global model behaviour. Single-parameter sensitivity studies can therefore be misleading. Given a single, illustrative scenario of CO₂ increase and fixing the polynomial coefficients governing the extremely simple radiation parameterisation, the spread of model predictions for global mean warming due solely to the transport parameters is around one degree after 100 years forcing, although in a typical 4,000-year ensemble-member simulation, the peak rate of warming in the deep Pacific occurs 400 years after the onset of the forcing. The corresponding uncertainty in Atlantic overturning after 100 years is around 5 Sv, with a small, but non-negligible, probability of a collapse in the long term.