ALBERT

Description: Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model Geoscientific Model Development Discussions, 6, 4563-4601, 2013 Author(s): F. Cossu and K. Hocke This study examines how different microphysical parameterization schemes influence orographically-induced precipitation and the distributions of hydrometeors and water vapour for mid-latitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m −2 in the downwelling shortwave radiation and up to 33 W m −2 in the downwelling longwave radiation.

Description: AnaWEGE: a weather generator based on analogues of atmospheric circulation Geoscientific Model Development Discussions, 6, 4745-4774, 2013 Author(s): P. Yiou This paper presents a stochastic weather generator based on analogues of circulation (AnaWEGE). Analogues of circulation have been a promising paradigm to analyse climate variability and its extremes. The weather generator uses precomputed analogues of sea-level pressure over the North Atlantic. The stochastic rules of the generator constrain the continuity in time of the simulations. The generator then simulates spatially coherent time series of a climate variable, drawn from meteorological observations. The weather generator is tested for European temperatures, and for winter and summer seasons. The biases in temperature quantiles and autocorrelation are rather small compared to observed variability. The ability of simulating extremely hot summers and cold winters is also assessed.

Description: Methodological aspects of a pattern-scaling approach to produce global fields of monthly means of daily maximum and minimum temperature Geoscientific Model Development Discussions, 6, 4833-4882, 2013 Author(s): S. Kremser, G. E. Bodeker, and J. Lewis A Climate Pattern-Scaling Model (CPSM) that simulates global patterns of climate change, for a prescribed emissions scenario, is described. A CPSM works by quantitatively establishing the statistical relationship between a climate variable at a specific location (e.g. daily maximum surface temperature, T max ) and one or more predictor time series (e.g. global mean surface temperature, T global ) – referred to as the "training" of the CPSM. This training uses a regression model to derive fit-coefficients that describe the statistical relationship between the predictor time series and the target climate variable time series. Once that relationship has been determined, and given the predictor time series for any greenhouse gas (GHG) emissions scenario, the change in the climate variable of interest can be reconstructed – referred to as the "application" of the CPSM. The advantage of using a CPSM rather than a typical atmosphere-ocean global climate model (AOGCM) is that the predictor time series required by the CPSM can usually be generated quickly using a simple climate model (SCM) for any prescribed GHG emissions scenario and then applied to generate global fields of the climate variable of interest. The training can be performed either on historical measurements or on output from an AOGCM. Using model output from 21st century simulations has the advantage that the climate change signal is more pronounced than in historical data and therefore a more robust statistical relationship is obtained. The disadvantage of using AOGCM output is that the CPSM training might be compromised by any AOGCM inadequacies. For the purposes of exploring the various methodological aspects of the CPSM approach, AOGCM output was used in this study to train the CPSM. These investigations of the CPSM methodology focus on monthly mean fields of daily temperature extremes ( T max and T min ). Key conclusions are: (1) overall, the CPSM trained on simulations based on the Representative Concentration Pathway (RCP) 8.5 emissions scenario is able to reproduce AOGCM simulations of T max and T min based on predictor time series from an RCP 4.5 emissions scenario; (2) access to hemisphere average land and ocean temperatures as predictors improves the variance that can be explained, particularly over the oceans; (3) regression model fit-coefficients derived from individual simulations based on the RCP 2.6, 4.5 and 8.5 emissions scenarios agree well over most regions of the globe (the Arctic is the exception); (4) training the CPSM on concatenated time series from an ensemble of simulations does not result in fit-coefficients that explain significantly more of the variance than an approach that weights results based on single simulation fits; and (5) the inclusion of a linear time dependence in the regression model fit-coefficients improves the variance explained, primarily over the oceans.

Description: Modeling different freeze/thaw processes in heterogeneous landscapes of the Arctic polygonal tundra using an ecosystem model Geoscientific Model Development Discussions, 6, 4883-4932, 2013 Author(s): given_name prefix surname suffix, S. Yi, K. Wischnewski, M. Langer, S. Muster, and J. Boike Freeze/thaw (F/T) processes can be quite different under the various land surface types found in the heterogeneous polygonal tundra of the Arctic. Proper simulation of these different processes is essential for accurate prediction of the release of greenhouse gases under a warming climate scenario. In this study we have modified the dynamic organic soil version of the Terrestrial Ecosystem Model (DOS-TEM) to simulate F/T processes beneath the polygon rims, polygon centers (with and without water), and lakes that are common features in Arctic lowland regions. We first verified the F/T algorithm in the DOS-TEM against analytical solutions, and then compared the results with in situ measurements from Samoylov Island, Siberia. In the final stage, we examined the different responses of the F/T processes for different water levels at the various land surface types. The simulations revealed that (1) the DOS-TEM was very efficient and its results compared very well with analytical solutions for idealized cases, (2) the simulations compared reasonably well with in situ measurements although there were a number of model limitations and uncertainties, (3) the DOS-TEM was able to successfully simulate the differences in F/T dynamics under different land surface types, and (4) permafrost beneath water bodies was found to respond highly sensitive to changes in water depths between 1 and 2 m. Our results indicate that water is very important in the thermal processes simulated by the DOS-TEM; the heterogeneous nature of the landscape and different water depths therefore need to be taken into account when simulating methane emission responses to a warming climate.

Description: A standard test case suite for two-dimensional linear transport on the sphere: results from a collection of state-of-the-art schemes Geoscientific Model Development Discussions, 6, 4983-5076, 2013 Author(s): P. H. Lauritzen, P. A. Ullrich, C. Jablonowski, P. A. Bosler, D. Calhoun, A. J. Conley, T. Enomoto, L. Dong, S. Dubey, O. Guba, A. B. Hansen, E. Kaas, J. Kent, J.-F. Lamarque, M. J. Prather, D. Reinert, V. V. Shashkin, W. C. Skamarock, B. Sørensen, M. A. Taylor, and M. A. Tolstykh Recently, a standard test case suite for 2-D linear transport on the sphere was proposed to assess important aspects of accuracy in geophysical fluid dynamics with a "minimal" set of idealized model configurations/runs/diagnostics. Here we present results from 19 state-of-the-art transport scheme formulations based on finite-difference/finite-volume methods as well as emerging (in the context of atmospheric/oceanographic sciences) Galerkin methods. Discretization grids range from traditional regular latitude-longitude grids to more isotropic domain discretizations such as icosahedral and cubed-sphere tessellations of the sphere. The schemes are evaluated using a wide range of diagnostics in idealized flow environments. Accuracy is assessed in single- and two-tracer configurations using conventional error norms as well as novel diagnostics designed for climate and climate-chemistry applications. In addition, algorithmic considerations that may be important for computational efficiency are reported on. The latter is inevitably computing platform dependent, The ensemble of results from a wide variety of schemes presented here helps shed light on the ability of the test case suite diagnostics and flow settings to discriminate between algorithms and provide insights into accuracy in the context of global atmospheric/ocean modeling. A library of benchmark results is provided to facilitate scheme intercomparison and model development. Simple software and data-sets are made available to facilitate the process of model evaluation and scheme intercomparison.

Description: On the analytic approximation of bulk collision rates of non-spherical hydrometeors Geoscientific Model Development Discussions, 6, 5077-5116, 2013 Author(s): A. Seifert, U. Blahak, and R. Buhr Analytic approximations of the binary collision rates of hydrometeors are derived for use in bulk microphysical parameterizations. Special attention is given to non-spherical hydrometeors like raindrops and snowflakes. The terminal fall velocity of these particles cannot be sufficiently well approximated by power law relations which are used in most microphysical parameterizations and therefore an improved formulation is needed. The analytic approximations of the bulk collision rates given in this paper are an alternative to look-up tables and can replace the Wisner approximation which is used in many atmospheric models.

Description: Ensemble initialization of the oceanic component of a coupled model through bred vectors at seasonal-to-interannual time scales Geoscientific Model Development Discussions, 6, 5189-5214, 2013 Author(s): J. Baehr and R. Piontek We evaluate the ensemble spread at seasonal-to-interannual time scales for two perturbation techniques implemented into the ocean component of a coupled model: (1) lagged initial conditions as commonly used for decadal predictions, (2) bred vectors as commonly used for weather and seasonal forecasting. We show that relative to an uninitialized reference simulation the implementation for bred vectors can improve the ensemble spread compared to lagged initialization at time scales from one months up to three years. As bred vectors have so far mostly been used at short time scales, we initially focus on the implementation of the bred vectors into the ocean component. We introduce a depth-dependent vertical rescaling norm, accounting for the vertical dependence of the variability, and extending the commonly used upper-ocean rescaling norm to the full water column. We further show that it is sufficient for the (sub-surface) ocean to breed temperature and salinity (i.e., scalar quantities), and rely on the governing physics to carry the temperature and salinity perturbations to the flow field. Using these bred vectors with a rescaling interval of 12 months, we initialize hindcast simulations and compare them to hindcast simulations initialized with lagged initial conditions. We quantify the ensemble spread by analyzing Talagrand diagrams and spread-error ratios. For both temperature and salinity, the lagged initialized ensemble is particularly under-dispersive for the first few months of predictable lead time. The ensemble initialized with bred vectors improves the spread for temperature and salinity for the 0–700 m and 1000–3500 m means, compared to the lagged ensemble at lead times of several months to one year. As the lead time increases to years, the differences between the two ensemble initialization techniques becomes more difficult to discern. While the results need to be confirmed in an initialized framework, the present analysis represents a first step towards an improved ensemble generation at the transition from seasonal-to-interannual time scales, in particular at lead times up to one year.

Description: Understanding the performance of the FLake model over the African Great Lakes Geoscientific Model Development Discussions, 6, 5141-5187, 2013 Author(s): W. Thiery, A. Martynov, F. Darchambeau, J.-P. Descy, P.-D. Plisnier, L. Sushama, and N. P. M. van Lipzig The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu, Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding Automatic Weather Stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration allowed to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and temperature structure. At Lake Kivu, mixolimnion temperatures predicted by FLake were found sensitive both to minimal variations in the external parameters (lake depth and water transparency) as to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may already lead to a regime switch from the correctly represented seasonal mixed layer deepening to either completely mixed (down to the model lake bottom) or permanently stratified (from ~10 m downwards) conditions. In contrast, model temperatures are found robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterize tropical lake water surface temperatures within atmospheric prediction models, but may be less appropriate, in its current form, to study complex limnological processes within tropical lakes. Furthermore, a study of different initial conditions showed that for tropical lakes lacking reliable initial data, a fully mixed, artificially warm initialisation is to be preferred, but only if the model is allowed to spin up until convergence is reached. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60 m during the main dry season is primarily due to enhanced lake evaporation and secondarily due to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Description: The mid-Pliocene climate simulated by FGOALS-g2 Geoscientific Model Development Discussions, 6, 2403-2428, 2013 Author(s): W. Zheng, Z. Zhang, L. Chen, and Y. Yu Within the framework of Pliocene Model Intercomparison Project (PlioMIP), the mid-Pliocene (3.264–3.025 Ma) climate simulated by the Flexible Global Ocean-Atmosphere-Land System model grid-point version 2 (FGOALS-g2) are analyzed in this study. Results show that the model reproduces the large-scale features of the global warming over the land and ocean. The simulated mid-Pliocene global annual mean surface air temperature (TAS) and sea surface temperature (SST) are 4.17 and 2.62°C warmer than the pre-Industrial simulation, respectively. In particular, the feature of larger warming over mid-high latitudes is well captured. In the simulated warm mid-Pliocene climate, the Atlantic Meridional Overturning Circulation (AMOC) and El Niño-Southern Oscillation (ENSO) become weaker.

Description: The potential of an observational data set for calibration of a computationally expensive computer model Geoscientific Model Development Discussions, 6, 2369-2401, 2013 Author(s): D. J. McNeall, P. G. Challenor, J. R. Gattiker, and E. J. Stone We measure the potential of an observational data set to constrain a set of inputs to a complex and computationally expensive computer model. We use each member in turn of an ensemble of output from a computationally expensive model, corresponding to some observable part of a modelled system, as a proxy for an observational data set. We argue that our ability to constrain uncertain parameter inputs to a model using its own output as data, provides a maximum bound for our ability to constrain the model inputs using observations of the real system. The ensemble provides a set of known parameter input and model output pairs, which we use to build a computationally efficient statistical proxy for the full computer model, termed an emulator. We use the emulator to find and rule out ''implausible" values for the inputs of held-out ensemble members, given the computer model output. As we know the true values of the inputs for the ensemble, we can compare our constraint of the model inputs with the true value of the input for any ensemble member. Measures of the quality of constraint have the potential to inform strategy for data collection campaigns, before any real-world data is collected, as well as acting as an effective sensitivity analysis. We use an ensemble of the ice sheet model Glimmer to demonstrate our measures of quality of constraint. The ensemble has 250 model runs with 5 uncertain input parameters, and an output variable representing the pattern of the thickness of ice over Greenland. We have an observation of historical ice sheet thickness that directly matches the output variable, and offers an opportunity to constrain the model. We show that different ways of summarising our output variable (ice volume, ice surface area and maximum ice thickness) offer different potential constraints on individual input parameters. We show that combining the observational data gives increased power to constrain the model. We investigate the impact of uncertainty in observations or in model biases on our measures, showing that even a modest uncertainty can seriously degrade the potential of the observational data to constrain the model.

Description: On searching for optimized set of physical parameterization schemes in a multi-physics land surface process model Geoscientific Model Development Discussions, 6, 4511-4530, 2013 Author(s): S. Hong, X. Yu, S. K. Park, Y.-S. Choi, and B. Myoung Optimization of land surface models has been very challenging due to the increasing complexity of such models. Typical parameter calibration techniques often limit the solution of the spatiotemporal discrepancy in the modeling performance levels especially for regional applications. Thus, in this study, an attempt was made to perform scheme-based model optimization by designing a framework for coupling a micro-genetic algorithm (micro-GA) with the Noah land surface model that has multiple physics options (Noah-MP). Micro-GA controls the scheme selections in 10 different land surface parameterization fields in Noah-MP in order to extract the optimal scheme combination for a certain region. This coupling framework was successfully applied to the optimization of the surface water partitioning in the Korean Peninsula, promising not only the effectiveness of the scheme-based optimization but also model diagnosis capability by exploring the scheme sensitivity during the micro-GA evolution process. Then, the method was applied to four different regions in East Asia that have different climatic characteristics. The results indicate that (1) the optimal scheme combinations vary with the regions, (2) schemes related to the surface water partitioning are important for the modeling accuracy, and (3) specialized post-parameter optimization for each region may be required.

Description: Atmospheric inverse modeling with known physical bounds: an example from trace gas emissions Geoscientific Model Development Discussions, 6, 4531-4562, 2013 Author(s): S. M. Miller, A. M. Michalak, and P. J. Levi Many inverse problems in the atmospheric sciences involve parameters with known physical constraints. Examples include non-negativity (e.g., emissions of some urban air pollutants) or upward limits implied by reaction or solubility constants. However, probabilistic inverse modeling approaches based on Gaussian assumptions cannot incorporate such bounds and thus often produce unrealistic results. The atmospheric literature lacks consensus on the best means to overcome this problem, and existing atmospheric studies rely on a limited number of the possible methods with little examination of the relative merits of each. This paper investigates the applicability of several approaches to bounded inverse problems and is also the first application of Markov chain Monte Carlo (MCMC) to estimation of atmospheric trace gas fluxes. The approaches discussed here are broadly applicable. A common method of data transformations is found to unrealistically skew estimates for the examined example application. The method of Lagrange multipliers and two MCMC methods yield more realistic and accurate results. In general, the examined MCMC approaches produce the most realistic result but can require substantial computational time. Lagrange multipliers offer an appealing alternative for large, computationally intensive problems when exact uncertainty bounds are less central to the analysis. A synthetic data inversion of US anthropogenic methane emissions illustrates the strengths and weaknesses of each approach.

Description: A database and tool for boundary conditions for regional air quality modeling: description and evaluation Geoscientific Model Development Discussions, 6, 4665-4704, 2013 Author(s): B. H. Henderson, F. Akhtar, H. O. T. Pye, S. L. Napelenok, and W. T. Hutzell Transported air pollutants receive increasing attention as regulations tighten and global concentrations increase. The need to represent international transport in regional air quality assessments requires improved representation of boundary concentrations. Currently available observations are too sparse vertically to provide boundary information, particularly for ozone precursors, but global simulations can be used to generate spatially and temporally varying Lateral Boundary Conditions (LBC). This study presents a public database of global simulations designed and evaluated for use as LBC for air quality models (AQMs). The database covers the contiguous United States (CONUS) for the years 2000–2010 and contains hourly varying concentrations of ozone, aerosols, and their precursors. The database is complimented by a tool for configuring the global results as inputs to regional scale models (e.g., Community Multiscale Air Quality or Comprehensive Air quality Model with extensions). This study also presents an example application based on the CONUS domain, which is evaluated against satellite retrieved ozone vertical profiles. The results show performance is largely within uncertainty estimates for the Tropospheric Emission Spectrometer (TES) with some exceptions. The major difference shows a high bias in the upper troposphere along the southern boundary in January. This publication documents the global simulation database, the tool for conversion to LBC, and the fidelity of concentrations on the boundaries. This documentation is intended to support applications that require representation of long-range transport of air pollutants.

Description: Effects of vegetation structure on biomass accumulation in a Balanced Optimality Structure Vegetation Model (BOSVM v1.0) Geoscientific Model Development Discussions, 6, 4603-4663, 2013 Author(s): Z. Yin, S. C. Dekker, B. J. J. M. van den Hurk, and H. A. Dijkstra A myriad of interactions exist between vegetation and local climate for arid and semi-arid regions. Vegetation function, structure and individual behavior have large impacts on carbon-water-energy balances, which consequently influence local climate variability that, in turn, feeds back to the vegetation. In this study, a conceptual vegetation structure scheme is formulated and tested in a new carbon-water-energy coupled model to explore the importance of vegetation structure and vegetation adaptation to water stress on equilibrium biomass states. Surface energy, water and carbon fluxes are simulated for a range of vegetation structures across a precipitation gradient in West Africa and optimal vegetation structures that maximizes biomass for each precipitation regime are determined. Two different strategies of vegetation adaptation to water stress are included. Under dry conditions vegetation tries to maximize the Water Use Efficiency and Leaf Area Index as it tries to maximize carbon gain. However, an important negative feedback mechanism is found as the vegetation also tries to minimize its cover to optimize the surrounding bare ground area from which water can be extracted, thereby forming patches of vertical vegetation. Under larger precipitation, a positive feedback mechanism is found in which vegetation tries to maximize its cover as it then can reduce water loss from bare soil while having maximum carbon gain due to a large Leaf Area Index. The competition between vegetation and bare soil determines a transition between a "survival" state to a "growing" state.

Description: A distributed computing approach to improve the performance of the Parallel Ocean Program (v2.1) Geoscientific Model Development Discussions, 6, 4705-4744, 2013 Author(s): B. van Werkhoven, J. Maassen, M. Kliphuis, H. A. Dijkstra, S. E. Brunnabend, M. van Meersbergen, F. J. Seinstra, and H. E. Bal The Parallel Ocean Program (POP) is used in many strongly eddying ocean circulation simulations. Ideally one would like to do thousand-year long simulations, but the current performance of POP prohibits this type of simulations. In this work, using a new distributed computing approach, two innovations to improve the performance of POP are presented. The first is a new block partitioning scheme for the optimization of the load balancing of POP such that it can be run efficiently in a multi-platform setting. The second is an implementation of part of the POP model code on Graphics Processing Units. We show that the combination of both innovations leads to a substantial performance increase also when running POP simultaneously over multiple computational platforms.

Description: Inherently mass-conservative version of the semi-Lagrangian Absolute Vorticity (SL-AV) atmospheric model dynamical core Geoscientific Model Development Discussions, 6, 4809-4832, 2013 Author(s): V. V. Shashkin and M. A. Tolstykh The semi-Lagrangian Absolute Vorticity (SL-AV) atmospheric model is the global semi-Lagrangian hydrostatic model used for operational medium-range and seasonal forecasts at Hydrometeorological centre of Russia. The distinct feature of SL-AV dynamical core is the semi-implicit semi-Lagrangian vorticity-divergence formulation on the unstaggered grid. Semi-implicit semi-Lagrangian approach allows for long time steps while violates the global and local mass-conservation. In particular, the total mass in simulations with semi-Lagrangian models can drift significantly if no aposteriori mass-fixing algorithms are applied. However, the global mass-fixing algorithms degrade the local mass conservation. The inherently mass-conservative version of SL-AV model dynamical core presented in the article ensures global and local mass conservation without mass-fixing algorithms. The mass conservation is achieved with the introduction of the finite-volume semi-Lagrangian discretization for continuity equation based on the 3-D extension of the conservative cascade semi-Lagrangian transport scheme (CCS). The numerical experiments show that the presented new version of SL-AV dynamical core combines the accuracy and stability of the standard SL-AV dynamical core with the mass-conservation properties. The results of the mountain induced Rossby wave test and baroclinic instability test for mass-conservative dynamical core are found to be in agreement with the results available in literature.

Description: A fast input/output library for high resolution climate models Geoscientific Model Development Discussions, 6, 4775-4807, 2013 Author(s): X. Huang, W. Wang, H. Fu, G. Yang, B. Wang, and C. Zhang We describe the design and implementation of Climate Fast Input/Output (CFIO), a fast input/output (I/O) library for high resolution climate models. CFIO provides a simple method for modelers to overlap the I/O phase with the computing phase automatically, so as to shorten the running time of numerical simulations. To minimize the code modifications required for porting, CFIO provides similar interfaces and features to Parallel network Common Data Form (PnetCDF), which is one of the most widely used I/O libraries in climate models. We deployed CFIO in three high resolution climate models, including two ocean models (POP and LICOM) and one sea ice model (CICE). The experimental results show that CFIO improves the performance of climate models significantly versus the original serial I/O approach. When running with CFIO at 0.1° resolution with about 1000 CPU cores, we managed to reduce the running time by factors of 7.9, 4.6 and 2.0 for POP, CICE, and LICOM respectively. We also compared the performance of CFIO against PnetCDF in different scenarios. For scenarios with both data output operations and computations, CFIO decreases the I/O overhead by a factor of 5.1 compared to PnetCDF.

Description: Application and evaluation of McICA scheme with new radiation code in BCC_AGCM2.0.1 Geoscientific Model Development Discussions, 6, 4933-4982, 2013 Author(s): given_name prefix surname suffix, H. Zhang, X. Jing, and J. Li This research incorporates the Monte Carlo Independent Column Approximation (McICA) scheme with the correlated k-distribution BCC-RAD radiation model into the climate model BCC_AGCM2.0.1 and examines the impacts on modeled climate through several simulations with variations in cloud structures. Results from experiments with consistent sub-grid cloud structures show that both clear-sky radiation fluxes and cloud radiative forcings (CRFs) calculated by the new scheme are mostly improved relative to those calculated from the original one. The modeled atmospheric temperature and specific humidity are also improved due to changes in the radiative heating rates. The vertical overlap of fractional clouds and horizontal distribution of cloud condensation are important for computing CRFs. The maximum changes in seasonal CRF using the general overlap assumption (GenO) with different decorrelation depths ( L cf ) are larger than 10 and 20 Wm 2 for longwave (LW) CRF and shortwave (SW) CRF, respectively, mostly located in the Tropics and mid-latitude storm tracks. Larger (smaller) L cf in the Tropics (mid-latitude storm tracks) yield better cloud fraction and CRF compared with observations. The inclusion of an observation-based horizontal inhomogeneity of cloud condensation has a distinct impact on LW CRF and SW CRF, with global means of ∼1.2 Wm −2 and ∼3.7 Wm −2 at the top of atmosphere, respectively, making these much closer to observations. These results prove the reliability of the new model configuration to be used in BCC_AGCM2.0.1 for climate simulations, and also indicate that more detailed real-world information on cloud structures should be obtained to constrain cloud settings in McICA in the future.

Description: A regional climate modelling projection ensemble experiment – NARCliM Geoscientific Model Development Discussions, 6, 5117-5139, 2013 Author(s): J. P. Evans, F. Ji, C. Lee, P. Smith, D. Argüeso, and L. Fita Including the impacts of climate change in decision making and planning processes is a challenge facing many regional governments including the New South Wales (NSW) and Australian Capital Territory (ACT) governments in Australia. NARCliM (NSW/ACT Regional Climate Modelling project) is a regional climate modelling project that aims to provide a comprehensive and consistent set of climate projections that can be used by all relevant government departments when considering climate change. To maximise end user engagement and ensure outputs are relevant to the planning process, a series of stakeholder workshops were run to define key aspects of the model experiment including spatial resolution, time slices, and output variables. As with all such experiments, practical considerations limit the number of ensembles members that can be simulated such that choices must be made concerning which Global Climate Models (GCMs) to downscale from, and which Regional Climate Models (RCMs) to downscale with. Here a methodology for making these choices is proposed that aims to sample the uncertainty in both GCMs and RCMs, as well as spanning the range of future climate projections present in the full GCM ensemble. The created ensemble provides a more robust view of future regional climate changes.

Description: Efficient performance of the Met Office Unified Model v8.2 on Intel Xeon partially used nodes Geoscientific Model Development Discussions, 7, 7395-7425, 2014 Author(s): I. Bermous The atmospheric Unified Model (UM) developed at the UK Met Office is used for weather and climate prediction by forecast teams at a number of international meteorological centres and research institutes on a wide variety of hardware and software environments. Over its 25 year history the UM sources have been optimised for a better application performance on a number of HPC systems including NEC SX vector architecture systems and recently the IBM Power6/Power7 platforms. Understanding the influence of the compiler flags, MPI libraries and run configurations is crucial to achieving the shortest elapsed times for a UM application on any particular HPC system. These aspects are very important for applications that must run within operational time frames. Driving the current study is the HPC industry trend since 1980 for processor arithmetic performance to increase at a faster rate than memory bandwidth. This gap has been growing especially fast for multicore processors in the past 10 years and it can have significant implication for the performance and performance scaling of memory bandwidth intensive applications, such as the UM. Analysis of partially used nodes on Intel Xeon clusters is provided in this paper for short and medium range weather forecasting systems using global and limited-area configurations. It is shown that on the Intel Xeon based clusters the fastest elapsed times and the most efficient system usage can be achieved using partially committed nodes.

Description: Simulation of groundwater and surface water over the continental US using a hyperresolution, integrated hydrologic model Geoscientific Model Development Discussions, 7, 7317-7349, 2014 Author(s): R. M. Maxwell, L. E. Condon, and S. J. Kollet Interactions between surface and groundwater systems are well-established theoretically and observationally. While numerical models that solve both surface and subsurface flow equations in a single framework (matrix) are increasingly being applied, computational limitations have restricted their use to local and regional studies. Regional or watershed, scale simulations have been effective tools in understanding hydrologic processes, however there are still many questions, such as the adaptation of water resources to anthropogenic stressors and climate variability, that need to be answered across large spatial extents at high resolution. In response to this "grand challenge" in hydrology, we present the results of a parallel, integrated hydrologic model simulating surface and subsurface flow at high spatial resolution (1 km) over much of continental North America (~ 6 300 000 or 6.3 million km 2 ). These simulations provide predictions of hydrologic states and fluxes, namely water table depth and streamflow, at unprecedented scale and resolution. The physically-based modeling approach used here requires limited parameterizations and relies only on more fundamental inputs, such as topography, hydrogeologic properties and climate forcing. Results are compared to observations and provide mechanistic insight into hydrologic process interaction. This study demonstrates both the feasibility of continental scale integrated models and their utility for improving our understanding of large-scale hydrologic systems; the combination of high resolution and large spatial extent facilitates novel analysis of scaling relationships using model outputs.

Description: EDDA: integrated simulation of debris flow erosion, deposition and property changes Geoscientific Model Development Discussions, 7, 7267-7316, 2014 Author(s): H. X. Chen and L. M. Zhang Debris flow material properties change during the initiation, transportation and deposition processes, which influences the runout characteristics of the debris flow. A quasi-three-dimensional depth-integrated numerical model, EDDA, is presented in this paper to simulate debris flow erosion, deposition and induced material property changes. The model considers changes in debris flow density, yield stress and dynamic viscosity during the flow process. The yield stress of debris flow mixture is determined at limit equilibrium using the Mohr–Coulomb equation, which is applicable to clear water flow, hyper-concentrated flow and fully developed debris flow. To assure numerical stability and computational efficiency at the same time, a variable time stepping algorithm is developed to solve the governing differential equations. Four numerical tests are conducted to validate the model. The first two tests involve a one-dimensional dam-break water flow and a one-dimensional debris flow with constant properties. The last two tests involve erosion and deposition, and the movement of multi-directional debris flows. The changes in debris flow mass and properties due to either erosion or deposition are shown to affect the runout characteristics significantly. The model is also applied to simulate a large-scale debris flow in Xiaojiagou Ravine to test the performance of the model in catchment-scale simulations. The results suggest that the model estimates well the volume, inundated area, and runout distance of the debris flow. The model is intended for use as a module in a real-time debris flow warning system.

Description: Northern Hemisphere storminess in the Norwegian Earth System Model (NorESM1-M) Geoscientific Model Development Discussions, 7, 8975-9015, 2014 Author(s): E. M. Knudsen and J. E. Walsh Metrics of storm activity in Northern Hemisphere high- and midlatitudes are evaluated from historical output and future projections by the Norwegian Earth System Model (NorESM1-M) coupled global climate model. The European Re-Analysis Interim (ERA-Interim) and the Community Climate System Model (CCSM4), a global climate model of the same vintage as NorESM1-M, provide benchmarks for comparison. The focus is on the autumn and early winter (September through December), the period when the ongoing and projected Arctic sea ice retreat is greatest. Storm tracks derived from a vorticity-based algorithm for storm identification are reproduced well by NorESM1-M, although the tracks are somewhat better resolved in the higher-resolution ERA-Interim and CCSM4. The tracks are projected to shift polewards in the future as climate changes under the Representative Concentration Pathway (RCP) forcing scenarios. Cyclones are projected to become generally more intense in the high-latitudes, especially over the Alaskan region, although in some other areas the intensity is projected to decrease. While projected changes in track density are less coherent, there is a general tendency towards less frequent storms in midlatitudes and more frequent storms in high-latitudes, especially the Baffin Bay/Davis Strait region. Autumn precipitation is projected to increase significantly across the entire high-latitudes. Together with the projected increases in storm intensity and sea level and the loss of sea ice, this increase in precipitation implies a greater vulnerability to coastal flooding and erosion, especially in the Alaskan region. The projected changes in storm intensity and precipitation (as well as sea ice and sea level pressure) scale generally linearly with the RCP value of the forcing and with time through the 21st century.

Description: The 1-way on-line coupled atmospheric chemistry model system MECO(n) – Part 1: The limited-area atmospheric chemistry model COSMO/MESSy Geoscientific Model Development Discussions, 4, 1305-1358, 2011 Author(s): A. Kerkweg and P. Jöckel The numerical weather prediction model of the Consortium for Small Scale Modelling (COSMO), maintained by the German weather service (DWD), is connected with the Modular Earth Submodel System (MESSy). This effort is undertaken in preparation of a~new, limited-area atmospheric chemistry model. This model is as consistent as possible, with respect to atmospheric chemistry and related processes, with a previously developed global atmospheric chemistry general circulation model: the ECHAM/MESSy Atmospheric Chemistry (EMAC) model. The combined system constitutes a new research tool, bridging the global to the meso-γ scale for atmospheric chemistry research. MESSy provides the infrastructure and includes, among others, the process and diagnostic submodels for atmospheric chemistry simulations. Furthermore, MESSy is highly flexible allowing model setups with tailor made complexity, depending on the scientific question. Here, the connection of the MESSy infrastructure to the COSMO model is documented. Previously published prototype submodels for simplified tracer studies are generalised to be plugged-in and used in the global and the limited-area model. They are used to evaluate the tracer transport characteristics of the new COSMO/MESSy model system, an important prerequisite for future atmospheric chemistry applications. A supplementary document with further details on the technical implementation of the MESSy interface into COSMO with a complete list of modifications to the COSMO code is provided.

Description: Carbon monoxide as a tracer for tropical troposphere to stratosphere transport in the Chemical Lagrangian Model of the Stratosphere (CLaMS) Geoscientific Model Development Discussions, 4, 1185-1211, 2011 Author(s): R. Pommrich, R. Müller, J.-U. Grooß, P. Konopka, G. Günther, H.-C. Pumphrey, S. Viciani, F. D'Amato, and M. Riese Variations in the mixing ratio of trace gases of tropospheric origin entering the stratosphere in the tropics are of interest for assessing both troposphere to stratosphere transport fluxes in the tropics and the impact on the composition of the tropical lower stratosphere of quasi-horizontal in-mixing into the tropical tropopause layer from the mid-latitude stratosphere. Here, we present a simplified chemistry scheme for the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the simulation, at comparatively low numerical cost, of CO, ozone, and long-lived trace substances (CH 4 , N 2 O, CCl 3 F, and CO 2 ) in the lower tropical stratosphere. The boundary conditions at the ground are represented for the long-lived trace substances CH 4 , N 2 O, CCl 3 F, and CO 2 based on ground-based measurements. The boundary condition for CO in the free troposphere is deduced from MOPITT measurements. We find that the zonally averaged tropical CO anomaly patterns simulated by this model version of CLaMS are in good agreement with observations. The introduction of a new scheme in the ECMWF integrated forecast system (Tompkins et al., 2007) for the ice supersaturation after September 2006, results in a somewhat less good agreement between observed and simulated CO patterns in the tropical lower stratosphere after this date.

Description: The Met Office Unified Model Global Atmosphere 3.0/3.1 and JULES Global Land 3.0/3.1 configurations Geoscientific Model Development Discussions, 4, 1213-1271, 2011 Author(s): D. N. Walters, M. J. Best, A. C. Bushell, D. Copsey, J. M. Edwards, P. D. Falloon, C. M. Harris, A. P. Lock, J. C. Manners, C. J. Morcrette, M. J. Roberts, R. A. Stratton, S. Webster, J. M. Wilkinson, M. R. Willett, I. A. Boutle, P. D. Earnshaw, P. G. Hill, C. MacLachlan, G. M. Martin, W. Moufouma-Okia, M. D. Palmer, J. C. Petch, G. G. Rooney, A. A. Scaife, and K. D. Williams We describe Global Atmosphere 3.0 (GA3.0): a configuration of the Met Office Unified Model (MetUM) developed for use across climate research and weather prediction activities. GA3.0 has been formulated by converging the development paths of the Met Office's weather and climate global atmospheric model components such that wherever possible, atmospheric processes are modelled or parametrized seamlessly across spatial resolutions and timescales. This unified development process will provide the Met Office and its collaborators with regular releases of a configuration that has been evaluated, and can hence be applied, over a variety of modelling régimes. We also describe Global Land 3.0 (GL3.0): a configuration of the JULES community land surface model developed for use with GA3.0. This paper provides a comprehensive technical and scientific description of the GA3.0 and GL3.0 (and related GA3.1 and GL3.1) configurations and presents the results of some initial evaluations of their performance in various applications. It is to be the first in a series of papers describing each subsequent Global Atmosphere release; this will provide a single source of reference for established users and developers as well as researchers requiring access to a current, but trusted, global MetUM setup.

Description: FAMOUS, faster: using parallel computing techniques to accelerate the FAMOUS/HadCM3 climate model with a focus on the radiative transfer algorithm Geoscientific Model Development Discussions, 4, 1273-1303, 2011 Author(s): P. Hanappe, A. Beurivé, F. Laguzet, L. Steels, N. Bellouin, O. Boucher, Y. H. Yamazaki, T. Aina, and M. Allen We have optimised the atmospheric radiation algorithm of the FAMOUS climate model on several hardware platforms. The optimisation involved translating the Fortran code to C and restructuring the algorithm around the computation of a single air column. A task queue and a thread pool are used to distribute the computation to several processors. Finally, four air columns are packed together in a single data structure and computed simultaneously using Single Instruction Multiple Data operations. The modified algorithm runs more than 50 times faster on the CELL's Synergistic Processing Elements than on its main PowerPC processing element. On Intel-compatible processors, the new radiation code runs 4 times faster and on graphics processors, using OpenCL, more than 2.5 times faster, as compared to the original code. Because the radiation code takes more than 60 % of the total CPU time, FAMOUS executes more than twice as fast. Our version of the algorithm returns bit-wise identical results, which demonstrates the robustness of our approach.

Description: Analyzing numerics of bulk microphysics schemes in Community models: warm rain processes Geoscientific Model Development Discussions, 4, 1403-1434, 2011 Author(s): I. Sednev and S. Menon In the last decade there has been only one study that discussed time integration scheme (TIS) applied to advance governing differential equations in bulk microphysics (BLK) schemes. Recently, Morrison and Gettelman (2008) examine numerical aspects of double-moment BLK scheme with diagnostic treatment of precipitating hydrometeors implemented into Community Atmosphere Model, version 3 (CAM) to find an acceptable level of accuracy and numerical stability. However, stability condition for their explicit non-positive definite TIS was not defined. It is conventionally thought that the Weather Research and Forecasting (WRF) model can be applied for a broad range of spatial scales from large eddy up to global scale simulations if time steps used for model integration satisfy to a certain limit imposed mainly by dynamics. However, numerics used in WRF BLK schemes has never been analyzed in detail. To improve creditability of BLK schemes we derive a general analytical stability and positive definiteness criteria for explicit Eulerian time integration scheme used to advanced finite-difference equations that govern warm rain formation processes in microphysics packages in Community models (CAM and WRF) and define well-behaved, conditionally well-behaved, and non-well-behaved Explicit Eulerian Bulk Microphysics Code (EEBMPC) classes. We highlight that source codes of BLK schemes, originally developed for use in cloud-resolving models, implemented in Community models belong to conditionally well-behaved EEBMPC class and exhibit better performance for finer spatial resolutions when time steps do not exceed seconds or tenths of seconds. For coarser spatial resolutions used in regional and global scale simulations time steps are usually increased from hundredths up to thousands of seconds. This might lead to a degradation of conditionally well-behaved EEBMPCs ability to calculate the amount of precipitation as well as its spatial and temporal distribution since both stability and positive definiteness conditions are not met in the TIS. The correction through the so called “mass conservation” technique commonly used in many models with bulk microphysics is a main characteristic of non-well-behaved EEBMPC, whose utilization leads to erroneous conclusions regarding relative importance of different microphysical processes. Moreover, surface boundary conditions for ocean, land, lake, and sea ice models are dependent on the precipitation and its spatial and temporal distribution. Uncertainties in calculations of temporal and spatial patterns of accumulated precipitation influence the global water cycle. In fact, numerics in non-well-behaved EEBMPCs, which are used in Community Earth System Model, act as a hidden climate forcing agent, if relatively long time steps are used for the host model integration. By analyzing numerics of warm rain processes in EEBMPCs implemented in Community models we provide general guidelines regarding appropriate choice of integration time steps for use in these models.

Description: The 1-way on-line coupled atmospheric chemistry model system MECO(n) – Part 2: On-line coupling Geoscientific Model Development Discussions, 4, 1359-1402, 2011 Author(s): A. Kerkweg and P. Jöckel A new, highly flexible model system for the seamless dynamical down-scaling of meteorological and chemical processes from the global to the meso-γ scale is presented. A global model and a cascade of an arbitrary number of limited-area model instances run concurrently in the same parallel environment, in which the coarser grained instances provide the boundary data for the finer grained instances. Thus, disk-space intensive and time consuming intermediate and pre-processing steps are entirely avoided and the time interpolation errors of common off-line nesting approaches are minimised. More specifically, the regional model COSMO of the German Weather Service (DWD) is nested on-line into the atmospheric general circulation model ECHAM5 within the Modular Earth Submodel System (MESSy) framework. ECHAM5 and COSMO have previously been equipped with the MESSy infrastructure, implying that the same process formulations (MESSy submodels) are available for both models. This guarantees the highest degree of achievable consistency, between both, the meteorological and chemical conditions at the domain boundaries of the nested limited-area model, and between the process formulations on all scales. The on-line nesting of the different models is established by a client-server approach with the newly developed Multi-Model-Driver (MMD), an additional component of the MESSy infrastructure. With MMD an arbitrary number of model instances can be run concurrently within the same message passing interface (MPI) environment, the respective coarser model (either global or regional) is the server for the nested finer (regional) client model, i.e., it provides the data required to calculate the initial and boundary fields to the client model. On-line nesting means that the coupled (client-server) models exchange their data via the computer memory, in contrast to the data exchange via files on disk in common off-line nesting approaches. MMD consists of a library (Fortran95 and some parts in C) which is based on the MPI standard and two new MESSy submodels, MMDSERV and MMDCLNT (both Fortran95) for the server and client models, respectively. MMDCLNT contains a further sub-submodel, INT2COSMO, for the interpolation of the coarse grid data provided by the server models (either ECHAM5/MESSy or COSMO/MESSy) to the grid of the respective client model (COSMO/MESSy). INT2COSMO is based on the off-line pre-processing tool INT2LM provided by the DWD. The new achievements allow the setup of model cascades for zooming (down-scaling) from the global scale to the lower edge of the meso-γ scale (≈1 km) with a very high degree of consistency between the different models and between the chemical and meteorological boundary conditions.

Description: Development and evaluation of an Earth-system model – HadGEM2 Geoscientific Model Development Discussions, 4, 997-1062, 2011 Author(s): W. J. Collins, N. Bellouin, M. Doutriaux-Boucher, N. Gedney, P. Halloran, T. Hinton, J. Hughes, C. D. Jones, M. Joshi, S. Liddicoat, G. Martin, F. O'Connor, J. Rae, C. Senior, S. Sitch, I. Totterdell, A. Wiltshire, and S. Woodward We describe here the development and evaluation of an Earth system model suitable for centennial-scale climate prediction. The principal new components added to the physical climate model are the terrestrial and ocean ecosystems and gas-phase tropospheric chemistry, along with their coupled interactions. The individual Earth system components are described briefly and the relevant interactions between the components are explained. Because the multiple interactions could lead to unstable feedbacks, we go through a careful process of model spin up to ensure that all components are stable and the interactions balanced. This spun-up configuration is evaluated against observed data for the Earth system components and is generally found to perform very satisfactorily. The reason for the evaluation phase is that the model is to be used for the core climate simulations carried out by the Met Office Hadley Centre for the Coupled Model Intercomparison Project (CMIP5), so it is essential that addition of the extra complexity does not detract substantially from its climate performance. Localised changes in some specific meteorological variables can be identified, but the impacts on the overall simulation of present day climate are slight. This model is proving valuable both for climate predictions, and for investigating the strengths of biogeochemical feedbacks.

Description: A multi-layer land surface energy budget model for implicit coupling with global atmospheric simulations Geoscientific Model Development Discussions, 7, 8649-8701, 2014 Author(s): J. Ryder, J. Polcher, P. Peylin, C. Ottlé, Y. Chen, E. van Gorsel, V. Haverd, M. J. McGrath, K. Naudts, J. Otto, A. Valade, and S. Luyssaert In Earth system modelling, a description of the energy budget of the vegetated surface layer is fundamental as it determines the meteorological conditions in the planetary boundary layer and as such contributes to the atmospheric conditions and its circulation. The energy budget in most Earth system models has long been based on a "big-leaf approach", with averaging schemes that represent in-canopy processes. Such models have difficulties in reproducing consistently the energy balance in field observations. We here outline a newly developed numerical model for energy budget simulation, as a component of the land surface model ORCHIDEE-CAN (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy). This new model implements techniques from single-site canopy models in a practical way. It includes representation of in-canopy transport, a multilayer longwave radiation budget, height-specific calculation of aerodynamic and stomatal conductance, and interaction with the bare soil flux within the canopy space. Significantly, it avoids iterations over the height of tha canopy and so maintains implicit coupling to the atmospheric model LMDz. As a first test, the model is evaluated against data from both an intensive measurement campaign and longer term eddy covariance measurements for the intensively studied Eucalyptus stand at Tumbarumba, Australia. The model performs well in replicating both diurnal and annual cycles of fluxes, as well as the gradients of sensible heat fluxes. However, the model overestimates sensible heat flux against an underestimate of the radiation budget. Improved performance is expected through the implementation of a more detailed calculation of stand albedo and a more up-to-date stomatal conductance calculation.

Description: iGen: the automated generation of a parameterisation of entrainment in marine stratocumulus Geoscientific Model Development Discussions, 4, 971-995, 2011 Author(s): D. F. Tang and S. Dobbie In a previous paper we described a new technique for automatically generating parameterisations using a program called iGen. iGen generates parameterisations by analysing the source code of a high resolution model that resolves the physics to be parameterised. In order to demonstrate that this technique scales up to deal with models of realistic complexity we have used iGen to generate a parameterisation of entrainment in marine stratocumulus. We present details of our technique in which iGen was used to analyse the source code of a cloud resolving model and generate a parameterisation of the mean and standard deviation of entrainment velocity in marine stratocumulus in terms of the large-scale state of the boundary layer. The parameterisation was tested against results from the DYCOMS-II intercomparison of cloud resolving models and iGen's parameterisation of mean entrainment velocity was found to be 5.27 × 10 −3 ± 0.62 × 10 −3 m s −1 compared to 5.2 × 10 −3 ± 0.8 × 10 −3 m s −1 for the DYCOMS-II ensemble of cloud resolving models.

Description: The HadGEM2-ES implementation of CMIP5 centennial simulations Geoscientific Model Development Discussions, 4, 689-763, 2011 Author(s): C. D. Jones, J. K. Hughes, N. Bellouin, S. C. Hardiman, G. S. Jones, J. Knight, S. Liddicoat, F. M. O'Connor, R. J. Andres, C. Bell, K.-O. Boo, A. Bozzo, N. Butchart, P. Cadule, K. D. Corbin, M. Doutriaux-Boucher, P. Friedlingstein, J. Gornall, L. Gray, P. R. Halloran, G. Hurtt, W. Ingram, J.-F. Lamarque, R. M. Law, M. Meinshausen, S. Osprey, E. J. Palin, L. Parsons Chini, T. Raddatz, M. Sanderson, A. A. Sellar, A. Schurer, P. Valdes, N. Wood, S. Woodward, M. Yoshioka, and M. Zerroukat The scientific understanding of the Earth's climate system, including the central question of how the climate system is likely to respond to human-induced perturbations, is comprehensively captured in GCMs and Earth System Models(ESM). Diagnosing the simulated climate response, and comparing responses across different models, is crucially dependent on transparent assumptions of how the GCM/ESM has been driven – especially because the implementation can involve subjective decisions and may differ between modelling groups performing the same experiment. This paper outlines the climate forcings and setup of the Met Office Hadley Centre ESM, HadGEM2-ES for the CMIP5 set of centennial experiments. We document the prescribed greenhouse gas concentrations, aerosol precursors, stratospheric and tropospheric ozone assumptions, as well as implementation of land-use change and natural forcings for the HadGEM2-ES historical and future experiments following the Representative Concentration Pathways. In addition, we provide details of how HadGEM2-ES ensemble members were initialised from the control run and how the palaeoclimate and AMIP experiments, as well as the "emission-driven" RCP experiments were performed.

Description: The HadGEM2 family of Met Office Unified Model Climate configurations Geoscientific Model Development Discussions, 4, 765-841, 2011 Author(s): The HadGEM2 Development Team: Martin, G. M., N. Bellouin, W. J. Collins, I. D. Culverwell, P. R. Halloran, S. C. Hardiman, T. J. Hinton, C. D. Jones, R. E. McDonald, A. J. McLaren, F. M. O'Connor, M. J. Roberts, J. M. Rodriguez, S. Woodward, M. J. Best, M. E. Brooks, A. R. Brown, N. Butchart, C. Dearden, S. H. Derbyshire, I. Dharssi, M. Doutriaux-Boucher, J. M. Edwards, P. D. Falloon, N. Gedney, L. J. Gray, H. T. Hewitt, M. Hobson, M. R. Huddleston, J. Hughes, S. Ineson, W. J. Ingram, P. M. James, T. C. Johns, C. E. Johnson, A. Jones, C. P. Jones, M. M. Joshi, A. B. Keen, S. Liddicoat, A. P. Lock, A. V. Maidens, J. C. Manners, S. F. Milton, J. G. L. Rae, J. K. Ridley, A. Sellar, C. A. Senior, I. J. Totterdell, A. Verhoef, P. L. Vidale, and A. Wiltshire We describe the HadGEM2 family of climate configurations of the Met Office Unified Model, MetUM. The concept of a model "family" comprises a range of specific model configurations incorporating different levels of complexity but with a common physical framework. The HadGEM2 family of configurations includes atmosphere and ocean components, with and without a vertical extension to include a well-resolved stratosphere, and an Earth-System (ES) component which includes dynamic vegetation, ocean biology and atmospheric chemistry. The HadGEM2 physical model includes improvements designed to address specific systematic errors encountered in the previous climate configuration, HadGEM1, namely Northern Hemisphere continental temperature biases and tropical sea surface temperature biases and poor variability. Targeting these biases was crucial in order that the ES configuration could represent important biogeochemical climate feedbacks. Detailed descriptions and evaluations of particular HadGEM2 family members are included in a number of other publications, and the discussion here is limited to a summary of the overall performance using a set of model metrics which compare the way in which the various configurations simulate present-day climate and its variability.

Description: Towards an online-coupled chemistry-climate model: evaluation of COSMO-ART Geoscientific Model Development Discussions, 4, 1809-1874, 2011 Author(s): C. Knote, D. Brunner, H. Vogel, J. Allan, A. Asmi, M. Äijälä, S. Carbone, H. D. van der Gon, J. L. Jimenez, A. Kiendler-Scharr, C. Mohr, L. Poulain, A. S. H. Prévôt, E. Swietlicki, and B. Vogel The online-coupled, regional chemistry transport model COSMO-ART is evaluated for periods in all seasons against several measurement datasets to assess its ability to represent gaseous pollutants and ambient aerosol characteristics over the European domain. Measurements used in the comparison include long-term station observations, satellite and ground-based remote sensing products, and complex datasets of aerosol chemical composition and number size distribution from recent field campaigns. This is the first time these comprehensive measurements of aerosol characteristics in Europe are used to evaluate a regional chemistry transport model. We show a detailed analysis of the simulated size-resolved chemical composition under different meteorological conditions. The model is able to represent trace gas concentrations with good accuracy and reproduces bulk aerosol properties rather well though with a clear tendency to underestimate both total mass (PM 10 and PM 2.5 ) and aerosol optical depth. We find indications of an overestimation of shipping emissions. Time evolution of aerosol chemical composition is captured, although some biases are found in relative composition. Nitrate aerosol components are on average overestimated, and sulfates underestimated. The accuracy of simulated organics depends strongly on season and location. While strongly underestimated during summer, organic mass is comparable in spring and autumn. We see indications for an overestimated fractional contribution of primary organic matter in urban areas and an underestimation of SOA at many locations. Aerosol number concentrations can be simulated well, size distributions are comparable. Our work sets the basis for subsequent studies of aerosol characteristics and climate impacts with COSMO-ART, and highlights areas where improvements are necessary for current regional modeling systems in general.

Description: Mapping the climate: guidance on appropriate techniques to map climate variables and their uncertainty Geoscientific Model Development Discussions, 4, 1875-1906, 2011 Author(s): N. R. Kaye, A. Hartley, and D. Hemming Maps are a crucial asset in communicating climate science to a diverse audience, and there is a wealth of software available to analyse and visualise climate information. However, this availability makes it easy to create poor maps as users often lack an underlying cartographic knowledge. Furthermore, communicating and visualising uncertainties in climate data and climate change projections, using for example ensemble based approaches, presents additional challenges for mapping that require careful consideration. This paper assesses a range of techniques for mapping uncertainties, comparing "intrinsic" approaches that use colour in much the same way as conventional thematic maps, and "extrinsic" approaches that incorporate additional geometry such as points or features. We proposes that, unlike traditional cartography, where many known standards allow maps to be interpreted easily, there is no standard mapping approach used to represent uncertainty (in climate or other information). Consequently, a wide range of techniques have been applied for this purpose, and users may spend unnecessary time trying to understand the mapping approach rather than interpreting the information presented. We use cartographic knowledge and lessons learned from mapping other information to propose a suitable mapping technique that represents both magnitude and uncertainty in climate data. This technique adjusts the hue of a small palette of colours to show the mean or median of a climate variable, and the saturation of the colour to illustrate a measure of uncertainty. It is designed to be easy to replicate, visible to colour blind people and intuitive to understand. This technique may be utilised to map a wide range of climate data, and it is proposed that it would provide a consistent approach suitable for mapping information for the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5).

Description: A new version of the CNRM Chemistry-Climate Model, CNRM-CCM: description and improvements from the CCMVal-2 simulations Geoscientific Model Development Discussions, 4, 1129-1183, 2011 Author(s): M. Michou, D. Saint-Martin, H. Teyssèdre, A. Alias, F. Karcher, D. Olivié, A. Voldoire, B. Josse, V.-H. Peuch, H. Clark, J. N. Lee, and F. Chéroux This paper presents a new version of the Météo-France CNRM Chemistry-Climate Model, so-called CNRM-CCM. It includes some fundamental changes from the previous version (CNRM-ACM) which was extensively evaluated in the context of the CCMVal-2 validation activity. The most notable changes concern the radiative code of the GCM, and the inclusion of the detailed stratospheric chemistry of our Chemistry-Transport model MOCAGE on-line within the GCM. A 47-yr transient simulation (1960–2006) is the basis of our analysis. CNRM-CCM generates satisfactory dynamical and chemical fields in the stratosphere. Several shortcomings of CNRM-ACM simulations for CCMVal-2 that resulted from an erroneous representation of the impact of volcanic aerosols as well as from transport deficiencies have been eliminated. Remaining problems concern the upper stratosphere (5 to 1 hPa) where temperatures are too high, and where there are biases in the NO 2 , N 2 O 5 and O 3 mixing ratios. In contrast, temperatures at the tropical tropopause are too cold. These issues are addressed through the implementation of a more accurate radiation scheme at short wavelengths. Despite these problems we show that this new CNRM CCM is a useful tool to study chemistry-climate applications.

Description: Pliocene Model Intercomparison Project: implementation strategy and mid-Pliocene Global climatology using GENESIS v3.0 GCM Geoscientific Model Development Discussions, 4, 2577-2603, 2011 Author(s): S. J. Koenig, R. M. DeConto, and D. Pollard The mid-Pliocene Warm Period (3.29 to 2.97 Ma BP) has been identified as an analogue for the future, with the potential to help understand climate processes in a warmer than modern world. Sets of climate proxies, combined to provide boundary conditions for Global Climate Model (GCM) simulations of the mid-Pliocene, form the basis for the international, data-driven Pliocene Model Intercomparison Project (PlioMIP). Here, we outline the strategy for implementing pre-industrial (modern) and mid-Pliocene forcings and boundary conditions into the GENESIS version 3 GCM, as part of PlioMIP. We describe the prescription of greenhouse gas concentrations and orbital parameters and the implementation of geographic boundary conditions such as land-ice-sea distribution, topography, sea surface temperatures, sea ice extent, vegetation, soils, and ice sheets. We further describe model-specific details including spin-up and integration times. In addition, the global climatology of the mid-Pliocene as simulated by the GENESIS v3 GCM is analyzed and compared to the pre-industrial control simulation. The simulated climate of the mid-Pliocene warm interval is found to differ considerably from pre-industrial. We identify model sensitivity to imposed forcings, and internal feedbacks that collectively affect both local and far-field responses. Our analysis points out the need to assess both the direct impacts of external forcings and the combined effects of indirect, internal feedbacks. This paper provides the basis for assessing model biases within the PlioMIP framework, and will be useful for comparisons with other studies of mid-Pliocene climates.

Description: Development of the high-order decoupled direct method in three dimensions for particulate matter: enabling advanced sensitivity analysis in air quality models Geoscientific Model Development Discussions, 4, 2605-2633, 2011 Author(s): W. Zhang, S. L. Capps, Y. Hu, A. Nenes, S. L. Napelenok, and A. G. Russell The high-order decoupled direct method in three dimensions for particular matter (HDDM-3D/PM) has been implemented in the Community Multiscale Air Quality (CMAQ) model to enable advanced sensitivity analysis. The major effort of this work is to develop high-order DDM sensitivity analysis of ISORROPIA, the inorganic aerosol module of CMAQ. A case-specific approach has been applied, and the sensitivities of activity coefficients and water content are explicitly computed. Stand-alone tests are performed for ISORROPIA by comparing the sensitivities (first- and second-order) computed by HDDM and the brute force (BF) approximations. Similar comparison has also been carried out for CMAQ results simulated using a week-long winter episode for a continental US domain. Second-order sensitivities of aerosol species (e.g., sulfate, nitrate, and ammonium) with respect to domain-wide SO 2 , NO x , and NH 3 emissions show agreement with BF results, yet exhibit less noise in locations where BF results are demonstrably inaccurate. Second-order sensitivity analysis elucidates nonlinear responses of secondary inorganic aerosols to their precursors and competing species that have not yet been well-understood with other approaches. Including second-order sensitivity coefficients in the Taylor series projection of the nitrate concentrations with a 50% reduction in domain-wide NO x emission shows a statistically significant improvement compared to the first-order Taylor series projection.

Description: Pliocene Ice Sheet Modelling Intercomparison Project (PLISMIP) – experimental design Geoscientific Model Development Discussions, 4, 2661-2686, 2011 Author(s): A. M. Dolan, S. J. Koenig, D. J. Hill, A. M. Haywood, and R. M. DeConto During the mid-Pliocene Warm Period (3.264 to 3.025 million yr ago), global mean temperature was similar to that predicted for the next century and atmospheric carbon dioxide concentrations were slightly higher. Sea level was also higher than today, implying a reduction in the extent of the ice sheets. Thus, the mid-Pliocene Warm Period provides a unique testing ground to investigate the stability of the Earth's ice sheets and their contribution to sea level in a warmer-than-modern world. Climate models and ice sheet models can be used to enhance our understanding of ice sheet stability, however, uncertainties associated with different ice-sheet modelling frameworks/approaches mean that a rigorous comparison of numerical ice sheet model simulations for the Pliocene is essential. As an extension to the Pliocene Model Intercomparison Project (PlioMIP; Haywood et al., 2010, 2011a), the Pliocene Ice Sheet Modelling Intercomparison Project (PLISMIP) will address these uncertainties. Here we outline the PLISMIP experimental design and initialisation conditions that have been adopted to simulate the Greenland and Antarctic ice sheets under present day and warm mid-Pliocene conditions. Not only will this project provide a new benchmark in the simulation of ice sheets in a past warm period, but the analysis of model sensitivity to various uncertainties could directly inform future predictions of ice sheet and sea level change.

Description: A web service based tool to plan atmospheric research flights Geoscientific Model Development Discussions, 4, 2123-2163, 2011 Author(s): M. Rautenhaus, G. Bauer, and A. Dörnbrack We present a web service based tool for the planning of atmospheric research flights. The tool provides online access to horizontal maps and vertical cross-sections of numerical weather prediction data and in particular allows the interactive design of a flight route in direct relation to the predictions. It thereby fills a crucial gap in the set of currently available tools for using data from numerical atmospheric models for research flight planning. A distinct feature of the tool is its lightweight, web service based architecture, requiring only commodity hardware and a basic Internet connection for deployment. Access to visualisations of prediction data is achieved by using an extended version of the Open Geospatial Consortium Web Map Service (WMS) standard, a technology that has gained increased attention in meteorology in recent years. With the WMS approach, we avoid the transfer of large forecast model output datasets while enabling on-demand generated visualisations of the predictions at campaign sites with limited Internet bandwidth. Usage of the Web Map Service standard also enables access to third-party sources of georeferenced data. We have implemented the software using the open-source programming language Python. In the present article, we describe the architecture of the tool. As an example application, we discuss a case study research flight planned for the scenario of the 2010 Eyjafjalla volcano eruption. Usage and implementation details are provided as Supplement.

Description: The CSIRO Mk3L climate system model version 1.0 – Part 2: Response to external forcings Geoscientific Model Development Discussions, 4, 3363-3435, 2011 Author(s): S. J. Phipps, L. D. Rotstayn, H. B. Gordon, J. L. Roberts, A. C. Hirst, and W. F. Budd The CSIRO Mk3L climate system model is a coupled general circulation model, designed primarily for millennial-scale climate simulation and palaeoclimate research. Mk3L includes components which describe the atmosphere, ocean, sea ice and land surface, and combines computational efficiency with a stable and realistic control climatology. It is freely available to the research community. This paper evaluates the response of the model to external forcings which correspond to past and future changes in the climate system. A simulation of the mid-Holocene climate is performed, in which changes in the seasonal and meridional distribution of incoming solar radiation are imposed. Mk3L correctly simulates increased summer temperatures at northern mid-latitudes and cooling in the tropics. However, it is unable to capture some of the regional-scale features of the mid-Holocene climate, with the precipitation over Northern Africa being deficient. The model simulates a reduction of between 7 and 15% in the amplitude of El Niño-Southern Oscillation, a smaller decrease than that implied by the palaeoclimate record. However, the realism of the simulated ENSO is limited by the model's relatively coarse spatial resolution. Transient simulations of the late Holocene climate are then performed. The evolving distribution of insolation is imposed, and an acceleration technique is applied and assessed. The model successfully captures the temperature changes in each hemisphere and the upward trend in ENSO variability. However, the lack of a dynamic vegetation scheme does not allow it to simulate an abrupt desertification of the Sahara. To assess the response of Mk3L to other forcings, transient simulations of the last millennium are performed. Changes in solar irradiance, atmospheric greenhouse gas concentrations and volcanic emissions are applied to the model. The model is again broadly successful at simulating larger-scale changes in the climate system. Both the magnitude and the spatial pattern of the simulated 20th century warming are consistent with observations. However, the model underestimates the magnitude of the relative warmth associated with the Mediaeval Climate Anomaly. Finally, three transient simulations are performed, in which the atmospheric CO 2 concentration is stabilised at two, three and four times the pre-industrial value. All three simulations exhibit ongoing surface warming, reduced sea ice cover, and a reduction in the rate of North Atlantic Deep Water formation followed by its gradual recovery. Antarctic Bottom Water formation ceases, with the shutdown being permanent for a trebling and quadrupling of the CO 2 concentration. The transient and equilibrium climate sensitivities of the model are determined. The short-term transient response to a doubling of the CO 2 concentration at 1% per year is a warming of 1.5 K, while the long-term equilibrium response is a warming of at least 3.9 K.

Description: The Nexus Land-Use model version 1.0, an approach articulating biophysical potentials and economic dynamics to model competition for land-use Geoscientific Model Development Discussions, 5, 571-638, 2012 Author(s): F. Souty, T. Brunelle, P. Dumas, B. Dorin, P. Ciais, R. Crassous, C. Müller, and A. Bondeau Interactions between food demand, biomass energy and forest preservation are driving both food prices and land-use changes, regionally and globally. This study presents a new model called Nexus Land-Use version 1.0 which describes these interactions through a generic representation of agricultural intensification mechanisms. The Nexus Land-Use model equations combine biophysics and economics into a single coherent framework to calculate crop yields, food prices, and resulting pasture and cropland areas within 12 regions inter-connected with each other by international trade. The representation of cropland and livestock production systems in each region relies on three components: (i) a biomass production function derived from the crop yield response function to inputs such as industrial fertilisers; (ii) a detailed representation of the livestock production system subdivided into an intensive and an extensive component, and (iii) a spatially explicit distribution of potential (maximal) crop yields prescribed from the Lund-Postdam-Jena global vegetation model for managed Land (LPJmL). The economic principles governing decisions about land-use and intensification are adapted from the Ricardian rent theory, assuming cost minimisation for farmers. The land-use modelling approach described in this paper entails several advantages. Firstly, it makes it possible to explore interactions among different types of biomass demand for food and animal feed, in a consistent approach, including indirect effects on land-use change resulting from international trade. Secondly, yield variations induced by the possible expansion of croplands on less suitable marginal lands are modelled by using regional land area distributions of potential yields, and a calculated boundary between intensive and extensive production. The model equations and parameter values are first described in details. Then, idealised scenarios exploring the impact of forest preservation policies or rising energy price on agricultural intensification are described, and their impacts on pasture and cropland areas are investigated.

Description: Modelling sub-grid wetland in the ORCHIDEE global land surface model: evaluation against river discharges and remotely sensed data Geoscientific Model Development Discussions, 5, 683-735, 2012 Author(s): B. Ringeval, B. Decharme, S. L. Piao, P. Ciais, F. Papa, N. de Noblet-Ducoudré, C. Prigent, P. Friedlingstein, I. Gouttevin, C. Koven, and A. Ducharne The quality of the global hydrological simulations performed by Land Surface Models (LSMs) strongly depends on processes that occur at unresolved spatial scales. Approaches such as TOPMODEL have been developed, which allow soil moisture redistribution within each grid-cell, based upon sub-grid scale topography. Moreover, the coupling between TOPMODEL and a LSM appears as a potential way to simulate wetland extent dynamic and its sensitivity to climate, a recently identified research problem for biogeochemical modelling, including methane emissions. Global evaluation of the coupling between TOPMODEL and an LSM is difficult, and prior attempts have been indirect, based on the evaluation of the simulated river flow. This study presents a new way to evaluate this coupling, within the ORCHIDEE LSM, using remote sensing data of inundated areas. Because of differences in nature between the satellite derived information – inundation extent – and the variable diagnosed by TOPMODEL/ORCHIDEE – area at maximum soil water content –, the evaluation focuses on the spatial distribution of these two quantities as well as on their temporal variation. Despite some difficulties in exactly matching observed localized inundated events, we obtain a rather good agreement in the distribution of these two quantities at a global scale. Floodplains are not accounted for in the model, and this is a major limitation. The difficulty of reproducing the year-to-year variability of the observed inundated area (for instance, the decreasing trend by the end of 90s) is also underlined. Classical indirect evaluation based on comparison between simulated and observed riverflow is also performed and underlines difficulties to simulate riverflow after coupling with TOPMODEL. The relationship between inundation and river flow at the basin scale in the model is analyzed, using both methods (evaluation against remote sensing data and riverflow). Finally, we discuss the potential of the TOPMODEL/LSM coupling to simulate wetland areas. A major limitation of the coupling for this purpose is linked to its ability to simulate a global wetland coverage consistent with the commonly used datasets. However, it seems to be a good opportunity to account for the wetland areas sensitivity to the climate and thus to simulate its temporal variability.

Description: Evaluation of a near-global eddy-resolving ocean model Geoscientific Model Development Discussions, 5, 4305-4354, 2012 Author(s): P. R. Oke, D. A. Griffin, A. Schiller, R. J. Matear, R. Fiedler, J. Mansbridge, A. Lenton, M. Cahill, M. A. Chamberlain, and K. Ridgway Analysis of the variability in an 18-yr run of a near-global, eddy-resolving ocean general circulation model coupled with biogeochemistry is presented. Comparisons between modelled and observed mean sea level (MSL), mixed-layer depth (MLD), sea-level anomaly (SLA), sea-surface temperature (SST), and Chlorophyll a indicate that the model variability is realistic. We find some systematic errors in the modelled MLD, with the model generally deeper than observations, that results in errors in the Chlorophyll a , owing to the strong biophysical coupling. We evaluate several other metrics in the model, including the zonally-averaged seasonal cycle of SST, meridional overturning, volume transports through key Straits and passages, zonal averaged temperature and salinity, and El Nino-related SST indices. We find that the modelled seasonal cycle in SST is 0.5–1.5 °C weaker than observed; volume transports of the Antarctic Circumpolar Current, the East Australian Current, and Indonesian Throughflow are in good agreement with observational estimates; and the correlation between the modelled and observed NINO SST indices exceed 0.91. Most aspects of the model circulation are realistic. We conclude that the model output is suitable for broader analysis to better understand ocean dynamics and ocean variability.

Description: The ICON-1.2 hydrostatic atmospheric dynamical core on triangular grids – Part 1: Formulation and performance of the baseline version Geoscientific Model Development Discussions, 6, 59-119, 2013 Author(s): H. Wan, M. A. Giorgetta, G. Zängl, M. Restelli, D. Majewski, L. Bonaventura, K. Fröhlich, D. Reinert, P. Rípodas, and L. Kornblueh A hydrostatic atmospheric dynamical core is developed for the purpose of global climate modelling. The model applies finite-difference methods to discretize the primitive equations on spherical icosahedral grids, using C-type staggering with triangles as control volumes for mass. This paper documents the numerical methods employed in the baseline version of the model, discusses their properties, and presents results from various idealized test cases. The evaluation shows that the new dynamical core is able to correctly represent the evolution of baroclinic eddies in the atmosphere as well as their role in heat and momentum transport. The simulations compare well with the reference solutions, and show a clear trend of convergence as the horizontal resolution increases. First results from two aqua-planet simulations are also presented, in which the equatorial wave spectra derived from tropical precipitation agree well with those simulated by a spectral transform model. The new dynamical core thus provides a good basis for further model development. Certain aspects of the model formulation that need further investigation and improvement are also pointed out.

Description: The GREENROOF module (v7.3) for modelling green roof hydrological and energetic performances within TEB Geoscientific Model Development Discussions, 6, 1127-1172, 2013 Author(s): C. S. de Munck, A. Lemonsu, R. Bouzouidja, V. Masson, and R. Claverie The need to prepare cities for climate change adaptation requests the urban modeller community to implement within their models sustainable adaptation strategies to be tested against specific city morphologies and scenarios. Greening city roofs is part of these strategies. In this context, a GREENROOF module for TEB (Town Energy Balance) has been developed to model the interactions between buildings and green roof systems at the scale of the city. This module allows one to describe an extensive green roof composed of four functional layers (vegetation – grasses or sedums, substrate, retention/drainage layers and artificial roof layers) and to model vegetation-atmosphere fluxes of heat, water and momentum, as well as the hydrological and thermal fluxes throughout the substrate and the drainage layers, and the thermal coupling with the structural building envelope. TEB-GREENROOF (v7.3) is therefore able to represent the impact of climate forcings on the functioning of the green roof vegetation and, conversely, the influence of the green roof on the local climate. A calibration exercise to adjust the model to the peculiar hydrological characteristics of the substrates and drainage layers commonly found on green roofs is performed for a case study located in Nancy (France) which consists of an extensive green roof with sedums. Model results for the optimum hydrological calibration show a good dynamics for the substrate water content which is nevertheless under-estimated but without impacting too much the green roof temperatures since they present a good agreement with observations. These results are encouraging with regard to modelling the impact of green roofs on thermal indoor comfort and energy consumption at the scale of cities, for which GREENROOF will be running with the building energy version of TEB, TEB-BEM. Moreover, the green roof studied for GREENROOF evaluation being a city-widespread type of extensive green roof, the hydrological characteristics derived through the evaluation exercise will be used as the standard configuration to model extensive green roofs at the scale of cities.

Description: The chemistry CATT–BRAMS model (CCATT–BRAMS 4.5): a regional atmospheric model system for integrated air quality and weather forecasting and research Geoscientific Model Development Discussions, 6, 1173-1222, 2013 Author(s): K. M. Longo, S. R. Freitas, M. Pirre, V. Marécal, L. F. Rodrigues, J. Panetta, M. F. Alonso, N. E. Rosário, D. S. Moreira, M. S. Gácita, J. Arteta, R. Fonseca, R. Stockler, D. M. Katsurayama, A. Fazenda, and M. Bela The Coupled Chemistry Aerosol-Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CCATT–BRAMS, version 4.5) is an online regional chemical transport model designed for local and regional studies of atmospheric chemistry from surface to the lower stratosphere suitable both for operational and research purposes. It includes gaseous/aqueous chemistry, photochemistry, scavenging and dry deposition. The CCATT–BRAMS model takes advantages of the BRAMS specific development for the tropics/subtropics and of the recent availability of preprocessing tools for chemical mechanisms and of fast codes for photolysis rates. BRAMS includes state-of-the-art physical parameterizations and dynamic formulations to simulate atmospheric circulations of scales down to meters. The online coupling between meteorology and chemistry allows the system to be used for simultaneous atmospheric weather and chemical composition forecasts as well as potential feedbacks between them. The entire system comprises three preprocessing software tools for chemical mechanism (which are user defined), aerosol and trace gases emission fields and atmospheric and chemistry fields for initial and boundary conditions. In this paper, the model description is provided along evaluations performed using observational data obtained from ground-based stations, instruments aboard of aircrafts and retrieval from space remote sensing. The evaluation takes into account model application on different scales from megacities and Amazon Basin up to intercontinental region of the Southern Hemisphere.

Description: Parallel algorithms for planar and spherical Delaunay construction with an application to centroidal Voronoi tessellations Geoscientific Model Development Discussions, 6, 1427-1466, 2013 Author(s): D. W. Jacobsen, M. Gunzburger, T. Ringler, J. Burkardt, and J. Peterson A new algorithm, featuring overlapping domain decompositions, for the parallel construction of Delaunay and Voronoi tessellations is developed. Overlapping allows for the seamless stitching of the partial Delaunay tessellations constructed by individual processors. The algorithm is then modified, by the addition of stereographic projections, to handle the parallel construction of spherical Delaunay and Voronoi tessellations. The algorithms are then embedded into algorithms for the parallel construction of planar and spherical centroidal Voronoi tessellations that require multiple constructions of Delaunay tessellations. Computational tests are used to demonstrate the efficiency and scalability of the algorithms for spherical Delaunay and centroidal Voronoi tessellations. Compared to serial versions of the algorithm and to the STRIPACK-based approaches, the new parallel algorithm results in significant speedups for the construction of spherical centroidal Voronoi tessellations.

Description: Evaluating a lightning parameterization based on cloud-top height for mesoscale numerical model simulations Geoscientific Model Development Discussions, 5, 3493-3531, 2012 Author(s): J. Wong, M. C. Barth, and D. Noone The Price and Rind lightning parameterization based on cloud-top height is a commonly used method for predicting flash rate in global chemistry models. As mesoscale simulations begin to implement flash rate predictions at resolutions that partially resolve convection, it is necessary to validate and understand the behavior of this method within such regime. In this study, we tested the flash rate parameterization, intra-cloud/cloud-to-ground (IC:CG) partitioning parameterization, and the associated resolution dependency "calibration factor" by Price and Rind using the Weather Research and Forecasting (WRF) model running at 36 km, 12 km, and 4 km grid spacings within the continental United States. Our results show that while the integrated flash count is consistent with observation when model biases in convection are taken into account, an erroneous frequency distribution is simulated. When the spectral characteristics of lightning flash rate is a concern, we recommend the use of prescribed IC:CG values. In addition, using cloud-top from convective parameterization, the "calibration factor" is also shown to be insufficient in reconciling the resolution dependency at the tested grid spacing used in this study. We recommend scaling by areal ratio relative to a base-case grid spacing determined by convective core density.

Description: Using model reduction to predict the soil-surface C 18 OO flux: an example of representing complex biogeochemical dynamics in a computationally efficient manner Geoscientific Model Development Discussions, 5, 3469-3491, 2012 Author(s): W. J. Riley Earth System Models (ESMs) must calculate large-scale interactions between the land and atmosphere while accurately characterizing fine-scale spatial heterogeneity in water, carbon, and nutrient dynamics. We present here a high-dimensional model representation (HDMR) approach that allows detailed process representation of a coupled carbon and water tracer (the δ 18 O value of the soil-surface CO 2 flux (δ F s )) in a computationally tractable manner. δ F s depends on the δ 18 O value of soil water, soil moisture, soil temperature, and soil CO 2 production (all of which are depth-dependent), and the δ 18 O value of above-surface CO 2 . We tested the HDMR approach over a growing season in a C 4 -dominated pasture using two vertical soil discretizations. The difference between the HDMR approach and the full model solution in the three-month integrated isoflux was less than 0.2% (0.5 mol m −2 ‰), and the approach is up to 100 times faster than the full numerical solution. This type of model reduction approach allows representation of complex coupled biogeochemical processes in regional and global climate models and can be extended to characterize subgrid-scale spatial heterogeneity.

Description: ECOCLIMAP-II/Europe: a twofold database of ecosystems and surface parameters at 1-km resolution based on satellite information for use in land surface, meteorological and climate models Geoscientific Model Development Discussions, 5, 3573-3620, 2012 Author(s): S. Faroux, A. T. Kaptué Tchuenté, J.-L. Roujean, V. Masson, E. Martin, and P. Le Moigne The overall objective of the present study is to introduce the new ECOCLIMAP-II database for Europe, which is an upgrade for this region of the former initiative, ECOCLIMAP-I, already implemented at global scale. The ECOCLIMAP programme is a dual database at 1-km resolution that includes an ecosystem classification and a coherent set of land surface parameters that are primarily mandatory in meteorological modelling (notably leaf area index and albedo). Hence, the aim of this innovative physiography is to enhance the quality of initialisation and impose some surface attributes within the scope of weather forecasting and climate related studies. The strategy for implementing ECOCLIMAP-II is to depart from prevalent land cover products such as CLC2000 (Corine Land Cover) and GLC2000 (Global Land Cover) by splitting existing classes into new classes that possess a better regional character by virtue of the climatic environment (latitude, proximity to the sea, topography). The leaf area index (LAI) from MODIS and NDVI from SPOT/Vegetation yield the two proxy variables that were considered here in order to perform a multi-year trimmed analysis between 1999 and 2005 using the K-means method. Further, meteorological applications require each land cover type to appear as a partition of fractions of 4 main surface types or tiles (nature, water bodies, sea, urban areas) and, inside the nature tile, fractions of 12 Plant Functional Types (PFTs) representing generic vegetation types – principally broadleaf forest, needleleaf forest, C3 and C4 crops, grassland and bare land – as incorporated by the SVAT model ISBA developed at Météo France. This landscape division also forms the cornerstone of a validation exercise. The new ECOCLIMAP-II can be verified with auxiliary land cover products at very fine and coarse resolutions by means of versatile land occupation nomenclatures.

Description: Intercomparison of temperature trends in IPCC CMIP5 simulations with observations, reanalyses and CMIP3 models Geoscientific Model Development Discussions, 5, 3621-3645, 2012 Author(s): J. Xu and A. M. Powell Jr. On the basis of the fifth Coupled Model Intercomparison Project (CMIP5) and the climate model simulations covering 1979 through 2005, the temperature trends and their uncertainties have been examined to note the similarities or differences compared to the radiosonde observations, reanalyses and the third Coupled Model Intercomparison Project (CMIP3) simulations. The results show noticeable discrepancies for the estimated temperature trends in the four data groups (Radiosonde, Reanalysis, CMIP3 and CMIP5) although similarities can be observed. Compared to the CMIP3 model simulations, the simulation in some of CMIP5 models were improved. The CMIP5 models displayed a negative temperature trend in the stratosphere closer to the strong negative trend seen in the observations. However, the positive tropospheric trend in the tropics is overestimated by the CMIP5 models relative to CMIP3 models. While some of the models produce temperature trend patterns more highly correlated with the observed patterns in CMIP5, the other models (such as CCSM4 and IPSL_CM5A-LR) exhibit the reverse tendency. The CMIP5 temperature trend uncertainty was significantly reduced in most areas, especially in the Arctic and Antarctic stratosphere, compared to the CMIP3 simulations. Similar to the CMIP3, the CMIP5 simulations overestimated the tropospheric warming in the tropics and Southern Hemisphere and underestimated the stratospheric cooling. The crossover point where tropospheric warming changes into stratospheric cooling occurred near 100 hPa in the tropics, which is higher than in the radiosonde and reanalysis data. The result is likely related to the overestimation of convective activity over the tropical areas in both the CMIP3 and CMIP5 models. Generally, for the temperature trend estimates associated with the numerical models including the reanalyses and global climate models, the uncertainty in the stratosphere is much larger than that in the troposphere, and the uncertainty in the Antarctic is the largest. In addition, note that the reanalyses show the largest uncertainty in the lower tropical stratosphere, and the CMIP3 simulations show the largest uncertainty in both the south and north polar regions.

Description: Present state of global wetland extent and wetland methane modelling: methodology of a model intercomparison project (WETCHIMP) Geoscientific Model Development Discussions, 5, 4071-4136, 2012 Author(s): R. Wania, J. R. Melton, E. L. Hodson, B. Poulter, B. Ringeval, R. Spahni, T. Bohn, C. A. Avis, G. Chen, A. V. Eliseev, P. O. Hopcroft, W. J. Riley, Z. M. Subin, H. Tian, V. Brovkin, P. M. van Bodegom, T. Kleinen, Z. C. Yu, J. S. Singarayer, S. Zürcher, D. P. Lettenmaier, D. J. Beerling, S. N. Denisov, C. Prigent, F. Papa, and J. O. Kaplan The Wetland and Wetland CH 4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH 4 ) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO 2 ) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2012). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO 2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extents and CH 4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH 4 fluxes from wetland systems. Even though modelling of wetland extents and CH 4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH 4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models.

Description: Improving the representation of secondary organic aerosol (SOA) in the MOZART-4 global chemical transport model Geoscientific Model Development Discussions, 5, 4187-4232, 2012 Author(s): A. Mahmud and K. C. Barsanti The secondary organic aerosol (SOA) module in the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) has been updated by replacing existing two-product (2p) parameters with those obtained from two-product volatility basis set (2p-VBS) fits, and by treating SOA formation from the following volatile organic compounds (VOCs): isoprene, propene and lumped alkenes. Strong seasonal and spatial variations in global SOA distributions were demonstrated, with significant differences in the predicted concentrations between the base-case and updated model versions. The base-case MOZART-4 predicted annual average SOA of 0.36 ± 0.50 μg m −3 in South America, 0.31 ± 0.38 μg m −3 in Indonesia, 0.09 ± 0.05 μg m −3 in the USA, and 0.12 ± 0.07 μg m −3 in Europe. Concentrations from the updated versions of the model showed a~marked increase in annual average SOA. Using the updated set of parameters alone (MZ4-v1) increased annual average SOA by ~8%, ~16%, ~56%, and ~108% from the base-case in South America, Indonesia, USA, and Europe, respectively. Treatment of additional parent VOCs (MZ4-v2) resulted in an even more dramatic increase of ~178–406% in annual average SOA for these regions over the base-case. The increases in predicted SOA concentrations further resulted in increases in corresponding SOA contributions to annual average total aerosol optical depth (AOD) by

Description: Modeling agriculture in the Community Land Model Geoscientific Model Development Discussions, 5, 4137-4185, 2012 Author(s): B. Drewniak, J. Song, J. Prell, V. R. Kotamarthi, and R. Jacob The potential impact of climate change on agriculture is uncertain. In addition, agriculture could influence above- and below-ground carbon storage. Development of models that represent agriculture is necessary to address these impacts. We have developed an approach to integrate agriculture representations for three crop types – maize, soybean, and spring wheat – into the coupled carbon-nitrogen version of the Community Land Model (CLM), to help address these questions. Here we present the new model, CLM-Crop, validated against observations from two AmeriFlux sites in the United States, planted with maize and soybean. Seasonal carbon fluxes compared well with field measurements. CLM-Crop yields were comparable with observations in some regions, although the generality of the crop model and its lack of technology and irrigation made direct comparison difficult. CLM-Crop was compared against the standard CLM3.5, which simulates crops as grass. The comparison showed improvement in gross primary productivity in regions where crops are the dominant vegetation cover. Crop yields and productivity were negatively correlated with temperature and positively correlated with precipitation. In case studies with the new crop model looking at impacts of residue management and planting date on crop yield, we found that increased residue returned to the litter pool increased crop yield, while reduced residue returns resulted in yield decreases. Using climate controls to signal planting date caused different responses in different crops. Maize and soybean had opposite reactions: when low temperature threshold resulted in early planting, maize responded with a loss of yield, but soybean yields increased. Our improvements in CLM demonstrate a new capability in the model – simulating agriculture in a realistic way, complete with fertilizer and residue management practices. Results are encouraging, with improved representation of human influences on the land surface and the potentially resulting climate impacts.

Description: A test of numerical instability and stiffness in the parametrizations of the ARPÉGE and ALADIN models Geoscientific Model Development Discussions, 5, 4233-4268, 2012 Author(s): M. Tudor Meteorological numerical weather prediction (NWP) models solve a system of partial differential equations in time and space. Semi-lagrangian advection scheme in the model dynamics allows for long time-steps. These longer time-steps can result in instabilities occurring in the model physics. A system of differential equations in which some solution components decay more rapidly than others is stiff. In this case it is stability rather than accuracy that restricts the time-step. The vertical diffusion parametrization can cause fast non-meteorological oscillations around the slowly evolving true solution (fibrillations). These are treated with an anti-fibrillation scheme. But small oscillations remain in an operational weather forecasts using ARPÉGE and ALADIN models. It is needed to test of the complete model formulation, as implemented in the operational forecast. In this paper, a simple test is designed. The test reveals if the formulation of particular physical parametrization is a stiff problem or potentially numerically unstable in combination with any other part of the model. When the test is applied to a stable scheme, the solution remains stable. But, applying the test to a potentially unstable scheme yields a solution with fibrillations of substantial amplitude. The parametrizations of a NWP model ARPÉGE were tested one by one to see which one may be the source of unstable model behaviour. The test has identified the stratiform precipitation scheme (a diagnostic Kessler type scheme) as a stiff problem, particularly the term that describes the evaporation of snow.

Description: PEATBOG: a biogeochemical model for analyzing coupled carbon and nitrogen dynamics in northern peatlands Geoscientific Model Development Discussions, 6, 1599-1688, 2013 Author(s): Y. Wu and C. Blodau Elevated nitrogen deposition and climate change alter the vegetation communities and carbon (C) and nitrogen (N) cycling in peatlands. To address this issue we developed a new process-oriented biogeochemical model (PEATBOG) for analyzing coupled carbon and nitrogen dynamics in northern peatlands. The model consists of four submodels, which simulate: (1) daily water table depth and depth profiles of soil moisture, temperature and oxygen levels; (2) competition among three plants functional types (PFTs), production and litter production of plants; (3) decomposition of peat; and (4) production, consumption, diffusion and export of dissolved C and N species in soil water. The model is novel in the integration of the C and N cycles, the explicit spatial resolution belowground, the consistent conceptualization of movement of water and solutes, the incorporation of stoichiometric controls on elemental fluxes and a consistent conceptualization of C and N reactivity in vegetation and soil organic matter. The model was evaluated for the Mer Bleue Bog, near Ottawa, Ontario, with regards to simulation of soil moisture and temperature and the most important processes in the C and N cycles. Model sensitivity was tested for nitrogen input, precipitation, and temperature, and the choices of the most uncertain parameters were justified. A simulation of nitrogen deposition over 40 yr demonstrates the advantages of the PEATBOG model in tracking biogeochemical effects and vegetation change in the ecosystem.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 3: A paleoperspective based on present-day data-model comparison for oxygen stable isotopes in carbonates Geoscientific Model Development Discussions, 6, 1527-1558, 2013 Author(s): T. Caley and D. M. Roche Oxygen stable isotopes ( 18 O) are among the most usual tools in paleoclimatology/paleoceanography. Simulation of oxygen stable isotopes allows testing how the past variability of these isotopes in water can be interpreted. By modelling the proxy directly in the model, the results can also be directly compared with the data. Water isotopes have been implemented in the global three-dimensional model of intermediate complexity i LOVECLIM allowing fully coupled atmosphere-ocean simulations. In this study, we present the validation of the model results for present day climate against global database for oxygen stable isotopes in carbonates. The limitation of the model together with the processes operating in the natural environment reveal the complexity of use the continental calcite 18 O signal of speleothems for a data-model comparison exercise. On the contrary, the reconstructed surface ocean calcite δ 18 O signal in i LOVECLIM does show a very good agreement with late Holocene database (foraminifers) at the global and regional scales. Our results indicate that temperature and the isotopic composition of the seawater are the main control on the fossil δ 18 O signal recorded in foraminifer shells and that depth habitat and seasonality play a role but have secondary importance. We argue that a data-model comparison for surface ocean calcite δ 18 O in past climate, such as the last glacial maximum (≈21 000 yr), could constitute an interesting tool for mapping the potential shifts of the frontal systems and circulation changes throughout time. Similarly, the potential changes in intermediate oceanic circulation systems in the past could be documented by a data (benthic foraminifers)-model comparison exercise whereas future investigations are necessary in order to quantitatively compare the results with data for the deep ocean.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 2: Evaluation of model results against observed δ 18 O in water samples Geoscientific Model Development Discussions, 6, 1495-1525, 2013 Author(s): D. M. Roche and T. Caley The H 2 18 O stable isotope was previously introduced in the three coupled components of the Earth System Model i LOVECLIM: atmosphere, ocean and vegetation. The results of a long (5000 yr) pre-industrial equilibrium simulation are presented and evaluated against measurement of H 2 18 O abundance in present-day water for the atmospheric and oceanic components. For the atmosphere, it is found that the model reproduces the observed spatial distribution and relationships to climate variables with some merit, though limitations following our approach are highlighted. Indeed, we obtain the main gradients with a robust representation of the Rayleigh distillation but caveats appear in Antarctica and around the Mediterranean region due to model limitation. For the oceanic component, the agreement between the modelled and observed distribution of water δ 18 O is found to be very good. Mean ocean surface latitudinal gradients are faithfully reproduced as well as the mark of the main intermediate and deep water masses. This opens large prospects for the applications in paleoclimatic context.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 1: Implementation and verification Geoscientific Model Development Discussions, 6, 1467-1494, 2013 Author(s): D. M. Roche A new 18 O stable water isotope scheme is developed for three components of the i LOVECLIM coupled climate model: atmospheric, oceanic and land surface. The equations required to reproduce the fractionation of stable water isotopes in the simplified atmospheric model ECBilt are developed consistently with the moisture scheme. Simplifications in the processes are made to account for the simplified vertical structure including only one moist layer. Implementation of these equations together with a passive tracer scheme for the ocean and a equilibrium fractionation scheme for the land surface leads to the closure of the (isotopic-)water budget in our climate system. Following the implementation, verification of the existence of usual δ 18 O to climatic relationships are performed for the Rayleigh distillation, the Dansgaard relationship and the δ 18 O–salinity relationship. Advantages and caveats of the approach taken are outlined.

Description: Capabilities and performance of Elmer/Ice, a new generation ice-sheet model Geoscientific Model Development Discussions, 6, 1689-1741, 2013 Author(s): O. Gagliardini, T. Zwinger, F. Gillet-Chaulet, G. Durand, L. Favier, B. de Fleurian, R. Greve, M. Malinen, C. Martín, P. Råback, J. Ruokolainen, M. Sacchettini, M. Schäfer, H. Seddik, and J. Thies The Fourth IPCC Assessment Report concluded that ice-sheet flow models are unable to forecast the current increase of polar ice sheet discharge and the associated contribution to sea-level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice-flow models, and as a result, a significant number of new ice-sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger scale problems, earning the status of an ice-sheet model. Here, we summarise almost 10 yr of development performed by different groups. We present the components already included in Elmer/Ice, its numerical performance, selected applications, as well as developments planned for the future.

Description: WRFv3.2-SPAv2: development and validation of a coupled ecosystem-atmosphere model, scaling from surface fluxes of CO 2 and energy to atmospheric profiles Geoscientific Model Development Discussions, 6, 1559-1598, 2013 Author(s): T. L. Smallman, J. B. Moncrieff, and M. Williams The Weather Research & Forecasting meteorological (WRF) model has been coupled to the Soil Plant Atmosphere (SPA) terrestrial ecosystem model, to produce WRF-SPA. SPA generates realistic land-atmosphere exchanges through fully coupled hydrological, carbon and energy cycles. The addition of a land surface model (SPA) capable of modelling biospheric CO 2 exchange allows WRF-SPA to be used for investigating the feedbacks between biosphere carbon balance, meteorology and land management/land use change. We have extensively validated WRF-SPA using multi-annual observations of air temperature, turbulent fluxes, net radiation and net ecosystem exchange of CO 2 at three sites, representing the dominant vegetation types in Scotland (forest, managed grassland and arable agriculture). WRF-SPA generates more realistic seasonal behaviour at the site level compared to an unmodified version of WRF, and produces realistic CO 2 exchanges. WRF-SPA is also able to realistically model atmospheric profiles of CO 2 over Scotland, spanning a 3 yr period (2004–2006), capturing both profile structure, indicating realistic transport, and magnitude indicating appropriate source sink distribution and CO 2 exchange. WRF-SPA makes use of CO 2 tracer pools and can therefore identify and quantify land surface contributions to the modelled atmospheric CO 2 signal at a specified location.

Description: MESMO 2: a mechanistic marine silica cycle and coupling to a simple terrestrial scheme Geoscientific Model Development Discussions, 5, 2999-3033, 2012 Author(s): K. Matsumoto, K. S. Tokos, A. Huston, and H. Joy-Warren Here we describe the second version of Minnesota Earth System Model for Ocean biogeochemistry (MESMO 2), an earth system model of intermediate complexity, which consists of a dynamical ocean, dynamic-thermodynamic sea ice, and energy moisture balanced atmosphere. The new version has more realistic land ice masks and is driven by seasonal winds. A major aim in version 2 is representing the marine silica cycle mechanistically in order to investigate climate-carbon feedbacks involving diatoms, a critically important class of phytoplankton in terms of carbon export production. This is achieved in part by including iron, on which phytoplankton uptake of silicic acid depends. Also, MESMO 2 is coupled to an existing terrestrial model, which allows for the exchange of carbon, water, and energy between land and the atmosphere. The coupled model, called MESMO 2E, is appropriate for more complete earth system simulations. The new version was calibrated with the goal of preserving reasonable interior ocean ventilation and various biological production rates in the ocean and land, while simulating key features of the marine silica cycle.

Description: A new method to diagnose the contribution of anthropogenic activities to temperature: temperature tagging Geoscientific Model Development Discussions, 5, 3183-3215, 2012 Author(s): V. Grewe This study presents a new methodology, called temperature tagging. It keeps track of the contributions of individual processes to temperature within a climate model simulation. As a first step and as a test bed a simple box climate model is regarded. The model consists of an atmosphere, which absorbs and emits radiation and of a surface, which reflects, absorbs and emits radiation. The tagging methodology is used to investigate the impact of the atmosphere on surface temperature. Four processes are investigated in more detail and their contribution to the surface temperature quantified: (i) shortwave influx and shortwave atmospheric absorption ("sw"), (ii) longwave atmospheric absorption due to non-CO 2 greenhouse gases ("nC"), (iii) due to a base case CO 2 concentration ("bC"), and (iv) due to an enhanced CO 2 concentration ("eC"). The differential equation for the temperature in the box climate model is decomposed into four equations for the tagged temperatures. This method is applied to investigate the contribution of longwave absorption to the surface temperature (greenhouse effect), which is calculated to be 68 K. This estimate contrasts an alternative calculation of the greenhouse effect of slightly more than 30 K based on the difference of the surface temperature with and without an atmosphere. The difference of the two estimates is due to a shortwave cooling effect and a reduced contribution of the shortwave to the total downward flux: The shortwave absorption of the atmosphere results in a reduced net shortwave flux at the surface of 192 W m −2 , leading to a cooling of the surface by 14 K. Introducing an atmosphere results in a downward longwave flux at the surface due to atmospheric absorption of 189 W m −2 , which roughly equals the net shortwave flux of 192 W m −2 . This longwave flux is a result of both, the radiation due to atmospheric temperatures and its longwave absorption. Hence the longwave absorption roughly accounts for 91 W m −2 out of a total of 381 W m −2 (roughly 25%) and therefore accounts for a temperature of 68 K. In a second experiment, the CO 2 concentration is doubled, which leads to an increase in surface temperature of 1.2 K, resulting from an temperature increase due to CO 2 of 1.9 K, due to non-CO 2 greenhouse gases of 0.6 K and a cooling of 1.3 K due to a reduced importance of the solar heating for the surface and atmospheric temperatures. These two experiments show the feasibility of temperature tagging and its potential as a diagnostic for climate simulations.

Description: Implementation of the Fast-JX Photolysis scheme into the UKCA component of the MetUM chemistry climate model Geoscientific Model Development Discussions, 5, 3217-3260, 2012 Author(s): P. J. Telford, N. L. Abraham, A. T. Archibald, P. Braesicke, M. Dalvi, O. Morgenstern, F. M. O'Connor, N. A. D. Richards, and J. A. Pyle Atmospheric chemistry is driven by photolytic reactions, making their modelling a crucial component of atmospheric models. We describe the implementation and validation of Fast-JX, a state of the art model of interactive photolysis, into the MetUM chemistry climate model. This allows for interactive photolysis frequencies to be calculated in the troposphere and augments the calculation of the frequencies in the stratosphere by accounting for clouds and aerosols in addition to ozone. In order to demonstrate the effectiveness of this new photolysis scheme we employ new methods of validating the model, including techniques for sampling the model to compare to flight track and satellite data.

Description: COSTRICE – three model online coupling using OASIS: problems and solutions Geoscientific Model Development Discussions, 5, 3261-3310, 2012 Author(s): H. T. M. Ho, B. Rockel, H. Kapitza, B. Geyer, and E. Meyer The coupled system COSTRICE is developed for the first time in order to reproduce the interactions and feedbacks between atmosphere, ocean and sea-ice in a two-way online coupled model system containing three component models for regional climate simulations over Baltic Sea and North Sea regions. The regional climate model CCLM 1 is coupled to the regional ocean model TRIMNP 1 and the sea ice model CICE 1 via the coupler OASIS3. In this study, CCLM is setup with a horizontal grid mesh size of 50 km and 32 vertical atmosphere layers and driven by the 6-h ERA-interim reanalysis data as initial and boundary conditions. TRIMNP is setup with a horizontal grid mesh size of 12.8 km and 50 vertical ocean levels. CICE calculates ice in 5 categories and runs with the same horizontal resolution as TRIMNP but only over the Baltic Sea and the Kattegat Bay of the North Sea. In a two-way online coupling process, CCLM is linked to TRIMNP through sea surface temperature (SST) as lower boundary condition every 3 h and TRIMNP is driven by 1-h atmospheric state variables and fluxes of CCLM. The data exchange processes between TRIMNP and CICE as well as from CCLM to CICE take place with an interval of 3 h. The coupled model is applied in a study for climate simulations over Baltic Sea and North Sea regions in 1997. The coupled system is set up to run in parallel on the super computing system IBM-power 6 at the German Climate Computing Center (DKRZ). 1 See Table A1.

Description: A generalized tagging method Geoscientific Model Development Discussions, 5, 3311-3324, 2012 Author(s): V. Grewe The understanding of causes of changes in climate-chemistry simulations is an important, but often challenging task. In atmospheric chemistry, one approach is to tag species according to their origin (e.g. emission categories) and to inherit these tags to other species during subsequent reactions. This concept was recently employed to calculate the contribution of atmospheric processes to temperature. Here a new concept for tagging any state variable is presented. This generalized tagging method results from a sensitivity analysis of the forcing terms of the right hand side of the governing differential equations. In a couple of examples, the consistency with previous approaches is shown. Since the method is based on a ratio describing relative sensitivities, singularities occur where the method is not applicable. For some applications, like in atmospheric chemistry, these singularities can easily be removed. However, one theoretical example is given, where this method is not applicable at all.

Description: Technical Note: Improving computational efficiency in large linear inverse problems: an example from carbon dioxide flux estimation Geoscientific Model Development Discussions, 5, 3325-3342, 2012 Author(s): V. Yadav and A. M. Michalak Addressing a variety of questions within Earth science disciplines entails the inference of the spatio-temporal distribution of parameters of interest based on observations of related quantities. Such estimation problems often represent inverse problems that are formulated as linear optimization problems. Computational limitations arise when the number of observations and/or the size of the discretized state space become large, especially if the inverse problem is formulated in a probabilistic framework and therefore aims to assess the uncertainty associated with the estimates. This work proposes two approaches to lower the computational costs and memory requirements for large linear space-time inverse problems, taking the Bayesian approach for estimating carbon dioxide (CO 2 ) emissions and uptake (a.k.a. fluxes) as a prototypical example. The first algorithm can be used to efficiently multiply two matrices, as long as one can be expressed as a Kronecker product of two smaller matrices, a condition that is typical when multiplying a sensitivity matrix by a covariance matrix in the solution of inverse problems. The second algorithm can be used to compute a posteriori uncertainties directly at aggregated spatio-temporal scales, which are the scales of most interest in many inverse problems. Both algorithms have significantly lower memory requirements and computational complexity relative to direct computation of the same quantities (O( n 2.5 ) vs. O( n 3 )). For an examined benchmark problem, the two algorithms yielded a three and six order of magnitude increase in computational efficiency, respectively, relative to direct computation of the same quantities. Sample computer code is provided for assessing the computational and memory efficiency of the proposed algorithms for matrices of different dimensions.

Description: Quality assessment concept of the World Data Center for Climate and its application to CMIP5 data Geoscientific Model Development Discussions, 5, 781-802, 2012 Author(s): M. Stockhause, H. Höck, F. Toussaint, and M. Lautenschlager The preservation of data in a high state of quality and suitable for interdisciplinary use is one of the most pressing and challenging current issues in long-term archiving. For high volume data such as climate model data, the data and data replica are no longer stored centrally but distributed over several local data repositories, e.g. the data of the Climate Model Intercomparison Project No. 5 (CMIP5). The most important part of the data is to be published as DOI according to the World Data Center for Climate's (WDCC) application of the DataCite regulations. The integrated part of WDCC's data publication process, the data quality assessment, was adapted to the requirements of a federated data infrastructure. A concept of a distributed and federated quality assessment procedure was developed, in which the work load and responsibility for quality control is shared between the three primary CMIP5 data centers: Program for Climate Model Diagnosis and Intercomparison (PCMDI), British Atmospheric Data Centre (BADC), and WDCC. This distributed quality control concept, its pilot implementation for CMIP5, and first experiences are presented.

Description: Sensitivities and uncertainties of modeled ground temperatures in mountain environments Geoscientific Model Development Discussions, 6, 791-840, 2013 Author(s): S. Gubler, S. Endrizzi, S. Gruber, and R. S. Purves Before operational use or for decision making, models must be validated, and the degree of trust in model outputs should be quantified. Often, model validation is performed at single locations due to the lack of spatially-distributed data. Since the analysis of parametric model uncertainties can be performed independently of observations, it is a suitable method to test the influence of environmental variability on model evaluation. In this study, the sensitivities and uncertainty of a physically-based mountain permafrost model are quantified within an artificial topography consisting of different elevations and exposures combined with six ground types characterized by their hydraulic properties. The analyses performed for all combinations of topographic factors and ground types allowed to quantify the variability of model sensitivity and uncertainty within mountain regions. We found that modeled snow duration considerably influences the mean annual ground temperature (MAGT). The melt-out day of snow (MD) is determined by processes determining snow accumulation and melting. Parameters such as the temperature and precipitation lapse rate and the snow correction factor have therefore a great impact on modeled MAGT. Ground albedo changes MAGT from 0.5 to 4°C in dependence of the elevation, the aspect and the ground type. South-exposed inclined locations are more sensitive to changes in ground albedo than north-exposed slopes since they receive more solar radiation. The sensitivity to ground albedo increases with decreasing elevation due to shorter snow cover. Snow albedo and other parameters determining the amount of reflected solar radiation are important, changing MAGT at different depths by more than 1°C. Parameters influencing the turbulent fluxes as the roughness length or the dew temperature are more sensitive at low elevation sites due to higher air temperatures and decreased solar radiation. Modeling the individual terms of the energy balance correctly is hence crucial in any physically-based permafrost model, and a separate evaluation of the energy fluxes could substantially improve the results of permafrost models. The sensitivity in the hydraulic properties change considerably for different ground types: rock or clay for instance are not sensitive while gravel or peat, accurate measurements of the hydraulic properties could significantly improve modeled ground temperatures. Further, the discretization of ground, snow and time have an impact on modeled MAGT that cannot be neglected (more than 1°C for several discretization parameters). We show that the temporal resolution should be at least one hour to ensure errors less than 0.2°C in modeled MAGT, and the uppermost ground layer should at most be 20 mm thick. Within the topographic setting, the total parametric output uncertainties expressed as the standard deviation of the Monte Carlo model simulations range from 0.1 to 0.5°C for clay, silt and rock, and from 0.1 to 0.8°C for peat, sand and gravel. These uncertainties are comparable to the variability of ground surface temperatures measured within 10 m × 10 m grids in Switzerland. The increased uncertainties for sand, peat and gravel is largely due to the high hydraulic conductivity.

Description: One-dimensional simulation of fire injection heights in contrasted meteorological scenarios with PRM and Meso-NH models Geoscientific Model Development Discussions, 6, 721-790, 2013 Author(s): S. Strada, S. R. Freitas, C. Mari, K. M. Longo, and R. Paugam Wild-fires release huge amounts of aerosol and hazardous trace gases in the atmosphere. The residence time and the dispersion of fire pollutants in the atmosphere can range from hours to days and from local to continental scales. These various scenarios highly depend on the injection height of smoke plumes. The altitude at which fire products are injected in the atmosphere is controlled by fire characteristics and meteorological conditions. Injection height however is still poorly accounted in chemistry transport models for which fires are sub-grid scale processes which need to be parametrised. Only recently, physically-based approaches for estimating the fire injection heights have been developed which consider both the convective updrafts induced by the release of fire sensible heat and the impact of background meteorological environment on the fire convection dynamics. In this work, two different models are used to simulate fire injection heights in contrasted meteorological scenarios: a Mediterranean arson fire and two Amazonian deforestation fires. A Eddy-Diffusivity/Mass-Flux approach, formerly developed to reproduce convective boundary layer in the non-hydrostatic meteorological model Meso-NH, is compared to the 1-D Plume Rise Model. For both models, radiosonde data and re-analyses from the European Center for Medium-Range Weather Forecasts (ECMWF) have been used as initial conditions to explore the sensitivity of the models responses to different meteorological forcings. The two models predict injection heights for the Mediterranean fire between 1.7 and 3.3 km with the Meso-NH/EDMF model systematically higher than the 1-D PRM model. Both models show a limited sensitivity to the meteorological forcings with a 20–30% difference in the injection height between radiosondes and ECMWF data for this case. Injection heights calculated for the two Amazonian fires ranges from 5 to 6.5 km for the 1-D PRM model and from 2 to 4 km for the Meso-NH/EDMF model. The difference of smoke plume heights between the two models can reach 3–4 km. A large difference is obtained for the windy-wet Amazonian fire by the 1-D PRM model with a injection height 1.5 km higher when ECMWF re-analyses are used compared to the run with the radiosonde forcing. For the Mediterranean case, both models forecast a plume injection height above the boundary layer, although there are evidences that this particular fire propagated near the surface, highlighting the current limitations of the two approaches.

Description: Representation of nucleation mode microphysics in global aerosol microphysics models Geoscientific Model Development Discussions, 6, 893-924, 2013 Author(s): Y. H. Lee, J. R. Pierce, and P. J. Adams In models, nucleation mode (1 nm 〈 D p 〈 10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively. We also investigate the impact of the time step for aerosol microphysical processes (a 10-min versus a 1-h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e. 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e. 1 nm boundary) with the 10-min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles ( J 10 ) and the burdens and lifetimes of ultrafine mode (10 nm 〈 D p 〈 70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e. 10 nm and 3 nm boundaries) result in higher J 10 and shorter coagulation lifetimes of ultrafine mode particles than the model with explicit dynamics (i.e. 1 nm boundary). The spatial distributions of CN10 ( D p 〉 10 nm) and CCN(0.2%) (i.e. CCN concentrations at 0.2% supersaturation) are moderately affected, especially CN10 predictions above ~ 700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J 10 or J 3 with the parameterized microphysics possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2%) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10–40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e. 1-h time step used in the 10 nm simulations compared to 10-min time step used in the benchmark simulations) but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

Description: An optimally tuned ensemble of the "eb_go_gs" configuration of GENIE: parameter sensitivity and bifurcations in the Atlantic overturning circulation Geoscientific Model Development Discussions, 6, 925-956, 2013 Author(s): R. Marsh, A. Sóbester, E. E. Hart, K. I. C. Oliver, N. R. Edwards, and S. J. Cox The key physical parameters of the "eb_go_gs" configuration of GENIE, an Earth system Model of Intermediate Complexity (EMIC), are tuned using a multi-objective genetic algorithm. An ensemble of 90 parameter sets is tuned using two ocean and two atmospheric state variables as targets. These are "Pareto-optimal", representing a range of trade-offs between the four tuning targets. For the leading five parameter sets, simulations are further evaluated alongside a simulation with untuned "default" parameters, comparing selected variables and diagnostics that describe the state of the atmosphere, ocean and sea ice. One of these parameter sets is selected for further analysis of the objective function (error) landscape in the vicinity of its tuned values. "Cliffs" along some dimensions motivate closer inspection of corresponding variations in the Atlantic meridional overturning circulation (AMOC). This reveals that bifurcations in the AMOC are highly sensitive to parameters that are not commonly associated with MOC stability. Specifically, the state of the AMOC is sensitive to parameters governing the wind-driven circulation and atmospheric heat transport. Five optimal parameter sets are recommended for future use of GENIE in the configuration presented here.

Description: An efficient method to generate a perturbed parameter ensemble of a fully coupled AOGCM without flux-adjustment Geoscientific Model Development Discussions, 6, 841-892, 2013 Author(s): P. J. Irvine, L Gregoire, D. J. Lunt, and P. J. Valdes We present a simple method to generate a perturbed parameter ensemble (PPE) of a fully-coupled atmosphere-ocean general circulation model (AOGCM), HadCM3, without requiring flux-adjustment. The aim was to produce an ensemble that samples parametric uncertainty in some key variables and displays a similar range of behavior as seen in multi-model ensembles (MMEs). Six atmospheric parameters, a sea-ice parameter and an ocean parameter were jointly perturbed within a reasonable range to generate an initial group of 200 members. To screen out implausible ensemble members, 20 yr pre-industrial control simulations were run and members whose temperature response to the parameter perturbations was projected to be outside the range of 13.6 ± 2°C, i.e. near to the observed pre-industrial global mean, were discarded. 21 members, including the standard unperturbed model, were accepted, covering almost the entire span of the eight parameters, challenging the argument that without flux-adjustment parameter ranges would be unduly restricted. This ensemble was used in 3 experiments; a 800 yr pre-industrial, a 150 yr quadrupled CO 2 , and a 150 yr 1% CO 2 rise per annum simulation. The behavior of the PPE for the pre-industrial control compared well to the CMIP3 ensemble for a number of surface and atmospheric column variables with the exception of a few members in the Tropics. However, we find that members of the PPE with low values of the entrainment rate coefficient show very large increases in upper tropospheric and stratospheric water vapor concentrations in response to elevated CO 2 and some show implausibly high climate sensitivities, and as such some of these members will be excluded from future experiments with this ensemble. The outcome of this study is a PPE of a fully-coupled AOGCM which samples parametric uncertainty with a range of behavior similar to the CMIP3 ensemble and a simple methodology which would be applicable to other GCMs.

Description: A cloud chemistry module for the 3-D cloud-resolving mesoscale model Meso-NH with application to idealized cases Geoscientific Model Development Discussions, 6, 957-1020, 2013 Author(s): M. Leriche, J.-P. Pinty, C. Mari, and D. Gazen A complete chemical module has been developed for use in the Meso-NH three-dimensional cloud resolving mesoscale model. This module includes gaseous and aqueous phase chemical reactions that are analysed by a pre-processor generating the Fortran90 code automatically. The kinetic solver is based on a Rosenbrock algorithm, which is robust and accurate for integrating stiff systems and especially multiphase chemistry. The exchange of chemical species between the gas phase and cloud droplets and raindrops is computed kinetically by mass transfers considering non-equilibrium between the gas and the condensed phases. Microphysical transfers of chemical species are considered for the various cloud microphysics schemes available, which are based on one-moment or two-moment schemes. The pH of the droplets and of the raindrops is diagnosed separately as the root of a high order polynomial equation. The chemical concentrations in the ice phase are modelled in a single phase encompassing the two categories of precipitating ice particles (snow and graupel) of the microphysical scheme. The only process transferring chemical species in ice is retention during freezing or riming of liquid hydrometeors. Three idealized simulations are reported, which highlight the sensitivity of scavenging efficiency to the choice of the microphysical scheme and the retention coefficient in the ice phase. A two-dimensional warm, shallow convection case is used to compare the impact of the microphysical schemes on the temporal evolution and rates of acid precipitation. Acid wet deposition rates are shown to be overestimated when a one-moment microphysics scheme is used compared to a two-moment scheme. The difference is induced by a better prediction of raindrop radius and raindrop number concentration in the latter scheme. A two-dimensional mixed-phase squall line and a three-dimensional mixed-phase supercell were simulated to test the sensitivity of cloud vertical transport to the retention efficiency of gases in the ice phase. The 2-D and 3-D simulations illustrate that the retention in ice of a moderately soluble gas such as formaldehyde substantially decreases its concentration in the upper troposphere. In these simulations, retention of highly soluble species in the ice phase significantly increased the wet deposition rates.

Description: GAPPARD: a computationally efficient method of approximating gap-scale disturbance in vegetation models Geoscientific Model Development Discussions, 6, 1021-1084, 2013 Author(s): M. Scherstjanoi, J. O. Kaplan, E. Thürig, and H. Lischke Models of vegetation dynamics that are designed for application at spatial scales larger than individual forest gaps suffer from several limitations. Typically, either a population average approximation is used that results in unrealistic tree allometry and forest stand structure, or models have a high computational demand because they need to simulate both a series of age-based cohorts and a number of replicate patches to account for stochastic gap-scale disturbances. The detail required by the latter method increases the number of calculations by two to three orders of magnitude compared to the less realistic population average approach. In an effort to increase the efficiency of dynamic vegetation models without sacrificing realism, and to explore patterns of spatial scaling in forests, we developed a new method for simulating stand-replacing disturbances that is both accurate and 10-50x faster than approaches that use replicate patches. The GAPPARD (approximating GAP model results with a Probabilistic Approach to account for stand Replacing Disturbances) method works by postprocessing the output of deterministic, undisturbed simulations of a cohort-based vegetation model by deriving the distribution of patch ages at any point in time on the basis of a disturbance probability. With this distribution, the expected value of any output variable can be calculated from the output values of the deterministic undisturbed run at the time corresponding to the patch age. To account for temporal changes in model forcing, e.g., as a result of climate change, GAPPARD performs a series of deterministic simulations and interpolates between the results in the postprocessing step. We integrated the GAPPARD method in the forest models LPJ-GUESS and TreeM-LPJ, and evaluated these in a series of simulations along an altitudinal transect of an inner-alpine valley. With GAPPARD applied to LPJ-GUESS results were insignificantly different from the output of the original model LPJ-GUESS using 100 replicate patches, but simulation time was reduced by approximately the factor 10. Our new method is therefore highly suited rapidly approximating LPJ-GUESS results, and provides the opportunity for future studies over large spatial domains, allows easier parameterization of tree species, faster identification of areas of interesting simulation results, and comparisons with large-scale datasets and forest models.

Description: A refined statistical cloud closure using double-Gaussian probability density functions Geoscientific Model Development Discussions, 6, 1085-1125, 2013 Author(s): A. K. Naumann, A. Seifert, and J. P. Mellado We introduce a probability density function (PDF) based scheme to parameterize cloud fraction, average liquid water and liquid water flux in large-scale models, that is developed from and tested against large-eddy simulations and observational data. Because the tails of the PDFs are crucial for an appropriate parameterization of cloud properties, we use a double-Gaussian distribution that is able to represent the observed, skewed PDFs properly. Introducing two closure equations, the resulting parameterization relies on the first three moments of the subgrid variability of temperature and moisture as input parameters. The parameterization is shown to be clearly superior to a single-Gaussian approach in diagnosing the cloud fraction and average liquid water profiles and improves existing double-Gaussian closures. We find that the error of the new parameterization is smallest for a horizontal resolution of about 5–20 km and also depends on the appearance of mesoscale structures that are accompanied by higher rain rates. In combination with simple autoconversion schemes that only depend on the liquid water, the error introduced by the new parameterization is orders of magnitude smaller than the difference between various autoconversion schemes. For the liquid water flux, we introduce a parameterization that is depending on the skewness of the subgrid variability of temperature and moisture and that reproduces the profiles of the liquid water flux well.

Description: Evaluation of WRF-SFIRE performance with field observations from the FireFlux experiment Geoscientific Model Development Discussions, 6, 121-169, 2013 Author(s): A. K. Kochanski, M. A. Jenkins, J. Mandel, J. D. Beezley, C. B. Clements, and S. Krueger This study uses in-situ measurements collected during the FireFlux field experiment to evaluate and improve the performance of coupled atmosphere-fire model WRF-SFIRE. The simulation of the experimental burn shows that WRF-SFIRE is capable of providing realistic head fire rate-of-spread and the vertical temperature structure of the fire plume, and, up to 10 m above ground level, fire-induced surface flow and vertical velocities within the plume. The model captured the changes in wind speed and direction before, during, and after fire front passage, along with arrival times of wind speed, temperature, and updraft maximae, at the two instrumented flux towers used in FireFlux. The model overestimated vertical velocities and underestimated horizontal wind speeds measured at tower heights above the 10 m, and it is hypothesized that the limited model resolution over estimated the fire front depth, leading to too high a heat release and, subsequently, too strong an updraft. However, on the whole, WRF-SFIRE fire plume behavior is consistent with the interpretation of FireFlux observations. The study suggests optimal experimental pre-planning, design, and execution of future field campaigns that are needed for further coupled atmosphere-fire model development and evaluation.

Description: SPITFIRE-2: an improved fire module for Dynamic Global Vegetation Models Geoscientific Model Development Discussions, 5, 2347-2443, 2012 Author(s): M. Pfeiffer and J. O. Kaplan Fire is the primary disturbance factor in many terrestrial ecosystems. Wildfire alters vegetation structure and composition, affects carbon storage and biogeochemical cycling, and results in the release of climatically relevant trace gases, including CO 2 , CO, CH 4 , NO x , and aerosols. Assessing the impacts of global wildfire on centennial to multi-millennial timescales requires the linkage of process-based fire modeling with vegetation modeling using Dynamic Global Vegetation Models (DGVMs). Here we present a new fire module, SPITFIRE-2, and an update to the LPJ-DGVM that includes major improvements to the way in which fire occurrence, behavior, and the effect of fire on vegetation is simulated. The new fire module includes explicit calculation of natural ignitions, the representation of multi-day burning and coalescence of fires and the calculation of rates of spread in different vegetation types, as well as a simple scheme to model crown fires. We describe a new representation of anthropogenic biomass burning under preindustrial conditions that distinguishes the way in which the relationship between humans and fire are different between hunter-gatherers, obligate pastoralists, and farmers. Where and when available, we evaluate our model simulations against remote-sensing based estimates of burned area. While wildfire in much of the modern world is largely influenced by anthropogenic suppression and ignitions, in those parts of the world where natural fire is still the dominant process, e.g. in remote areas of the boreal forest, our results demonstrate a significant improvement in simulated burned area over previous models. With its unique properties of being able to simulate preindustrial fire, the new module we present here is particularly well suited for the investigation of climate-human-fire relationships on multi-millennial timescales.

Description: The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): overview and description of models, simulations and climate diagnostics Geoscientific Model Development Discussions, 5, 2445-2502, 2012 Author(s): J.-F. Lamarque, D. T. Shindell, B. Josse, P. J. Young, I. Cionni, V. Eyring, D. Bergmann, P. Cameron-Smith, W. J. Collins, R. Doherty, S. Dalsoren, G. Faluvegi, G. Folberth, S. J. Ghan, L. W. Horowitz, Y. H. Lee, I. A. MacKenzie, T. Nagashima, V. Naik, D. Plummer, M. Righi, S. Rumbold, M. Schulz, R. B. Skeie, D. S. Stevenson, S. Strode, K. Sudo, S. Szopa, A. Voulgarakis, and G. Zeng The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) consists of a series of timeslice experiments targeting the long-term changes in atmospheric composition between 1850 and 2100, with the goal of documenting radiative forcing and the associated composition changes. Here we introduce the various simulations performed under ACCMIP and the associated model output. The ACCMIP models have a wide range of horizontal and vertical resolutions, vertical extent, chemistry schemes and interaction with radiation and clouds. While anthropogenic and biomass burning emissions were specified for all time slices in the ACCMIP protocol, it is found that the natural emissions lead to a significant range in emissions, mostly for ozone precursors. The analysis of selected present-day climate diagnostics (precipitation, temperature, specific humidity and zonal wind) reveals biases consistent with state-of-the-art climate models. The model-to-model comparison of changes in temperature, specific humidity and zonal wind between 1850 and 2000 and between 2000 and 2100 indicates mostly consistent results, but with outliers different enough to possibly affect their representation of climate impact on chemistry.

Description: Setup of the PMIP3 paleoclimate experiments conducted using an Earth System Model, MIROC-ESM Geoscientific Model Development Discussions, 5, 2527-2569, 2012 Author(s): T. Sueyoshi, R. Ohgaito, A. Yamamoto, M. O. Chikamoto, T. Hajima, H. Okajima, M. Yoshimori, M. Abe, R. O'ishi, F. Saito, S. Watanabe, M. Kawamiya, and A. Abe-Ouchi The importance of climate model evaluation using paleoclimate simulations for better future climate projections has been recognized by the Intergovernmental Panel on Climate Change. In recent years, Earth System Models (ESMs) were developed to investigate carbon-cycle climate feedback, as well as to project the future climate. Paleoclimate events, especially those associated with the variations in atmospheric CO 2 level or land vegetation, provide suitable benchmarks to evaluate ESMs. Here we present implementations of the paleoclimate experiments proposed by the Coupled Model Intercomparison Project phase 5/Paleoclimate Modelling Intercomparison Project phase 3 (CMIP5/PMIP3) using an Earth System Model, MIROC-ESM. In this paper, experimental settings and procedures of the mid-Holocene, the Last Glacial Maximum, and the Last Millennium experiments are explained. The first two experiments are time slice experiments and the last one is a transient experiment. The complexity of the model requires various steps to correctly configure the experiments. Several basic outputs are also shown.

Description: Aerosol-climate interactions in the Norwegian Earth System Model – NorESM Geoscientific Model Development Discussions, 5, 2599-2685, 2012 Author(s): A. Kirkevåg, T. Iversen, Ø. Seland, C. Hoose, J. E. Kristjánsson, H. Struthers, A. M. L. Ekman, S. Ghan, J. Griesfeller, E. D. Nilsson, and M. Schulz The objective of this study is to document and evaluate recent changes and updates to the module for aerosols and aerosol-cloud-radiation interactions in the atmospheric module CAM4-Oslo of the Norwegian Earth System Model (NorESM). Particular attention is paid to the role of natural organics, sea salt, and mineral dust in determining the gross aerosol properties as well as the anthropogenic contribution to these properties and the associated direct and indirect radiative forcing. The aerosol module is extended from earlier versions that have been published, and includes life-cycling of sea-salt, mineral dust, particulate sulphate, black carbon, and primary and secondary organics. The impacts of most of the numerous changes since previous versions are thoroughly explored by sensitivity experiments. The most important changes are: modified prognostic sea salt emissions; updated treatment of precipitation scavenging and gravitational settling; inclusion of biogenic primary organics and methane sulphonic acid (MSA) from oceans; almost doubled production of land-based biogenic secondary organic aerosols (SOA); and increased ratio of organic matter to organic carbon (OM / OC) for biomass burning aerosols from 1.4 to 2.6. Compared with in-situ measurements and remotely sensed data, the new treatments of sea salt and dust aerosols give smaller biases in near surface mass concentrations and aerosol optical depth than in the earlier model version. The model biases for mass concentrations are approximately unchanged for sulphate and BC. The enhanced levels of modeled OM yield improved overall statistics, even though OM is still underestimated in Europe and over-estimated in North America. The global direct radiative forcing (DRF) at the top of the atmosphere has changed from a small positive value to −0.08 W m −2 in CAM4-Oslo. The sensitivity tests suggest that this change can be attributed to the new treatment of biomass burning aerosols and gravitational settling. Although it has not been a goal in this study, the new DRF estimate is closer both to the median model estimate from the AeroCom inter-comparison and the best estimate in IPCC AR4. Estimated DRF at the ground surface has increased by ca. 60%, to −1.89 W m −2 . We show that this can be explained by new emission data and omitted mixing of constituents between updrafts and downdrafts in convective clouds. The increased abundance of natural OM and the introduction of a cloud droplet spectral dispersion formulation are the most important contributions to a considerably decreased estimate of the indirect radiative forcing (IndRF). The IndRF is also found to be sensitive to assumptions about the coating of insoluble aerosols by sulphate and OM. The IndRF of −1.2 W m −2 , which is closer to the IPCC AR4 estimates than the previous estimate of −1.9 W m −2 , has thus been obtained without imposing unrealistic artificial lower bounds on cloud droplet number concentrations.

Description: Inclusion of Ash and SO 2 emissions from volcanic eruptions in WRF-CHEM: development and some applications Geoscientific Model Development Discussions, 5, 2571-2597, 2012 Author(s): M. Stuefer, S. R. Freitas, G. Grell, P. Webley, S. Peckham, and S. A. McKeen We describe a new functionality within the Weather Research and Forecasting model with coupled Chemistry (WRF-Chem) that allows simulating emission, transport, dispersion, transformation and sedimentation of pollutants released during volcanic activities. Emissions from both an explosive eruption case and relatively calm degassing situation are considered using the most recent volcanic emission databases. A preprocessor tool provides emission fields and additional information needed to establish the initial three-dimensional cloud umbrella/vertical distribution within the transport model grid, as well as the timing and duration of an eruption. From this source condition, the transport, dispersion and sedimentation of the ash-cloud can be realistically simulated by WRF-Chem using its own dynamics, physical parameterization as well as data assimilation. Examples of model validation include a comparison of tephra fall deposits from the 1989 eruption of Mount Redoubt (Alaska), and the dispersion of ash from the 2010 Eyjafjallajökull eruption in Iceland. Both model applications show good coincidence between WRF-Chem and observations.

Description: A web-based software tool to estimate unregulated daily streamflow at ungauged rivers Geoscientific Model Development Discussions, 5, 2503-2526, 2012 Author(s): S. A. Archfield, P. A. Steeves, J. D. Guthrie, and K. G. Ries III Streamflow information is critical for solving any number of hydrologic problems. Often times, streamflow information is needed at locations which are ungauged and, therefore, have no observations on which to base water management decisions. Furthermore, there has been increasing need for daily streamflow time series to manage rivers for both human and ecological functions. To facilitate negotiation between human and ecological demands for water, this paper presents the first publically-available, map-based, regional software tool to interactively estimate daily streamflow time series at any user-selected ungauged river location. The map interface allows users to locate and click on a river location, which then returns estimates of daily streamflow for the location selected. For the demonstration region in the northeast United States, daily streamflow was shown to be reliably estimated by the software tool, with efficiency values computed from observed and estimated streamflows ranging from 0.69 to 0.92. The software tool provides a general framework that can be applied to other regions for which daily streamflow estimates are needed.

Description: Using multi-model averaging to improve the reliability of catchment scale nitrogen predictions Geoscientific Model Development Discussions, 5, 2289-2310, 2012 Author(s): J.-F. Exbrayat, N. R. Viney, H.-G. Frede, and L. Breuer Hydro-biogeochemical models are used to foresee the impact of mitigation measures on water quality. Usually, scenario-based studies rely on single model applications. This is done in spite of the widely acknowledged advantage of ensemble approaches to cope with structural model uncertainty issues. As an attempt to demonstrate the reliability of such multi-model efforts in the hydro-biogeochemical context, this methodological contribution proposes an adaptation of the Reliability Ensemble Averaging (REA) philosophy to nitrogen losses predictions. A total of 4 models are used to predict the total nitrogen (TN) losses from the well-monitored Ellen Brook catchment in Western Australia. Simulations include re-predictions of current conditions and a set of straightforward management changes targeting fertilization scenarios. Results show that, in spite of good calibration metrics, one of the models provides a very different response to management changes. This behaviour leads the simple average of the ensemble members to also predict reductions in TN export that are not in agreement with the other models. However, considering the convergence of model predictions in the more sophisticated REA approach assigns more weight to previously less well calibrated models that are more in agreement with each other. This method also avoids having to disqualify any of the ensemble members, which is always sensible.

Description: Development of high resolution land surface parameters for the Community Land Model Geoscientific Model Development Discussions, 5, 1435-1481, 2012 Author(s): Y. Ke, L. R. Leung, M. Huang, A. M. Coleman, H. Li, and M. S. Wigmosta There is a growing need for high-resolution land surface parameters as land surface models are being applied at increasingly higher spatial resolution offline as well as in regional and global models. The default land surface parameters for the most recent version of the Community Land Model (i.e. CLM 4.0) are at 0.5° or coarser resolutions, released with the model from the National Center for Atmospheric Research (NCAR). Plant Functional Types (PFTs), vegetation properties such as Leaf Area Index (LAI), Stem Area Index (SAI), and non-vegetated land covers were developed using remotely-sensed datasets retrieved in late 1990's and the beginning of this century. In this study, we developed new land surface parameters for CLM 4.0, specifically PFTs, LAI, SAI and non-vegetated land cover composition, at 0.05° resolution globally based on the most recent MODIS land cover and improved MODIS LAI products. Compared to the current CLM 4.0 parameters, the new parameters produced a decreased coverage by bare soil and trees, but an increased coverage by shrub, grass, and cropland. The new parameters result in a decrease in global seasonal LAI, with the biggest decrease in boreal forests; however, the new parameters also show a large increase in LAI in tropical forest. Differences between the new and the current parameters are mainly caused by changes in the sources of remotely sensed data and the representation of land cover in the source data. The new high-resolution land surface parameters have been used in a coupled land-atmosphere model (WRF-CLM) applied to the western US to demonstrate their use in high-resolution modeling. Future work will include global offline CLMsimulations to examine the impacts of source data resolution and subsequent land parameter changes on simulated land surface processes.

Description: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions Geoscientific Model Development Discussions, 5, 1503-1560, 2012 Author(s): A. B. Guenther, X. Jiang, C. L. Heald, T. Sakulyanontvittaya, T. Duhl, L. K. Emmons, and X. Wang The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1) is a modeling framework for estimating fluxes of 147 biogenic compounds between terrestrial ecosystems and the atmosphere using simple mechanistic algorithms to account for the major known processes controlling biogenic emissions. It is available as an offline code and has also been coupled into land surface models and atmospheric chemistry models. MEGAN2.1 is an update from the previous versions including MEGAN2.0 for isoprene emissions and MEGAN2.04, which estimates emissions of 138 compounds. Isoprene comprises about half of the estimated total global biogenic volatile organic compound (BVOC) emission of 1 Pg (1000 Tg or 10 15 g). Another 10 compounds including methanol, ethanol, acetaldehyde, acetone, α-pinene, β-pinene, t −β-ocimene, limonene, ethene, and propene together contribute another 30% of the estimated emission. An additional 20 compounds (mostly terpenoids) are associated with another 17% of the total emission with the remaining 3% distributed among 125 compounds. Emissions of 41 monoterpenes and 32 sesquiterpenes together comprise about 15% and 3%, respectively, of the total global BVOC emission. Tropical trees cover about 18% of the global land surface and are estimated to be responsible for 60% of terpenoid emissions and 48% of other VOC emissions. Other trees cover about the same area but are estimated to contribute only about 10% of total emissions. The magnitude of the emissions estimated with MEGAN2.1 are within the range of estimates reported using other approaches and much of the differences between reported values can be attributed to landcover and meteorological driving variables. The offline version of MEGAN2.1 source code and driving variables is available from http://acd.ucar.edu/~guenther/MEGAN/MEGAN.htm and the version integrated into the Community Land Model version 4 (CLM4) can be downloaded from http://www.cesm.ucar.edu/ .

Description: A simulation study of the ensemble-based data assimilation of satellite-borne lidar aerosol observations Geoscientific Model Development Discussions, 5, 1877-1947, 2012 Author(s): T. T. Sekiyama, T. Y. Tanaka, and T. Miyoshi A four-dimensional ensemble-based data assimilation system was assessed by observing system simulation experiments (OSSEs), in which the CALIPSO satellite was emulated via simulated satellite-borne lidar aerosol observations. Its performance over athree-month period was validated according to the Method for Object-based Diagnostic Evaluation (MODE), using aerosol optical thickness (AOT) distributions in East Asia as the objects of analysis. Consequently, this data assimilation system demonstrated the ability to produce better analyses of sulfate and dust aerosols in comparison to a free-running simulation model. For example, the mean centroid distance (from the truth) over a three-month collection period of aerosol plumes was improved from 2.15 grids (≈ 600 km) to 1.45 grids (≈ 400 km) for sulfate aerosols and from 2.59 grids (≈ 750 km) to 1.14 grids (≈ 330 km) for dust aerosols; the mean area ratio (to the truth) over a three-month collection period of aerosol plumes was improved from 0.49 to 0.76 for sulfate aerosols and from 0.51 to 0.72 for dust aerosols. The satellite-borne lidar data assimilation successfully improved the aerosol plume analysis and the dust emission estimation in the OSSEs. These results present great possibilities for the beneficial use of lidar data, whose distribution is vertically/temporally dense but horizontally sparse, when coupled with a four-dimensional data assimilation system. In addition, sensitivity tests were conducted, and their results indicated that the degree of freedom to control the aerosol variables was probably limited in the data assimilation because the meteorological field in the system was constrained to weather reanalysis using Newtonian relaxation. Further improvements to the aerosol analysis can be performed through the simultaneous assimilation of aerosol observations with meteorological observations. The OSSE results strongly suggest that the use of real CALIPSO data will have a beneficial effect on obtaining more accurate sulfate and dust aerosol analyses. Furthermore, the use of the same OSSE technique will allow us to perform a prior assessment of the next-generation lidar satellite EarthCARE, which will be launched in 2015.

Description: Sensitivity analysis and calibration of a soil carbon model (SoilGen2) in two contrasting loess forest soils Geoscientific Model Development Discussions, 5, 1817-1849, 2012 Author(s): Y. Y. Yu, P. A. Finke, H. B. Wu, and Z. T. Guo To accurately estimate past terrestrial carbon pools is the key to understand the global carbon cycle and its relationship with the climate system. SoilGen2 is a useful tool to obtain aspects of soil properties (including carbon content) by simulating soil formation processes; thus it offers an opportunity for past soil carbon pool reconstruction. In order to apply it to various environmental conditions, parameters related to carbon cycle process in SoilGen2 are calibrated based on 6 soil pedons from two typical loess deposition regions (Belgium and China). Sensitivity analysis using Morris' method shows that decomposition rate of humus ( k HUM ), fraction of incoming plant material as leaf litter (fr ecto ) and decomposition rate of resistant plant material ( k RPM ) are 3 most sensitive parameters that would cause the greatest uncertainty in simulated change of soil organic carbon in both regions. According to the principle of minimizing the difference between simulated and measured organic carbon by comparing quality indices, the suited values of k HUM , fr ecto and k RPM in the model are deduced step by step. The difference of calibrated parameters between Belgium and China may be attributed to their different vegetation types and climate conditions. This calibrated model is improved for better simulation of carbon change in the whole pedon and has potential for future modeling of carbon cycle in paleosols.

Description: Activation of the operational ecohydrodynamic model (3-D CEMBS) – the hydrodynamic part Geoscientific Model Development Discussions, 5, 1851-1875, 2012 Author(s): L. Dzierzbicka-Głowacka, J. Jakacki, M. Janecki, and A. Nowicki The paper presents a description of the hydrodynamic part of the coupled ice-ocean model that also includes ecosystem predictive model for evaluation of the condition of the marine environment and the Baltic ecosystem, as well as a preliminary empirical verification of the operational hydrodynamic model based on the POP code in order to determine the consistence between the results obtained from the model and experimental results for the sea surface temperature. The current Baltic Sea model is based on the Community Earth System Model (CESM from NCAR – National Center for Atmospheric Research). CESM was adopted for the Baltic Sea as a coupled sea-ice model. It consists of the Community Ice Code (CICE model, version 4.0) and the Parallel Ocean Program (POP, version 2.1). The models are coupled through the coupler (CPL7), which is based on the Model Coupling Toolkit (MCT) routines. The current horizontal resolution is about 2 km (1/48 degrees). The ocean model has 21 vertical levels. The driver time step is 1440 s and it is also coupling the time step. The ocean model time step is about 480 s (8 min). Currently, the model is forced by fields from the European Center for Medium Weather Forecast. In the operational mode, 48-h atmospheric forecasts are used, which are supplied by the UM model of the Interdisciplinary Centre for Mathematical and Computational Modelling of the Warsaw University. The model of the marine ecosystem is the right tool for monitoring the state and bioproductivity of the marine ecosystem and forecasting the physical and ecological changes in the studied basin.

Description: Tagged ozone mechanism for MOZART-4, CAM-chem, and other chemical transport models Geoscientific Model Development Discussions, 5, 1949-1985, 2012 Author(s): L. K. Emmons, P. G. Hess, J.-F. Lamarque, and G. G. Pfister A procedure for tagging ozone produced from NO sources through updates to an existing chemical mechanism is described, and results from its implementation in the Model for Ozone and Related chemical Tracers (MOZART-4), a global chemical transport model, are presented. Artificial tracers are added to the mechanism, thus not affecting the standard chemistry. The results are linear in the troposphere, i.e., the sum of ozone from individual tagged sources equals the ozone from all sources to within 3% in zonal mean monthly averages. The stratospheric ozone contribution to the troposphere determined from the difference between total ozone and ozone from all tagged sources is significantly less than estimates using a traditional stratospheric ozone tracer (8 vs 20 ppbv at the surface). The commonly used technique of perturbing NO emissions by 20% in a region to determine its ozone contribution is compared to the tagging technique, showing that the tagged ozone is 2–4 times the ozone contribution that was deduced from perturbing emissions.

Description: Coupling technologies for Earth System Modelling Geoscientific Model Development Discussions, 5, 1987-2006, 2012 Author(s): S. Valcke, V. Balaji, A. Craig, C. DeLuca, R. Dunlap, R. W. Ford, R. Jacob, J. Larson, R. O'Kuinghttons, G. D. Riley, and M. Vertenstein This paper presents a review of the software currently used in climate modelling in general and in CMIP5 in particular to couple the numerical codes representing the different components of the Earth system. The coupling technologies presented show common features, such as the ability to communicate and regrid data, but also offer different functions and implementations. Design characteristics of the different approaches are discussed as well as future challenges arising from the increasing complexity of scientific problems and computing platforms.

Description: Downscale cascades in tracer transport test cases: an intercomparison of the dynamical cores in the Community Atmosphere Model CAM5 Geoscientific Model Development Discussions, 5, 1781-1816, 2012 Author(s): J. Kent, C. Jablonowski, J. P. Whitehead, and R. B. Rood The accurate modelling of cascades to unresolved scales is an important part of the tracer transport component of dynamical cores of weather and climate models. This paper aims to investigate the ability of the advection schemes in the National Center for Atmospheric Research's Community Atmosphere Model version 5 (CAM5) to model this cascade. In order to quantify the effects of the different advection schemes in CAM5, four two-dimensional tracer transport test cases are presented. Three of the tests stretch the tracer below the scale of coarse resolution grids to ensure the downscale cascade of tracer variance. These results are compared with a high resolution reference solution, which is simulated on a resolution fine enough to resolve the tracer during the test. The fourth test has two separate flow cells, and is designed so that any tracer in the Western Hemisphere should not pass into the Eastern Hemisphere. This is to test whether the diffusion in transport schemes, often in the form of explicit hyper-diffusion terms or implicit through monotonic limiters, contains unphysical mixing. An intercomparison of three of the dynamical cores of the National Center for Atmospheric Research's Community Atmosphere Model version 5 is performed. The results show that the finite-volume (CAM-FV) and spectral element (CAM-SE) dynamical cores model the downscale cascade of tracer variance better than the semi-Lagrangian transport scheme of the Eulerian spectral transform core (CAM-EUL). Each scheme tested produces unphysical mass in the Eastern Hemisphere of the separate cells test.

Description: CLM4-BeTR, a generic biogeochemical transport and reaction module for CLM4: model development, evaluation, and application Geoscientific Model Development Discussions, 5, 2705-2744, 2012 Author(s): J. Tang, W. J. Riley, C. D. Koven, and Z. M. Subin To improve regional and global biogeochemistry modeling and climate predictability, we have developed a generic reactive transport module for the land model CLM4 (called CLM4-BeTR (Biogeochemical Transport and Reactions)). CLM4-BeTR represents the transport, interactions, and biotic and abiotic transformations of an arbitrary number of tracers (aka chemical species) in an arbitrary number of phases (e.g. dissolved, gaseous, sorbed, aggregate). An operator splitting approach was employed and consistent boundary conditions were derived for each modeled sub-process. Tracer fluxes, associated with hydrological processes such as surface run-on and run-off, belowground drainage, and ice to liquid conversion were also computed consistently with the bulk water fluxes calculated by the soil physics module in CLM4. The transport code was evaluated and found be in good agreement with several analytical test cases. The model was then applied at the Harvard Forest site with a representation of depth-dependent belowground biogeochemistry. The results indicated that, at this site, (1) CLM4-BeTR was able to simulate soil-surface CO 2 effluxes and soil CO 2 profiles accurately; (2) the transient surface CO 2 effluxes calculated based on the tracer transport mechanism were in general not equal to the belowground CO 2 production rates and that their differences varied according to the seasonal cycle of soil physics and biogeochemistry; (3) losses of CO 2 through processes other than surface gas efflux were less than 1% of the overall soil respiration; and (4) the contributions of root respiration and heterotrophic respiration have distinct temporal signals in surface CO 2 effluxes and soil CO 2 concentrations. The development of CLM4-BeTR will allow detailed comparisons between ecosystem observations and predictions and insights to the modeling of terrestrial biogeochemistry.

Description: The Norwegian Earth System Model, NorESM1-M – Part 2: Climate response and scenario projections Geoscientific Model Development Discussions, 5, 2933-2998, 2012 Author(s): T. Iversen, M. Bentsen, I. Bethke, J. B. Debernard, A. Kirkevåg, Ø. Seland, H. Drange, J. E. Kristjánsson, I. Medhaug, M. Sand, and I. A. Seierstad The NorESM1-M simulation results for CMIP5 ( http://cmip-pcmdi.llnl.gov/cmip5/index.html ) are described and discussed. Together with the accompanying paper by Bentsen et al. (2012), this paper documents that NorESM1-M is a valuable global climate model for research and for providing complementary results to the evaluation of possible man made climate change. NorESM is based on the model CCSM4 operated at NCAR on behalf of many contributors in USA. The ocean model is replaced by a developed version of MICOM and the atmospheric model is extended with on-line calculations of aerosols, their direct effect, and their indirect effect on warm clouds. Model validation is presented in a companion paper (Bentsen et al., 2012). NorESM1-M is estimated to have equilibrium climate sensitivity slightly smaller than 2.9 K, a transient climate response just below 1.4 K, and is less sensitive than most other models. Cloud feedbacks damp the response, and a strong AMOC reduces the heat fraction available for increasing near surface temperatures, for evaporation, and for melting ice. The future projections based on RCP scenarios yield global surface air temperature increase almost one standard deviation lower than a 15-model average. Summer sea-ice is projected to decrease considerably by 2100, and completely for RCP8.5. The AMOC is projected to reduce by 12%, 15–17%, and 32% for the RCP2.6, 4.5, 6.0 and 8.5 respectively. Precipitation is projected to increase in the tropics, decrease in the subtropics and in southern parts of the northern extra-tropics during summer, and otherwise increase in most of the extra-tropics. Changes in the atmospheric water cycle indicate that precipitation events over continents will become more intense and dry spells more frequent. Extra-tropical storminess in the Northern Hemisphere is projected to shift northwards. There are indications of more frequent spring and summer blocking in the Euro-Atlantic sectors and that ENSO events weaken but appear more frequent. These indications are uncertain because of biases in the model's representation of present-day conditions. There are indications that positive phase PNA and negative phase NAO become less frequent under the RCP8.5 scenario, but also this result is considered uncertain. Single-forcing experiments indicate that aerosols and greenhouse gases produce similar geographical patterns of response for near surface temperature and precipitation. These patterns tend to have opposite sign, with important exceptions for precipitation at low latitudes. The asymmetric aerosol effects between the two hemispheres leads to a southward displacement of ITCZ. Both forcing agents thus tend to reduce northern hemispheric subtropical precipitation.

Description: The Norwegian Earth System Model, NorESM1-M – Part 1: Description and basic evaluation Geoscientific Model Development Discussions, 5, 2843-2931, 2012 Author(s): M. Bentsen, I. Bethke, J. B. Debernard, T. Iversen, A. Kirkevåg, Ø. Seland, H. Drange, C. Roelandt, I. A. Seierstad, C. Hoose, and J. E. Kristjánsson The core version of the Norwegian Climate Center's Earth System Model, named NorESM1-M, is presented. The NorESM-family of models are based on the Community Climate System Model version 4 (CCSM4) of the University Corporation for Atmospheric Research, but differs from the latter by, in particular, an isopycnic coordinate ocean model and advanced chemistry-aerosol-cloud-radiation interaction schemes. NorESM1-M has a horizontal resolution of approximately 2° for the atmosphere and land components and 1° for the ocean and ice components. NorESM is also available in a lower resolution version (NorESM1-L) and a version that includes prognostic biogeochemical cycling (NorESM1-ME). The latter two model configurations are not part of this paper. Here, a first-order assessment of the model stability, the mean model state and the internal variability based on the model experiments made available to CMIP5 are presented. Further analysis of the model performance is provided in an accompanying paper (Iversen et al., 2012), presenting the corresponding climate response and scenario projections made with NorESM1-M.

Description: Simulations of the Mid-Pliocene Warm Period using the NASA/GISS ModelE2-R Earth System Model Geoscientific Model Development Discussions, 5, 2811-2842, 2012 Author(s): M. A. Chandler, L. E. Sohl, J. A. Jonas, and H. J. Dowsett Climate reconstructions of the mid-Pliocene Warm Period (mPWP) bear many similarities to aspects of future global warming as projected by the Intergovernmental Panel on Climate Change. In particular, marine and terrestrial paleoclimate data point to high latitude temperature amplification, with associated decreases in sea ice and land ice and altered vegetation distributions that show expansion of warmer climate biomes into higher latitudes. NASA GISS climate models have been used to study the Pliocene climate since the USGS PRISM project first identified that the mid-Pliocene North Atlantic sea surface temperatures were anomalously warm. Here we present the most recent simulations of the Pliocene using the AR5/CMIP5 version of the GISS Earth System Model known as ModelE2-R. These simulations constitute the NASA contribution to the Pliocene Model Intercomparison Project (PlioMIP) Experiment 2. Many findings presented here corroborate results from other PlioMIP multi-model ensemble papers, but we also emphasize features in the ModelE2-R simulations that are unlike the ensemble means. We provide discussion of features that show considerable improvement compared with simulations from previous versions of the NASA GISS models, improvement defined here as simulation results that more closely resemble the ocean core data as well as the PRISM3D reconstructions of the mid-Pliocene climate. In some regions even qualitative agreement between model results and paleodata are an improvement over past studies, but the dramatic warming in the North Atlantic and Greenland-Iceland-Norwegian Sea in these new simulations is by far the most accurate portrayal ever of this key geographic region by the GISS climate model. Our belief is that continued development of key physical routines in the atmospheric model, along with higher resolution and recent corrections to mixing parameterizations in the ocean model, have led to an Earth System Model that will produce more accurate projections of future climate.

Description: PORT, a CESM tool for the diagnosis of radiative forcing Geoscientific Model Development Discussions, 5, 2687-2704, 2012 Author(s): A. J. Conley, J.-F. Lamarque, F. Vitt, W. D. Collins, and J. Kiehl The Parallel Offline Radiative Transfer (PORT) model is a tool for diagnosing radiative forcing. It isolates the radiation code from the Community Atmosphere Model (CAM4) in the Community Earth System Model (CESM1). The computation of radiative forcing from doubling of carbon dioxide and from the change of ozone concentration from year 1850 to 2000 illustrates the use of PORT.

Description: A methodology for estimating seasonal cycles of atmospheric CO 2 resulting from terrestrial net ecosystem exchange (NEE) fluxes using the Transcom T3L2 pulse-response functions Geoscientific Model Development Discussions, 5, 2789-2809, 2012 Author(s): C. D. Nevison, D. F. Baker, and K. R. Gurney We present a method for translating modeled terrestrial net ecosystem exchange (NEE) fluxes of carbon into the corresponding seasonal cycles in atmospheric CO 2 . The method is based on the pulse-response functions from the Transcom 3 Level 2 (T3L2) atmospheric tracer transport model (ATM) intercomparison. The new pulse-response method is considerably faster than a full forward ATM simulation, allowing CO 2 seasonal cycles to be computed in seconds, rather than the days or weeks required for a forward simulation. Further, the results provide an estimate of the range of transport uncertainty across 13 different ATMs associated with the translation of surface NEE fluxes into an atmospheric signal. We evaluate the method against the results of archived forward ATM simulations from T3L2. The latter are also used to estimate the uncertainties associated with oceanic and fossil fuel influences. We present a regional breakdown at selected monitoring sites of the contribution to the atmospheric CO 2 cycle from the 11 different T3L2 land regions. A test case of the pulse-response code, forced by NEE fluxes from the Community Land Model, suggests that for many terrestrial models, discrepancies between model results and observed atmospheric CO 2 cycles will be large enough to clearly transcend ATM uncertainties.