ALBERT

Description: Pore formation during dehydration of a polycrystalline gypsum sample observed and quantified in a time-series synchrotron X-ray micro-tomography experiment Solid Earth, 3, 71-86, 2012 Author(s): F. Fusseis, C. Schrank, J. Liu, A. Karrech, S. Llana-Fúnez, X. Xiao, and K. Regenauer-Lieb We conducted an in-situ X-ray micro-computed tomography heating experiment at the Advanced Photon Source (USA) to dehydrate an unconfined 2.3 mm diameter cylinder of Volterra Gypsum. We used a purpose-built X-ray transparent furnace to heat the sample to 388 K for a total of 310 min to acquire a three-dimensional time-series tomography dataset comprising nine time steps. The voxel size of 2.2 μm 3 proved sufficient to pinpoint reaction initiation and the organization of drainage architecture in space and time. We observed that dehydration commences across a narrow front, which propagates from the margins to the centre of the sample in more than four hours. The advance of this front can be fitted with a square-root function, implying that the initiation of the reaction in the sample can be described as a diffusion process. Novel parallelized computer codes allow quantifying the geometry of the porosity and the drainage architecture from the very large tomographic datasets (2048 3 voxels) in unprecedented detail. We determined position, volume, shape and orientation of each resolvable pore and tracked these properties over the duration of the experiment. We found that the pore-size distribution follows a power law. Pores tend to be anisotropic but rarely crack-shaped and have a preferred orientation, likely controlled by a pre-existing fabric in the sample. With on-going dehydration, pores coalesce into a single interconnected pore cluster that is connected to the surface of the sample cylinder and provides an effective drainage pathway. Our observations can be summarized in a model in which gypsum is stabilized by thermal expansion stresses and locally increased pore fluid pressures until the dehydration front approaches to within about 100 μm. Then, the internal stresses are released and dehydration happens efficiently, resulting in new pore space. Pressure release, the production of pores and the advance of the front are coupled in a feedback loop.

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: The regulation of the air: a hypothesis Solid Earth, 3, 87-96, 2012 Author(s): E. G. Nisbet, C. M. R. Fowler, and R. E. R. Nisbet We propose the hypothesis that natural selection, acting on the specificity or preference for CO 2 over O 2 of the enzyme rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase), has controlled the CO 2 :O 2 ratio of the atmosphere since the evolution of photosynthesis and has also sustained the Earth's greenhouse-set surface temperature. Rubisco works in partnership with the nitrogen-fixing enzyme nitrogenase to control atmospheric pressure. Together, these two enzymes control global surface temperature and indirectly the pH and oxygenation of the ocean. Thus, the co-evolution of these two enzymes may have produced clement conditions on the Earth's surface, allowing life to be sustained.

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: Constraining fault interpretation through tomographic velocity gradients: application to northern Cascadia Solid Earth, 3, 53-61, 2012 Author(s): K. Ramachandran Spatial gradients of tomographic velocities are seldom used in interpretation of subsurface fault structures. This study shows that spatial velocity gradients can be used effectively in identifying subsurface discontinuities in the horizontal and vertical directions. Three-dimensional velocity models constructed through tomographic inversion of active source and/or earthquake traveltime data are generally built from an initial 1-D velocity model that varies only with depth. Regularized tomographic inversion algorithms impose constraints on the roughness of the model that help to stabilize the inversion process. Final velocity models obtained from regularized tomographic inversions have smooth three-dimensional structures that are required by the data. Final velocity models are usually analyzed and interpreted either as a perturbation velocity model or as an absolute velocity model. Compared to perturbation velocity model, absolute velocity models have an advantage of providing constraints on lithology. Both velocity models lack the ability to provide sharp constraints on subsurface faults. An interpretational approach utilizing spatial velocity gradients applied to northern Cascadia shows that subsurface faults that are not clearly interpretable from velocity model plots can be identified by sharp contrasts in velocity gradient plots. This interpretation resulted in inferring the locations of the Tacoma, Seattle, Southern Whidbey Island, and Darrington Devil's Mountain faults much more clearly. The Coast Range Boundary fault, previously hypothesized on the basis of sedimentological and tectonic observations, is inferred clearly from the gradient plots. Many of the fault locations imaged from gradient data correlate with earthquake hypocenters, indicating their seismogenic nature.

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: 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: Triplicated P-wave measurements for waveform tomography of the mantle transition zone Solid Earth, 3, 339-354, 2012 Author(s): S. C. Stähler, K. Sigloch, and T. Nissen-Meyer Triplicated body waves sample the mantle transition zone more extensively than any other wave type, and interact strongly with the discontinuities at 410 km and 660 km. Since the seismograms bear a strong imprint of these geodynamically interesting features, it is highly desirable to invert them for structure of the transition zone. This has rarely been attempted, due to a mismatch between the complex and band-limited data and the (ray-theoretical) modelling methods. Here we present a data processing and modelling strategy to harness such broadband seismograms for finite-frequency tomography. We include triplicated P-waves (epicentral distance range between 14 and 30°) across their entire broadband frequency range, for both deep and shallow sources. We show that is it possible to predict the complex sequence of arrivals in these seismograms, but only after a careful effort to estimate source time functions and other source parameters from data, variables that strongly influence the waveforms. Modelled and observed waveforms then yield decent cross-correlation fits, from which we measure finite-frequency traveltime anomalies. We discuss two such data sets, for North America and Europe, and conclude that their signal quality and azimuthal coverage should be adequate for tomographic inversion. In order to compute sensitivity kernels at the pertinent high body wave frequencies, we use fully numerical forward modelling of the seismic wavefield through a spherically symmetric Earth.

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: 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: Heat-flow and subsurface temperature history at the site of Saraya (eastern Senegal) Solid Earth, 3, 213-224, 2012 Author(s): F. Lucazeau and F. Rolandone New temperature measurements from eight boreholes in the West African Craton (WAC) reveal superficial perturbations down to 100 m below the alteration zone. These perturbations are both related to a recent increase in the surface air temperature (SAT) and to the site effects caused by fluid circulations and/or the lower conduction in the alterites. The ground surface temperature (GST), inverted from the boreholes temperatures, increased slowly in the past (~0.4 °C from 1700 to 1940) and then, more importantly, in recent years (~1.5 °C from 1940 to 2010). This recent trend is consistent with the increase of the SAT recorded at two nearby meteorological stations (Tambacounda and Kedougou), and more generally in the Sahel with a coeval rainfall decrease. Site effects are superimposed to the climatic effect and interpreted by advective (circulation of fluids) or conductive (lower conductivity of laterite and of high-porosity sand) perturbations. We used a 1-D finite differences thermal model and a Monte-Carlo procedure to find the best estimates of these site perturbations: all the eight borehole temperature logs can be interpreted with the same basal heat-flow and the same surface temperature history, but with some realistic changes of thermal conductivity and/or fluid velocity. The GST trend observed in Senegal can be confirmed by two previous borehole measurements made in 1983 in other locations of West Africa, the first one in an arid zone of northern Mali and the second one in a sub-humid zone in southern Mali. Finally, the background heat-flow is low (31±2 mW m −2 ), which makes this part of the WAC more similar with the observations in the southern part (33±8 mW m −2 ) rather than with those in the northern part and in the Pan-African domains where the surface heat-flow is 15–20 mW m −2 higher.

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: 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: Plate kinematics in the Cantabrian domain of the Pyrenean orogen Solid Earth, 3, 265-292, 2012 Author(s): S. Tavani The Cantabrian domain represents the western portion of the Pyrenean orogen, in the area where the Iberian continental lithosphere was subducted toward the north underneath the transitional to oceanic lithosphere of the Bay of Biscay. There, the about 100 km of orogenic convergence have been mostly accommodated in the northern portion of the orogen (i.e. the retro wedge) developed in the Bay of Biscay abyssal plain, while only crustal-scale folding with limited internal deformation occurred in the Cantabrian southern wedge (pro-wedge). Integrated meso- and macrostructural analyses and a reappraisal of available information from the transitional area between the Pyrenean and Cantabrian domains are presented in this work, allowing to set geometric and kinematic constraints on the entire Meso-Cenozoic history of the northern portion of the Iberian Plate, including subduction initiation and evolution in the western portion of the Pyrenean orogen. The structural record of the Late Jurassic to Early Cretaceous deformation stage, which was associated with rifting and seafloor spreading in the Bay of Biscay, indicates a ridge perpendicular (NNE-SSW oriented) extension, with no evidence of relevant strike-slip components during rifting. A Cenozoic NNW-SSE oriented shortening stage followed, related to the limited (about 100 km) north-directed subduction of the Iberian continental lithosphere underneath the transitional to oceanic lithosphere of the Bay of Biscay. Subduction led to the formation of the poorly-developed Cantabrian pro-wedge, which is laterally juxtaposed to the well-developed Pyrenean pro-wedge to the east. During this convergence stage, the structural framework in the Cantabrian pro-wedge, and particularly along its transition with the Pyrenean wedge to the east, was severely complicated by the reactivation of Paleozoic and Mesozoic inherited structures. Data presented in this work fully support the development of the Cantabrian Mountains as related to indentation and consequent thickening of the Bay of Biscay transitional lower crust during north-directed subduction of Iberian continental lithosphere. In essence, the Cantabrian pro-wedge is a lithospheric south-verging fault-propagation anticline developing above the subduction plane. The structural record in the area indicates that a lithospheric fault-propagation folding stage was predated, during the very early stages of orogenic shortening, by the development of a lithospheric-scale open syncline overlying the nucleation point of lithosphere sinking. Such a syncline is today partially preserved and represents one of the few natural examples of subduction initiation.

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: Seismicity from February 2006 to September 2007 at the Rwenzori Mountains, East African Rift: earthquake distribution, magnitudes and source mechanisms Solid Earth, 3, 251-264, 2012 Author(s): M. Lindenfeld, G. Rümpker, A. Batte, and A. Schumann We have analysed the microseismic activity within the Rwenzori Mountains area in the western branch of the East African Rift. Seismogram recordings from a temporary array of up to 27 stations reveal approximately 800 events per month with local magnitudes ranging from –0.5 to 5.1. The earthquake distribution is highly heterogeneous. The majority of located events lie within faults zones to the east and west of the Rwenzoris with the highest seismic activity observed in the northeastern area, where the mountains are in contact with the rift shoulders. The hypocentral depth distribution exhibits a pronounced peak of seismic energy release at 15 km depth. The maximum extent of seismicity ranges from 20 to 32 km and correlates well with Moho depths that were derived from teleseismic receiver functions. We observe two general features: (i) beneath the rift shoulders, seismicity extends from the surface down to ca. 30 km depth; (ii) beneath the rift valley, seismicity is confined to depths greater than 10 km. From the observations there is no indication for a crustal root beneath the Rwenzori Mountains. The magnitude frequency distribution reveals a b -value of 1.1, which is consistent with the hypothesis that part of the seismicity is caused by magmatic processes within the crust. Fault plane solutions of 304 events were derived from P-polarities and SV/P amplitude ratios. More than 70 % of the source mechanisms exhibit pure or predominantly normal faulting. T-axis trends are highly uniform and oriented WNW–ESE, which is perpendicular to the rift axis and in good agreement with kinematic rift models. At the northernmost part of the region we observe a rotation of the T-axis trends to NEN–SWS, which may be indicative of a local perturbation of the regional stress field.

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: Shallow water carbonate platforms (Late Aptian–Early Albian, Southern Apennines) in the context of supraregional to global changes: re-appraisal of palaeoecological events as reflectors of carbonate factory response Solid Earth, 3, 225-249, 2012 Author(s): A. Raspini This paper discusses the palaeoenvironmental significance of the "Orbitolina Level", the microbial carbonates and the Salpingoporella dinarica -rich deposits encased in the Aptian/Albian shallow water carbonate platform strata of Monte Tobenna and Monte Faito (Southern Italy). These facies show a peculiar field appearance due to their color and/or fossil content. In the shallow water carbonate strata, the Late Aptian "Orbitolina Level" was formed during a period of decreasing accommodation space. Microbial carbonates occur in different levels in the composite section. They reach their maximum thickness around the sequence boundaries just above the "Orbitolina Level" and close to the Aptian–Albian transition, and were not deposited during maximum flooding. S. dinarica -rich deposits occur in the lower part of the Monte Tobenna-Monte Faito composite section, in both restricted and more open lagoonal sediments. S. dinarica has its maximum abundance below the "Orbitolina Level" and disappears 11 m above this layer. On the basis of δ 13 C and δ 18 O values recorded at Tobenna-Faito, the succession has been correlated to global sea-level changes and to the main volcanic and climatic events during the Aptian. Deterioration of the inner lagoon environmental conditions was related to high trophic levels triggered by volcano-tectonic activity. Microbial carbonates were deposited especially in periods of third-order sea level lowering. In such a scenario, periods of increased precipitation during the Gargasian induced the mobilization of clay during flooding of the exposed platform due to high-frequency sea-level changes, with consequent terrigenous input to the lagoon. This and the high nutrient levels made the conditions unsuitable for the principle carbonate producers, and an opportunistic biota rich in orbitolinids ( Mesorbitolina texana and M. parva ) populated the platform. In the more open marine domain, the increased nutrient input enhanced the production of organic matter and locally led to the formation of black shales (e.g. the Niveau Fallot in the Vocontian Basin). It is argued that the concomitant low Mg/Ca molar ratio and high concentration of calcium in seawater could have favoured the development of the low-Mg calcite skeleton of the S. dinarica green algae. During third-order sea-level rise, no or minor microbial carbonates formed in the shallowlagoonal settings and S. dinarica disappeared. Carbonate neritic ecosystems were not influenced by the environmental changes inferred to have been induced by the mid-Cretaceous volcanism. The "Orbitolina Level", the microbial carbonates and the Salpingoporella dinarica -rich deposits in the studied Aptian/Albian shallow water carbonate strata are interpreted to be the response to environmental and oceanographic changes in shallow-water and deeper-marine ecosystems.

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: High resolution reflection seismic profiling over the Tjellefonna fault in the Møre-Trøndelag Fault Complex, Norway Solid Earth, 3, 175-188, 2012 Author(s): E. Lundberg, C. Juhlin, and A. Nasuti The Møre-Trøndelag Fault Complex (MTFC) is one of the most prominent fault zones of Norway, both onshore and offshore. In spite of its importance, very little is known of the deeper structure of the individual fault segments comprising the fault complex. Most seismic lines have been recorded offshore or focused on deeper structures. This paper presents results from two reflection seismic profiles, located on each side of the Tingvollfjord, acquired over the Tjellefonna fault in the southeastern part of the MTFC. Possible kilometer scale vertical offsets, reflecting large scale northwest-dipping normal faulting, separating the high topography to the southeast from lower topography to the northwest have been proposed for the Tjellefonna fault or the Baeverdalen lineament. In this study, however, the Tjellefonna fault is interpreted to dip approximately 50–60° towards the southeast to depths of at least 1.3 km. Travel-time modeling of reflections associated with the fault was used to establish the geometry of the fault structure at depth, while detailed analysis of first P-wave arrivals in shot gathers, together with resistivity profiles, were used to define the near surface geometry of the fault zone. A continuation of the structure on the northeastern side of the Tingvollfjord is suggested by correlation of an in strike direction P-S converted reflection (generated by a fracture zone) seen on the reflection data from that side of the Tingvollfjord. The reflection seismic data correlate well with resistivity profiles and recently published near surface geophysical data. A highly reflective package forming a gentle antiform structure was also identified on both seismic profiles. This structure could be related to the folded amphibolite lenses seen on the surface or possibly by an important boundary within the gneissic basement rocks of the Western Gneiss Region. The fold hinge line of the structure is parallel with the Tjellefonna fault trace suggesting that the folding and faulting may have been related.

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.

Description: Modeling atmospheric ammonia and ammonium using a backward-in-time stochastic Lagrangian air quality model (STILT-Chem v0.7) Geoscientific Model Development Discussions, 5, 2745-2788, 2012 Author(s): D. Wen, J. C. Lin, L. Zhang, R. Vet, and M. D. Moran A new chemistry module of atmospheric ammonia (NH 3 ) and ammonium (NH 4 + ) was incorporated into a backward-in-time stochastic Lagrangian air quality model (STILT-Chem) that was originally developed to simulate the concentrations of a variety of gas-phase species at receptors. STILT-Chem simulates the transport of air parcels backward in time using ensembles of fictitious particles with stochastic motions, while simulating emissions, deposition and chemical transformation forward in time along trajectories identified by the backward-in-time simulations. The incorporation of the new chemistry module allows the model to simulate not only gaseous species, but also multi-phase species involving NH 3 and NH 4 + . The model was applied to simulate concentrations of NH 3 and particulate NH 4 + at six sites in the Canadian province of Ontario for a six-month period in 2006. The model-predicted concentrations of NH 3 and particulate NH 4 + were compared with observations, which show broad agreement between simulated concentrations and observations. Since the model is based on back trajectories, the influence of each major process such as emission, deposition and chemical conversion on the concentration of a modeled species at a receptor can be determined for every upstream location at each time step. This makes it possible to quantitatively investigate the upstream processes affecting receptor concentrations. The modeled results suggest that the concentrations of NH 3 at those sites were significantly and frequently affected by southwestern Ontario, northern Ohio, and nearby areas. NH 3 is mainly contributed by emission sources whereas particulate NH 4 + is mainly contributed by the gas-to-aerosol chemical conversion of NH 3 . Dry deposition is the largest removal process for both NH 3 and particulate NH 4 + . This study revealed the contrast between agricultural versus forest sites. Not only were emissions of NH 3 higher, but removal mechanisms (especially chemical loss for NH 3 and dry deposition for NH 4 + ) were less efficient for agricultural sites. This combination explains the significantly higher concentrations of NH 3 and particulate NH 4 + observed at agricultural sites.

Description: How should sparse in situ measurements be compared to continuous model data? Geoscientific Model Development Discussions, 5, 2311-2345, 2012 Author(s): L. de Mora, M. Butenschön, and J. I. Allen This work demonstrates the importance of an adequate method to sub-sample model results when comparing with in situ measurements. A test of model skill was performed by comparing a multi-decadal hindcast against a sparse, unevenly distributed historic in situ dataset. The comparison was performed using a point-to-point method. The point-to-point method masked out all hindcast cells that did not have a corresponding in situ measurement in order to compare each in situ measurement against its most similar cell from the model. The application of the point-to-point method showed that the model was successful at reproducing many inter-annual trends. Furthermore, this success was not immediately apparent using the previous comparison methods, which compared model and measurements aggregated to regional averages. Time series, data density and target diagrams were employed to illustrate the impact of switching from the previous method to the point-to-point method. The comparison based on regional averages gave significantly different and sometimes contradicting results that could lead to erroneous conclusions on the model performance. We therefore recommend that researchers take into account for the limitations of the in situ datasets, process the model to resemble the data as much as possible, and we advocate greater transparency in the publication of methodology.

Description: The chemical transport model Oslo CTM3 Geoscientific Model Development Discussions, 5, 1561-1626, 2012 Author(s): O. A. Søvde, M. J. Prather, I. S. A. Isaksen, T. K. Berntsen, F. Stordal, X. Zhu, C. D. Holmes, and J. Hsu We present here the global chemical transport model Oslo CTM3, an update of the Oslo CTM2. The update comprises a faster transport scheme, an improved wet scavenging scheme for large scale rain, updated photolysis rates and a new lightning parameterization. Oslo CTM3 is better parallelized and allows for stable, large time steps for advection, enabling more complex or high resolution simulations. Thorough comparisons between the Oslo CTM3, Oslo CTM2 and measurements are performed, and in general the Oslo CTM3 is found to reproduce measurements well. Inclusion of tropospheric sulfur chemistry and nitrate aerosols in CTM3 is shown to be important to reproduce tropospheric O 3 , OH and the CH 4 lifetime well. Using the same meteorology to drive the two models, shows that some features related to transport are better resolved by the CTM3, such as polar cap transport, while features like transport close to the vortex edge are resolved better in the Oslo CTM2 due to its required shorter transport time step. The longer transport time steps in CTM3 result in larger errors e.g. near the jets, and when necessary, this can be remedied by using a shorter time step. An additional, more accurate and time consuming, treatment of polar cap transport is presented, however, both perform acceptably. A new treatment of the horizontal distribution of lightning is presented and found to compare well with measurements. Vertical distributions of lighting are updated, and tested against the old vertical distribution. The new profiles are found to produce more NO x in the tropical middle troposphere, and less at the surface and at high altitudes.

Description: Corrigendum to "Floating stones off El Hierro, Canary Islands: xenoliths of pre-island sedimentary origin in the early products of the October 2011 eruption" published in Solid Earth, 3, 97–110, 2012 Solid Earth, 3, 189-189, 2012 Author(s): V. R. Troll, A. Klügel, M.-A. Longpré, S. Burchardt, F. M. Deegan, J. C. Carracedo, S. Wiesmaier, U. Kueppers, B. Dahren, L. S. Blythe, T. H. Hansteen, C. Freda, D. A. Budd, E. M. Jolis, E. Jonsson, F. C. Meade, C. Harris, S. E. Berg, L. Mancini, M. Polacci, and K. Pedroza No abstract available.

Description: Unified parameterization of the planetary boundary layer and shallow convection with a higher-order turbulence closure in the community atmosphere model: single column experiments Geoscientific Model Development Discussions, 5, 1743-1780, 2012 Author(s): P. A. Bogenschutz, A. Gettelman, H. Morrison, V. E. Larson, D. P. Schanen, N. R. Meyer, and C. Craig This paper describes the coupling of the Community Atmosphere Model (CAM) version 5 with a unified multi-variate probability density function (PDF) parameterization, Cloud Layers Unified by Binormals (CLUBB). CLUBB replaces the planetary boundary layer (PBL), shallow convection, and cloud macrophysics schemes in CAM5 with a higher-order turbulence closure based on an assumed PDF. Comparisons of single-column versions of CAM5 and CAM-CLUBB are provided in this paper for several boundary layer regimes. As compared to Large Eddy Simulations (LES), CAM-CLUBB and CAM5 simulate marine stratocumulus regimes with similar accuracy. For shallow convective regimes, CAM-CLUBB improves the representation of cloud cover and liquid water path (LWP). In addition, for shallow convection CAM-CLUBB offers better fidelity for sub-grid scale vertical velocity, which is an important input for aerosol activation. Finally, CAM-CLUBB results are more robust to changes in vertical and temporal resolution when compared to CAM5.

Description: Implementation of the chemistry module MECCA (v2.5) in the modal aerosol version of the Community Atmosphere Model component (v3.6.33) of the Community Earth System Model Geoscientific Model Development Discussions, 5, 1483-1501, 2012 Author(s): M. S. Long, W. C. Keene, R. Easter, R. Sander, A. Kerkweg, D. Erickson, X. Liu, and S. Ghan A coupled atmospheric chemistry and climate system model was developed using the modal aerosol version of the National Center for Atmospheric Research Community Atmosphere Model (modal-CAM) and the Max Planck Institute for Chemistry's Module Efficiently Calculating the Chemistry of the Atmosphere (MECCA) to provide enhanced resolution of multiphase processes, particularly those involving inorganic halogens, and associated impacts on atmospheric composition and climate. Three Rosenbrock solvers (Ros-2, Ros-3, RODAS-3) were tested in conjunction with the basic load balancing options available to modal CAM (1) to establish an optimal configuration of the implicitly-solved multiphase chemistry module that maximizes both computational speed and repeatability of Ros-2 and RODAS-3 results versus Ros-3, and (2) to identify potential implementation strategies for future versions of this and similar coupled systems. RODAS-3 was faster than Ros-2 and Ros-3 with good reproduction of Ros-3 results, while Ros-2 was both slower and substantially less reproducible relative to Ros-3 results. Modal-CAM with MECCA chemistry was a factor of 15 slower than modal-CAM using standard chemistry. MECCA chemistry integration times demonstrated a systematic frequency distribution for all three solvers, and revealed that the change in run-time performance was due to a change in the frequency distribution chemical integration times; the peak frequency was similar for all solvers. This suggests that efficient chemistry-focused load-balancing schemes can be developed that rely on the parameters of this frequency distribution.

Description: Describing Earth System Simulations with the Metafor CIM Geoscientific Model Development Discussions, 5, 1669-1689, 2012 Author(s): B. N. Lawrence, V. Balaji, P. Bentley, S. Callaghan, C. DeLuca, S. Denvil, G. Devine, M. Elkington, R. W. Ford, E. Guilyardi, M. Lautenschlager, M. Morgan, M.-P. Moine, S. Murphy, C. Pascoe, H. Ramthun, P. Slavin, L. Steenman-Clark, F. Toussaint, A. Treshansky, and S. Valcke The Metafor project has developed a Common Information Model (CIM) using the ISO1900 series formalism to describe the sorts of numerical experiments carried out by the earth system modelling community, the models they use, and the simulations that result. Here we describe the mechanism by which the CIM was developed, and its key properties. We introduce the conceptual and application versions and the controlled vocabularies developed in the context of supporting the fifth Coupled Model Intercomparison Project (CMIP5). We describe how the CIM has been used in experiments to describe model coupling properties and describe the near term expected evolution of the CIM.

Description: Lidar signal simulation for the evaluation of aerosols in chemistry-transport models Geoscientific Model Development Discussions, 5, 1691-1741, 2012 Author(s): S. Stromatas, S. Turquety, L. Menut, H. Chepfer, J. C. Péré, G. Cesana, and B. Bessagnet We present an adaptable tool, the OPTSIM (OPTical properties SIMulation) software, for the simulation of optical properties and lidar attenuated backscattered profiles (β ' ) from aerosol concentrations calculated by chemistry-transport models (CTM). It was developed to support model evaluation using an original approach based on the lidar Level 1 observations (β ' ), avoiding the use of Level 2 aerosol retrievals which include specific assumptions on aerosol types that may not be in agreement with the CTM. In addition to an evaluation of the aerosol loading and optical properties, active remote sensing allows the analysis of aerosols' vertical structures. An academic case study for two different species (black carbon and dust) is presented and shows the consistency of the simulator. Illustrations are then given through the analysis of dust events in the Mediterranean region during the summer 2007. These are based on simulations by the CHIMERE regional CTM and observations from the CALIOP space-based lidar, and highlight the potential of this approach to evaluate the concentration, size and vertical structure of the aerosol plumes.

Description: Better constraints on the sea-ice state using global sea-ice data assimilation Geoscientific Model Development Discussions, 5, 1627-1667, 2012 Author(s): P. Mathiot, C. König Beatty, T. Fichefet, H. Goosse, F. Massonnet, and M. Vancoppenolle Short-term and decadal sea-ice prediction systems need a realistic initial state, generally obtained using ice-ocean model simulations with data assimilation. However, only sea-ice concentration and velocity data are currently assimilated. In this work, an Ensemble Kalman Filter system is used to assimilate observed ice concentration and freeboard (i.e. thickness of emerged sea ice) data into a global coupled ocean–sea-ice model. The impact and effectiveness of our data assimilation system is assessed in two steps: firstly, through the assimilation of synthetic data (i.e., model-generated data) and, secondly, through the assimilation of satellite data. While ice concentrations are available daily, freeboard data used in this study are only available during six one-month periods spread over 2005–2007. Our results show that the simulated Arctic and Antarctic sea-ice extents are improved by the assimilation of synthetic ice concentration data. Assimilation of synthetic ice freeboard data improves the simulated sea-ice thickness field. Using real ice concentration data enhances the model realism in both hemispheres. Assimilation of ice concentration data significantly improves the total hemispheric sea-ice extent all year long, especially in summer. Combining the assimilation of ice freeboard and concentration data leads to better ice thickness, but does not further improve the ice extent. Moreover, the improvements in sea-ice thickness due to the assimilation of ice freeboard remain visible well beyond the assimilation periods.

Description: Modeling wet deposition of inorganics over Northeast Asia with MRI-PM/c and the effects of super large sea salt droplets at near-the-coast stations Geoscientific Model Development Discussions, 5, 1341-1379, 2012 Author(s): M. Kajino, M. Deushi, T. Maki, N. Oshima, Y. Inomata, K. Sato, T. Ohizumi, and H. Ueda We conducted a regional-scale simulation (with grid spacing = 60 km) over Northeast Asia for the entire year of 2006 by using an aerosol chemical transport model, the lateral and upper boundary concentrations of which we predicted with a global stratospheric and tropospheric chemistry-climate model, with a horizontal resolution of T42 (grid spacing ~300 km) and a time resolution of 1 h. The present one-way nested global-through-regional-scale model is called the Meteorological Research Institute – Passive-tracers Model system for atmospheric Chemistry (MRI-PM/c). We evaluated the model performance with respect to the major inorganic components in rain and snow measured by stations of the Acid Deposition Monitoring Network in East Asia (EANET). Through statistical analysis, we show that the model successfully reproduced the regional-scale processes of emission, transport, transformation, and wet deposition of major inorganic species derived from anthropogenic and natural sources, including SO 4 2− , NH 4 + , NO 3 − , Na + and Ca 2+ . Interestingly, the only exception was Na + in precipitation at near-coastal stations (where the distance from the coast was from 150 to 700 m), concentrations of which were significantly underestimated by the model, by up to a factor of 30. This result suggested that the contribution of short-lived, super-large sea salt droplets (SLSD; D 〉 10–100 μm) was substantial in precipitation samples at stations near the coast of Japan; thus samples were horizontally representative only within the traveling distances of SLSD (from 1 to 10 km). Nevertheless, the calculated effect of SLSD on precipitation pH was very low, a change of about +0.014 on average, even if the ratio of SLSD to all sea salt in precipitation was assumed to be 90%.

Description: New developments in the analysis of column-collapse pyroclastic density currents through numerical simulations of multiphase flows Solid Earth, 3, 161-173, 2012 Author(s): S. Lepore and C. Scarpati A granular multiphase model has been used to evaluate the action of differently sized particles on the dynamics of fountains and associated pyroclastic density currents. The model takes into account the overall disequilibrium conditions between a gas phase and several solid phases, each characterized by its own physical properties. The dynamics of the granular flows (fountains and pyroclastic density currents) has been simulated by adopting a Reynolds-averaged Navier-Stokes model for describing the turbulence effects. Numerical simulations have been carried out by using different values for the eruptive column temperature at the vent, solid particle frictional concentration, turbulent kinetic energy, and dissipation. The results obtained provide evidence of the multiphase nature of the model and describe several disequilibrium effects. The low concentration (≤5 × 10 −4 ) zones lie in the upper part of the granular flow, above the fountain, and above the tail and body of pyroclastic density current as thermal plumes. The high concentration zones, on the contrary, lie in the fountain and at the base of the current. Hence, pyroclastic density currents are assimilated to granular flows constituted by a low concentration suspension flowing above a high concentration basal layer (boundary layer), from the proximal regions to the distal ones. Interactions among the solid particles in the boundary layer of the granular flow are controlled by collisions between particles, whereas the dispersal of particles in the suspension is determined by the dragging of the gas phase. The simulations describe well the dynamics of a tractive boundary layer leading to the formation of stratified facies during Strombolian to Plinian eruptions.

Description: Mid-Pliocene climate modelled using the UK Hadley Centre Model: PlioMIP Experiments 1 and 2 Geoscientific Model Development Discussions, 5, 837-871, 2012 Author(s): F. J. Bragg, D. J. Lunt, and A. M. Haywood The Pliocene Model Intercomparison Project (PlioMIP) project is a sub-project of the Paleoclimate Modelling Intercomparison Project (PMIP) whose objective is to compare predictions of the mid-Pliocene climate from the widest possible range of general circulation models. The mid-Pliocene (3.3–3.0 Ma) is the most recent sustained period of greater warmth and atmospheric carbon dioxide concentration than the pre-industrial times and as such has potential to inform predictions of our warming climate in the coming century. This paper describes the UK contribution to PlioMIP using the Hadley Centre Model both in atmosphere-only mode (HadAM3, PlioMIP Experiment 1) and atmosphere-ocean coupled mode (HadCM3, PlioMIP Experiment 2). The coupled model predicts a greater overall warming (3.3 °C) relative to the control than the atmosphere-only (2.5 °C). The Northern Hemisphere latitudinal temperature gradient is greater in the coupled model with a warmer equator and colder Arctic than the atmosphere-only model, which is constrained by sea surface temperatures from Pliocene proxy reconstructions. The atmosphere-only model predicts a reduction in equatorial precipitation and south Asian monsoon intensity whereas the coupled models shows and increase in the intensity of these systems. Sensitivity studies using alternative boundary conditions for both the Pliocene and the control simulations are presented, which indicate the sensitivity of the mid-Pliocene warming to uncertainties in both pre-industrial and mid-Pliocene climate.

Description: The fate of fluids released from subducting slab in northern Cascadia Solid Earth, 3, 121-129, 2012 Author(s): K. Ramachandran and R. D. Hyndman Large amounts of water carried down in subduction zones are driven upward into the overlying forearc upper mantle and crust as increasing temperatures and pressure dehydrate the subducting crust. Through seismic tomography velocities we show (a) the overlying forearc mantle in northern Cascadia is hydrated to serpentinite, and (b) there is low Poisson's ratio at the base of the forearc lower crust that may represent silica deposited from the rising fluids. From the velocities observed in the forearc mantle, the volume of serpentinite estimated is ∼30 %. This mechanically weak hydrated forearc region has important consequences in limits to great earthquakes and to collision tectonics. An approximately 10 km thick lower crustal layer of low Poisson's ratio (σ = 0.22) in the forearc is estimated to represent a maximum addition of ∼14 % by volume of quartz (σ = 0.09). If this quartz is removed from rising silica-saturated fluids over long times, it represents a significant addition of silica to the continental crust and an important contributor to its average composition.

Description: Performance of McRAS-AC in the GEOS-5 AGCM: aerosol-cloud-microphysics, precipitation, cloud radiative effects, and circulation Geoscientific Model Development Discussions, 5, 1381-1434, 2012 Author(s): Y. C. Sud, D. Lee, L. Oreopoulos, D. Barahona, A. Nenes, and M. J. Suarez A revised version of the Microphysics of clouds with Relaxed Arakawa-Schubert and Aerosol-Cloud interaction scheme (McRAS-AC) including, among others, the Barahona and Nenes ice nucleation parameterization, is implemented in the GEOS-5 AGCM. Various fields from a 10-yr long integration of the AGCM with McRAS-AC were compared with their counterparts from an integration of the baseline GEOS-5 AGCM using satellite data as observations. Generally McRAS-AC simulations have smaller biases in cloud fields and cloud radiative effects over most of the regions of the Earth than the baseline GEOS-5 AGCM. Two systematic biases are identified in the McRAS-AC runs: one under-prediction of cloud particles around 40° S–60° S, and one over-prediction of cloud water path during Northern Hemisphere summer over the Gulf Stream and North Pacific. Sensitivity analyses show that these biases potentially originate from biases in the aerosol input. The first bias is largely eliminated in a sensitivity test using 50% smaller sea-salt aerosol particles, while the second bias is much reduced when interactive aerosol chemistry was turned on. The main drawback of McRAS-AC is dearth of low-level marine stratus clouds, probably due to lack of boundary-layer clouds that is an outcome of explicit dry-convection not yet implemented into the cloud model. Nevertheless, McRAS-AC simulates realistic clouds and their optical properties that can further improve with better aerosol-input. Thereby, McRAS-AC has the potential to be a valuable tool for climate modeling research because of its superior simulation capabilities that physically couple aerosols, cloud microphysics, cloud macrophysics, and cloud-radiation interaction for all clouds.

Description: Modelling mid-Pliocene climate with COSMOS Geoscientific Model Development Discussions, 5, 917-966, 2012 Author(s): C. Stepanek and G. Lohmann In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized community earth system models (COSMOS) and document the procedures which we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping Project (PRISM) mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo and preindustrial (PI) time-slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene as simulated with our COSMOS-setup and PRISM boundary conditions is both warmer and wetter than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.

Description: TopoSUB: a tool for efficient large area numerical modelling in complex topography at sub-grid scales Geoscientific Model Development Discussions, 5, 1041-1076, 2012 Author(s): J. Fiddes and S. Gruber Mountain regions are highly sensitive to global climate change. However, large scale assessments of mountain environments remain problematic due to the high resolution required of model grids to capture strong lateral variability. To alleviate this, tools are required to bridge the scale gap between gridded climate datasets (climate models and re-analyses) and unresolved (by coarse grids) sub-grid mountain topography. We address this problem with a sub-grid method. It relies on sampling the most important aspects of land surface heterogeneity through a lumped scheme, allowing for the application of numerical land surface models (LSM) over large areas in mountain regions. This is achieved by including the effect of mountain topography on these processes at the sub-grid scale using a multidimensional informed sampling procedure together with a 1-D lumped model that can be driven by gridded climate datasets. This paper provides a description of this sub-grid scheme, TopoSUB, as well as assessing its performance against a distributed model. We demonstrate the ability of TopoSUB to approximate results simulated by a distributed numerical LSM at around 10 4 less computations. These significant gains in computing resources allow for: (1) numerical modelling of processes at fine grid resolutions over large areas; (2) extremely efficient statistical descriptions of sub-grid behaviour; (3) a "sub-grid aware" aggregation of simulated variables to course grids; and (4) freeing of resources for treatment of uncertainty in the modelling process.

Description: Impact of a time-dependent background error covariance matrix on air quality analysis Geoscientific Model Development Discussions, 5, 873-916, 2012 Author(s): E. Jaumouillé, S. Massart, A. Piacentini, D. Cariolle, and V.-H. Peuch In this article we study the influence of different characteristics of our assimilation system on the surface ozone analyses over Europe. Emphasis is placed on the evaluation of the background error covariance matrix (BECM). Data assimilation systems require a BECM in order to obtain an optimal representation of the physical state. A posteriori diagnostics are an efficient way to check the consistency of the used BECM. In this study we derived a diagnostic to estimate the BECM. On the other hand an increasingly used approach to obtain such a covariance matrix is to estimate it from an ensemble of perturbed assimilation experiments. We applied this method, combined with variational assimilation, while analysing the surface ozone distribution over Europe. We first show that the resulting covariance matrix is strongly time (hourly and seasonally) and space dependent. We then built several configurations of the background error covariance matrix with none, one or two of its components derived from the ensemble estimation. We used each of these configurations to produce surface ozone analyses. All the analyses are compared between themselves and compared to assimilated data or data from independent validation stations. The configurations are very well correlated with the validation stations, but with varying regional and seasonal characteristics. The largest correlation is obtained with the experiments using time and space dependent correlation of the background errors. Results show that our assimilation process is efficient in bringing the model assimilations closer to the observations than the direct simulation, but we cannot conclude which BECM configuration is the best. The impact of the background error covariances configuration on four-days forecasts is also studied. Although mostly positive, the impact depends on the season and lasts longer during the winter season.

Description: Numerical uncertainty at mesoscale in a Lagrangian model in complex terrain Geoscientific Model Development Discussions, 5, 967-991, 2012 Author(s): J. Brioude, W. M. Angevine, S. A. McKeen, and E.-Y. Hsie Recently, it has been shown that mass conservation in Lagrangian models is improved by using time-average winds out of Eulerian models. In this study, we evaluate the mass conservation and trajectory uncertainties in complex terrain at mesoscale using the FLEXPART Lagrangian particle dispersion model coupled with the WRF mesoscale model. The specific form of vertical wind used is found to have a large effect. Time average wind with time average sigma dot (σ · ), instantaneous wind with geometric cartesian vertical wind ( w ) and instantaneous wind with σ · are used to simulate mixing ratios of a passive tracer in forward and backward runs using different time interval outputs and horizontal resolutions in California. Mass conservation in the FLEXPART model was not an issue when using time-average wind or instantaneous wind with σ · . However, mass was poorly conserved using instantaneous wind with w , with a typical variation of 25% within 24 h. Uncertainties in surface residence time (a backtrajectory product commonly used in source receptor studies or inverse modeling) calculated for each backtrajectory run were also analyzed. The smallest uncertainties were systematically found when using time-average wind. Uncertainties using instantaneous wind with σ · were slightly larger, as long as the time interval of output was sufficiently small. The largest uncertainties were found when using instantaneous wind with w . Those uncertainties were found to be linearly correlated with the local average gradient of orography. Differences in uncertainty were much smaller when trajectories were calculated over flat terrain. For a typical run at mesoscale in complex terrain, 4 km horizontal resolution and 1 h time interval output, the average uncertainty and bias in surface residence time is, respectively 8.4% and −2.5% using time-average wind, and 13% and −3.7% using instantaneous wind with σ · in complex terrain. The corresponding values for instantaneous wind with cartesian w were 24% and −11%. While the use of time-average wind systematically improves uncertainty in FLEXPART, the improvements are small, and therfore a systematic use of time-average wind in Lagrangian models is not necessarily required. Use of cartesian vertical wind in complex terrain, however, should be avoided.

Description: Models of soil organic matter decomposition: the SOILR package, version 1.0 Geoscientific Model Development Discussions, 5, 993-1039, 2012 Author(s): C. A. Sierra, M. Müller, and S. E. Trumbore Organic matter decomposition is a very important process within the Earth System because it controls the rates of mineralization of carbon and other biogeochemical elements, determining their flux to the atmosphere and the hydrosphere. SOILR is a modeling framework that contains a library of functions and tools for modeling soil organic matter decomposition under the R environment for computing. It implements a variety of model structures and tools to represent carbon storage and release from soil organic matter. In SOILR organic matter decomposition is represented as a linear system of ordinary differential equations that generalizes the structure of most compartment-based decomposition models. A variety of functions is also available to represent environmental effects on decomposition rates. This document presents the conceptual basis for the functions implemented in the package. It is complementary to the help pages released with the software.

Description: A semi-analytical solution to accelerate spin-up of a coupled carbon and nitrogen land model to steady state Geoscientific Model Development Discussions, 5, 803-836, 2012 Author(s): J. Xia, Y. Luo, Y.-P. Wang, E. Weng, and O. Hararuk The spin-up of land models to steady state of coupled carbon-nitrogen processes is computationally so costly that it becomes a~bottleneck issue for global analysis. In this study, we introduced a semi-analytical solution (SAS) for the spin-up issue. SAS is fundamentally based on the analytic solution to a set of equations that describe carbon transfers within ecosystems over time. SAS is implemented by three steps: (1) having an initial spin-up with prior pool-size values until net primary productivity (NPP) reaches steady state, (2) calculating quasi steady-state pool sizes by letting fluxes of the equations equal zero, and (3) having a final spin-up to meet the criterion of steady state. Step 2 is enabled by averaged time-varying variables over one period of repeated driving forcings. SAS was applied to both site-level and global scale spin-up of the Australian Community Atmosphere Biosphere Land Exchange (CABLE) model. For the carbon-cycle-only simulations, SAS saved 95.7% and 92.4% of computational time for site-level and global spin-up, respectively, in comparison with the traditional method. For the carbon-nitrogen-coupled simulations, SAS reduced computational cost by 84.5% and 86.6% for site-level and global spin-up, respectively. The estimated steady-state pool sizes represent the ecosystem carbon storage capacity, which was 12.1 kg C m −2 with the coupled carbon-nitrogen global model, 14.6% lower than that with the carbon-only model. The nitrogen down-regulation in modeled carbon storage is partly due to the 4.6% decrease in carbon influx (i.e., net primary productivity) and partly due to the 10.5% reduction in residence times. This steady-state analysis accelerated by the SAS method can facilitate comparative studies of structural differences in determining the ecosystem carbon storage capacity among biogeochemical models. Overall, the computational efficiency of SAS potentially permits many global analyses that are impossible with the traditional spin-up methods, such as ensemble analysis of land models against parameter variations.

Description: A community diagnostic tool for Chemistry Climate Model Validation Geoscientific Model Development Discussions, 5, 1229-1261, 2012 Author(s): A. Gettelman, V. Eyring, C. Fischer, H. Shiona, I. Cionni, M. Neish, O. Morgenstern, S. W. Wood, and Z. Li This technical note presents an overview of the Chemistry-Climate Model Validation Diagnostic (CCMVal-Diag) tool for model evaluation. The CCMVal-Diag tool is a flexible and extensible open source package that facilitates the complex evaluation of global models. Models can be compared to other models, ensemble members (simulations with the same model), and/or many types of observations. The tool can also compute quantitative performance metrics. The initial construction and application is to coupled Chemistry-Climate Models (CCMs) participating in CCMVal, but the evaluation of climate models that submitted output to the Coupled Model Intercomparison Project (CMIP) is also possible. The package has been used to assist with analysis of simulations for the 2010 WMO/UNEP Scientific Ozone Assessment and the SPARC Report on the Evaluation of CCMs. The CCMVal-Diag tool is described and examples of how it functions are presented, along with links to detailed descriptions, instructions and source code. The CCMVal-Diag tool is supporting model development as well as quantifying model improvements, both for different versions of individual models and for different generations of community-wide collections of models used in international assessments. The code allows further extensions by different users for different applications and types, e.g. to other components of the Earth System. User modifications are encouraged and easy to perform with a minimum of coding.

Description: Description of a hybrid ice sheet-shelf model, and application to Antarctica Geoscientific Model Development Discussions, 5, 1077-1134, 2012 Author(s): D. Pollard and R. M. DeConto The formulation of a 3-D ice sheet-shelf model is described. The model is designed for long-term continental-scale applications, and has been used mostly in paleoclimatic studies. It uses a hybrid combination of the scaled Shallow Ice and Shallow Shelf Approximations for ice flow. Floating ice shelves and grounding-line migration are included, with parameterized ice fluxes at grounding lines that allows relatively coarse resolutions to be used. All significant components and parameterizations of the model are described in some detail. Basic results for modern Antarctica are compared with observations, and simulations over the last 5 million yr are shown to be similar to previously published results using an earlier model version. The sensitivity of ice retreat during the last deglaciation to basal sliding coefficients is discussed.

Description: A new marine ecosystem model for the University of Victoria Earth system climate model Geoscientific Model Development Discussions, 5, 1135-1201, 2012 Author(s): D. P. Keller, A. Oschlies, and M. Eby Earth system climate models (ESCMs) are valuable tools that can be used to gain a better understanding of the climate system, global biogeochemical cycles, and how anthropogenically-driven changes may affect them. Here we describe improvements made to the marine biogeochemical ecosystem component of the University of Victoria's ESCM (version 2.9). Major changes include corrections to the code and equations describing phytoplankton light limitation and zooplankton grazing, the implementation of a more realistic zooplankton growth and grazing model, and the implementation of an iron limitation scheme to constrain phytoplankton growth. The new model is evaluated after a 10 000-yr spin-up and compared to both the previous version and observations. For the majority of biogeochemical tracers and ecosystem processes the new model shows significant improvements when compared to the previous version and evaluated against observations. Many of the improvements are due to better simulation of seasonal changes in higher latitude ecosystems and the effect that this has on ocean biogeochemistry. This improved model is intended to provide a basic new ESCM model component, which can be used as is or expanded upon (i.e., the addition of new tracers), for climate change and biogeochemical cycling research.

Description: Development of a parameterization of black carbon aging for use in general circulation models Geoscientific Model Development Discussions, 5, 1263-1293, 2012 Author(s): N. Oshima and M. Koike A parameterization of black carbon (BC) aging in the atmosphere is developed for use in general circulation models (GCMs) that separately treats size distributions of hydrophobic BC and hydrophilic BC using lognormal modes. The rate of BC aging is expressed as the conversion rate from hydrophobic BC to hydrophilic BC modes (i.e., inverse of the e-folding time of the conversion, 1/ τ BC ). In this study, the conversion rates are estimated using results of detailed calculations by a size and mixing state resolved aerosol box model with numerous initial conditions. We introduce a new concept, the hydrophobic-BC-mass-normalized coating rate ( V BC ), defined as the rate of increase of the mass concentration of condensed materials on hydrophobic BC normalized by the hydrophobic BC mass concentration. Although the conversion rate largely varies depending on the atmospheric conditions and the concentrations of chemical species, we find that the variations of the conversion rate are generally expressed well by a unique function of V BC for given lognormal size distributions of hydrophobic BC. The parameterized conversion rate is expressed as a function of V BC , which enables the representation of diurnal and seasonal variations of the BC aging rate and its spatial differences in polluted and clean air, while other widely used constant conversion rates cannot. Application of our newly developed parameterization to GCMs will provide more reliable estimates of the spatial distribution of BC and its radiative effects at regional and global scales.

Description: Pre-industrial and mid-Pliocene simulations with NorESM-L – AGCM simulations Geoscientific Model Development Discussions, 5, 1203-1227, 2012 Author(s): Z. Zhang and Q. Yan In the Pliocene Model Intercomparison Project (PlioMIP), two sets of experiments are suggested. One includes a reference and a mid-Pliocene experiment run with atmosphere general circulation models (AGCM experiments, referred to as Experiments I), the other includes a pre-industrial and a mid-Pliocene experiment run with coupled ocean-atmosphere general circulation models (AOGCM experiments, referred to as Experiment II). In this paper, we describe the AGCM experiments with the atmosphere model in the low resolution version of the Norwegian Earth System Model (NorESM-L), and also assess the potential uncertainties in analyzing mid-Pliocene climate anomalies, due to choosing SST fields for the reference experiment. We carry out a mid-Pliocene experiment, a control experiment forced by the modern SST fields, and a pre-industrial experiment forced by the monthly SST fields from HadISST averaged between 1879 and 1900. Our experiments illustrate that the simulated mid-Pliocene global annual mean SAT is 17.1 °C. It is 2.5 °C warmer than the control experiment, but 2.7 °C warmer than the pre-industrial experiment. We find that the uncertainties in analyses of mid-Pliocene climate anomalies are small on a global scale, but still large on a regional scale. On the regional scale, these uncertainties should be noticed and assessed in future PlioMIP studies.

Description: Geomagnetic jerks characterization via spectral analysis Solid Earth, 3, 131-148, 2012 Author(s): B. Duka, A. De Santis, M. Mandea, A. Isac, and E. Qamili In this study we have applied spectral techniques to analyze geomagnetic field time-series provided by observatories, and compared the results with those obtained from analogous analyses of synthetic data estimated from models. Then, an algorithm is here proposed to detect the geomagnetic jerks in time-series, mainly occurring in the eastern component of the geomagnetic field. Applying such analysis to time-series generated from global models has allowed us to depict the most important space-time features of the geomagnetic jerks all over the globe, since the beginning of XXth century. Finally, the spherical harmonic power spectrum of the third derivative of the main geomagnetic field has been computed from 1960 to 2002.5, bringing new insights to understand the spatial evolution of these rapid changes of the geomagnetic field.

Description: Effect of glacial-interglacial sea-level changes on the displacement and stress field in the forearc and along the plate interface of subduction zones Solid Earth, 3, 63-70, 2012 Author(s): T. Li and A. Hampel Combined seismological, space-geodetic and numerical studies have shown that the seismicity at subduction zones may be modulated by tides and glacier fluctuations on timescales of 1–100 a, because these changes in loads on Earth's surface are able to alter the stress field in the upper plate and along the plate interface. Here we use a two-dimensional finite-element model of a subduction zone to investigate how glacial-interglacial sea-level changes affect the forearc region and the plate interface. The model results show that a sea-level fall by 125 m over 100 ka causes up to 0.7 m of vertical displacement, with the maximum uplift occurring between the trench and the coast. The uplift signal induced by the sea-level fall decreases to zero ~20 km landward of the coastline. A subsequent sea-level rise by 125 m over 20 ka causes subsidence, which is again most pronounced in the submarine part of the forearc. The sea-level changes cause horizontal displacements of up to 0.12 m, which are directed seaward during sea-level fall and landward during sea-level rise. With respect to the stress field, the sea-level changes lead to variations in the vertical stress and the shear stress of up to 1.23 MPa and 0.4 MPa, respectively. The shear stress variations are highest beneath the coast, i.e. in the area where the sea-level changes cause the strongest flexure. The resulting Coulomb stress changes on the plate interface are of the order of 0.2–0.5 MPa and indicate that earthquakes are promoted during sea-level fall and delayed during sea-level rise. Our findings imply that eustatic sea-level changes during glacial-interglacial periods may have induced displacements and stress changes that were large enough to affect the seismic cycle of subduction thrusts.

Description: Seasonal leaf dynamics for tropical evergreen forests in a process based global ecosystem model Geoscientific Model Development Discussions, 5, 639-681, 2012 Author(s): M. De Weirdt, H. Verbeeck, F. Maignan, P. Peylin, B. Poulter, D. Bonal, P. Ciais, and K. Steppe The influence of seasonal phenology in tropical humid forests on canopy photosynthesis remains poorly understood and its representation in global vegetation models highly simplified, typically with no seasonal variability of canopy leaf area properties taken into account. However, recent flux tower and remote sensing studies suggest that seasonal phenology in tropical rainforests exerts a large influence over carbon and water fluxes, with feedbacks that can significantly influence climate dynamics. A more realistic description of the underlying mechanisms that drive seasonal tropical forest photosynthesis and phenology could improve the correspondence of global vegetation model outputs with the wet-dry season biogeochemical patterns measured at flux tower sites. Here, we introduce a leaf Net Primary Production (NPP) based canopy dynamics scheme for evergreen tropical forests in the global terrestrial ecosystem model ORCHIDEE and validated the new scheme against in-situ carbon flux measurements. Modelled Gross Primary Productivity (GPP) patterns are analyzed in details for a flux tower site in French Guiana, in a forest where the dry season is short and where the vegetation is considered to have developed adaptive mechanisms against drought stress. By including leaf litterfall seasonality and a coincident light driven leaf flush and seasonal change in photosynthetic capacity in ORCHIDEE, modelled carbon and water fluxes more accurately represent the observations. The fit to GPP flux data was substantially improved and the results confirmed that by modifying canopy dynamics to benefit from increased light conditions, a better representation of the seasonal carbon flux patterns was made.

Description: The regional aerosol-climate model REMO-HAM Geoscientific Model Development Discussions, 5, 737-779, 2012 Author(s): J.-P. Pietikäinen, D. O'Donnell, C. Teichmann, U. Karstens, S. Pfeifer, J. Kazil, R. Podzun, S. Fiedler, H. Kokkola, W. Birmili, C. O'Dowd, U. Baltensperger, E. Weingartner, R. Gehrig, G. Spindler, M. Kulmala, J. Feichter, D. Jacob, and A. Laaksonen REMO-HAM is a new regional aerosol-climate model. It is based on the REMO regional climate model and includes all of the major aerosol processes. The structure for aerosol is similar to the global aerosol-climate model ECHAM5-HAM, for example the aerosol module HAM-M7 has been coupled with a two-moment stratiform cloud scheme. In this work, we have evaluated the model and compared the results against ECHAM5-HAM and measurements. Four different measurement sites was chosen for the comparison of total number concentrations, size distributions and gas phase sulfur dioxide concentrations: Hyytiälä in Finland, Melpitz in Germany, Mace Head in Ireland and Jungfraujoch in Switzerland. REMO-HAM is run with two different resolutions: 50×50 km 2 and 10×10 km 2 . Based on our simulations, REMO-HAM can represent the measured values reasonably well. The total number concentrations are slightly underestimated, which is probably due to the missing boundary layer nucleation and online secondary organic aerosol model. The differences in the total number concentrations between REMO-HAM and ECHAM5-HAM can be mainly explained by the difference in the nucleation mode. From the meteorological point of view, REMO-HAM represents the precipitation fields and 2 m temperature profile very well compared to measurement. Overall, we have shown that REMO-HAM is a functional aerosol-climate model, which will be used in further studies.

Description: DInSAR Coseismic Deformation of the May 2011 M w 5.1 Lorca Earthquake (southeastern Spain) Solid Earth, 3, 111-119, 2012 Author(s): T. Frontera, A. Concha, P. Blanco, A. Echeverria, X. Goula, R. Arbiol, G. Khazaradze, F. Pérez, and E. Suriñach The coseismic superficial deformation at the region of Lorca (Murcia, southeastern Spain) due to the M w 5.1 earthquake on 11 May 2011 was characterized by a multidisciplinary team, integrating information from DInSAR, GPS and numerical modelling techniques. Despite the moderate magnitude of the event, quantitative information was obtained from the interferometric study of a pair of TerraSAR-X images. The DinSAR results defined the trace of the fault plane and evidenced uplift of the hanging wall block in agreement with the estimated deformation obtained through an elastic rupture dislocation numerical model. Meanwhile for the footwall block, interferometric results showed that tectonic deformation is masked by an important subsidence related to groundwater extraction previously identified at the area of study. Horizontal crustal deformation rates and velocity vectors, obtained from GPS stations existent at the area, were also coherent with the tectonic setting of the southern margin of the Iberian Peninsula and with the focal mechanism calculated for the Lorca event. The analysis of a continuous GPS site in Lorca showed good agreement with the horizontal N–S direction component relative to the numerical model and tectonics of the region. This is the first time at this seismic active area that a multi-technique analysis has been performed immediately after the occurrence of a seismic event, comparing the existing deformation data with a theoretical numerical model based on estimated seismic rupture dislocation.

Description: Floating stones off El Hierro, Canary Islands: xenoliths of pre-island sedimentary origin in the early products of the October 2011 eruption Solid Earth, 3, 97-110, 2012 Author(s): V. R. Troll, A. Klügel, M.-A. Longpré, S. Burchardt, F. M. Deegan, J. C. Carracedo, S. Wiesmaier, U. Kueppers, B. Dahren, L. S. Blythe, T. H. Hansteen, C. Freda, D. A. Budd, E. M. Jolis, E. Jonsson, F. C. Meade, C. Harris, S. E. Berg, L. Mancini, M. Polacci, and K. Pedroza A submarine eruption started off the south coast of El Hierro, Canary Islands, on 10 October 2011 and continues at the time of this writing (February 2012). In the first days of the event, peculiar eruption products were found floating on the sea surface, drifting for long distances from the eruption site. These specimens, which have in the meantime been termed "restingolites" (after the close-by village of La Restinga), appeared as black volcanic "bombs" that exhibit cores of white and porous pumice-like material. Since their brief appearance, the nature and origin of these "floating stones" has been vigorously debated among researchers, with important implications for the interpretation of the hazard potential of the ongoing eruption. The "restingolites" have been proposed to be either (i) juvenile high-silica magma (e.g. rhyolite), (ii) remelted magmatic material (trachyte), (iii) altered volcanic rock, or (iv) reheated hyaloclastites or zeolite from the submarine slopes of El Hierro. Here, we provide evidence that supports yet a different conclusion. We have analysed the textures and compositions of representative "restingolites" and compared the results to previous work on similar rocks found in the Canary Islands. Based on their high-silica content, the lack of igneous trace element signatures, the presence of remnant quartz crystals, jasper fragments and carbonate as well as wollastonite (derived from thermal overprint of carbonate) and their relatively high oxygen isotope values, we conclude that "restingolites" are in fact xenoliths from pre-island sedimentary layers that were picked up and heated by the ascending magma, causing them to partially melt and vesiculate. As they are closely resembling pumice in appearance, but are xenolithic in origin, we refer to these rocks as "xeno-pumice". The El Hierro xeno-pumices hence represent messengers from depth that help us to understand the interaction between ascending magma and crustal lithologies beneath the Canary Islands as well as in similar Atlantic islands that rest on sediment-covered ocean crust (e.g. Cape Verdes, Azores). The occurrence of "restingolites" indicates that crustal recycling is a relevant process in ocean islands, too, but does not herald the arrival of potentially explosive high-silica magma in the active plumbing system beneath El Hierro.

Description: Carbon-nitrogen feedbacks in the UVic ESCM Geoscientific Model Development Discussions, 5, 67-118, 2012 Author(s): R. Wania, K. J. Meissner, M. Eby, V. Arora, I. Ross, and A. J. Weaver A representation of the terrestrial nitrogen cycle is introduced into the UVic Earth System Climate Model (UVic ESCM). The UVic ESCM now contains five terrestrial carbon pools and seven terrestrial nitrogen pools: soil, litter, leaves, stem and roots for both elements and ammonium and nitrate in the soil for nitrogen. Nitrogen cycles through plant tissue, litter, soil and the mineral pools before being taken up again by the plant. Biological N 2 fixation and nitrogen deposition represent the external input and loss from the plant-soil system can occur via leaching. Simulated carbon and nitrogen pools and fluxes are in the range of other models and data. Gross primary production (GPP) for the 1990s in the CN-coupled version is 129.6 Pg C a −1 and net C uptake is 0.83 Pg C a −1 , whereas the C-only version results in a GPP of 133.1 Pg C a −1 and a net C uptake of 1.57 Pg C a −1 . At the end of a transient experiment for the years 1800–2000, where temperature is held constant but CO 2 fertilisation for vegetation is allowed to happen, the CN-coupled version shows an enhanced net C uptake of 1.05 Pg C a −1 , whereas in the experiment where CO 2 is held constant and temperature is transient the land turns into a C source of 0.60 Pg C a −1 by the 1990s. The arithmetic sum of the temperature and CO 2 effects results in 0.45 Pg C a −1 , which is 0.38 Pg C a −1 lower than seen in the fully forced model, suggesting a strong non-linearity in the CN-coupled version. Anthropogenic N deposition has a positive effect on Net Ecosystem Production of 0.35 Pg C a −1 . Overall, the UVic CN-coupled version shows similar characteristics in terms of C and N pools and fluxes to other CN-coupled Earth System Models.

Description: Pre-industrial and mid-Pliocene simulations with NorESM-L Geoscientific Model Development Discussions, 5, 119-148, 2012 Author(s): Z. S. Zhang, K. Nisancioglu, M. Bentsen, J. Tjiputra, I. Bethke, Q. Yan, B. Risebrobakken, C. Andersson, and E. Jansen The mid-Pliocene period (3.3 to 3.0 Ma) is known as a warm climate with atmospheric greenhouse gas levels similar to the present. As the climate at this time was in equilibrium with the greenhouse forcing, it is a valuable test case to better understand the long term response to high levels of atmospheric greenhouse gases. In this study, we use the low resolution version of the Norwegian Earth System Model (NorESM-L) to simulate the pre-industrial and the mid-Pliocene climate. Comparison of the simulation with observations demonstrates that NorESM-L simulates a realistic pre-industrial climate. The simulated mid-Pliocene global mean surface air temperature is 16.7 °C, which is 3.2 °C warmer than the pre-industrial. The simulated mid-Pliocene global mean sea surface temperature is 19.1 °C, which is 2.0 °C warmer than the pre-industrial. The warming is relatively uniform globally, except for a strong amplification at high latitudes.

Description: Implementation and evaluation of online gas-phase chemistry within a regional climate model (RegCM-CHEM4) Geoscientific Model Development Discussions, 5, 149-188, 2012 Author(s): A. K. Shalaby, A. S. Zakey, A. B. Tawfik, F. Solmon, F. Giorgi, F. Stordal, S. Sillman, R. A. Zaveri, and A. L. Steiner The RegCM-CHEM4 is a new online climate-chemistry model based on the International Centre for Theoretical Physics (ICTP) regional climate model (RegCM4). Tropospheric gas-phase chemistry is integrated into the climate model using the condensed version of the Carbon Bond Mechanism (CBM-Z; Zaveri and Peters, 1999) with a fast solver based on radical balances. We evaluate the model over Continental Europe for two different time scales: (1) an event-based analysis of the ozone episode associated with the heat wave of August 2003 and (2) a climatological analysis of a six-year simulation (2000–2005). For the episode analysis, model simulations show good agreement with European Monitoring and Evaluation Program (EMEP) observations of hourly ozone over different regions in Europe and capture ozone concentrations during and after the August 2003 heat wave event. For long-term climate simulations, the model captures the seasonal cycle of ozone concentrations with some over prediction of ozone concentrations in non-heat wave summers. Overall, the ozone and ozone precursor evaluation shows the feasibility of using RegCM-CHEM4 for decadal-length simulations of chemistry-climate interactions.

Description: Bio-chemostratigraphy of the Barremian-Aptian shallow-water carbonates of the southern Apennines (Italy): pinpointing the OAE1a in a Tethyan carbonate platform Solid Earth, 3, 1-28, 2012 Author(s): M. Di Lucia, A. Trecalli, M. Mutti, and M. Parente Low biostratigraphic resolution and lack of chronostratigraphic calibration hinder precise correlations between platform carbonates and coeval deep-water successions. These are the main obstacle when studying the record of Mesozoic oceanic anoxic events in carbonate platforms. In this paper carbon and strontium isotope stratigraphy are used to produce the first chronostratigraphic calibration of the Barremian-Aptian biostratigraphy of the Apenninic carbonate platform of southern Italy. According to this calibration, the segment of decreasing δ 13 C values, leading to the negative peak that is generally taken as the onset of the Selli event, starts a few metres above the last occurrence of Palorbitolina lenticularis and Voloshinoides murgensis . The following rise of δ 13 C values, corresponding to the interval of enhanced accumulation of organic matter in deep-water sections, ends just below the first acme of Salpingoporella dinarica , which roughly corresponds to the segment of peak δ 13 C values. The whole carbon isotope excursion associated with the oceanic anoxic event 1a is bracketed in the Apenninic carbonate platform between the last occurrence of Voloshinoides murgensis and the "Orbitolina level", characterized by the association of Mesorbitolina parva and Mesorbitolina texana . Since these bioevents have been widely recognized beyond the Apenninic platform, the calibration presented in this paper can be used to pinpoint the interval corresponding to the Early Aptian oceanic anoxic event in other carbonate platforms of central and southern Tethys. This calibration will be particularly useful to interpret the record of the Selli event in carbonate platform sections for which a reliable carbon isotope stratigraphy is not available.

Description: Downscaling the climate change for oceans around Australia Geoscientific Model Development Discussions, 5, 425-458, 2012 Author(s): M. A. Chamberlain, C. Sun, R. J. Matear, M. Feng, and S. J. Phipps At present, global climate models used to project changes in climate do not resolve mesoscale ocean features such as boundary currents and eddies. These missing features may be important to realistically project the marine impacts of climate change. Here we present a framework for dynamically downscaling coarse climate change projections utilising a global ocean model that resolves these features in the Australian region. The downscaling model used here is ocean-only. The ocean feedback on the air-sea fluxes is explored by restoring to surface temperature and salinity, as well as a calculated feedback to wind stress. These feedback approximations do not replace the need for fully coupled models, but they allow us to assess the sensitivity of the ocean in downscaled climate change simulations. Significant differences are found in sea surface temperature, salinity, stratification and transport between the downscaled projections and those of the climate model. While the magnitude of the climate change differences may vary with the feedback parameterisation used, the patterns of the climate change differences are consistent and develop rapidly indicating they are mostly independent of feedback that ocean differences may have on the air-sea fluxes. Until such a time when it is feasible to regularly run a global climate model with eddy resolution, our framework for ocean climate change downscaling provides an attractive way to explore how climate change may affect the mesoscale ocean environment.

Description: Mid-Pliocene global climate simulation with MRI-CGCM2.3: set-up and initial results of PlioMIP Experiments 1 and 2 Geoscientific Model Development Discussions, 5, 383-423, 2012 Author(s): Y. Kamae and H. Ueda The mid-Pliocene (3.3 to 3.0 million yr ago), a globally warm period before the Quaternary, is recently attracting attention as a new target for paleoclimate modelling and data-model synthesis. This paper reports set-ups and results of experiments proposed in Pliocene Model Intercomparison Project (PlioMIP) using with a global climate model, MRI-CGCM2.3. We conducted pre-industrial and mid-Pliocene runs by using of the coupled atmosphere-ocean general circulation model (AOGCM) and its atmospheric component (AGCM) for the PlioMIP Experiments 2 and 1, respectively. In addition, we conducted two types of integrations in AOGCM simulation, with and without flux adjustments on sea surface. General characteristics of differences in the simulated mid-Pliocene climate relative to the pre-industrial in the three integrations are compared in this study. Generally, difference of simulated surface climate between AGCM and AOGCM is larger than that between the two AOGCM runs, with and without flux adjustments. The simulated climate shows different pattern between AGCM and AOGCM particularly over low latitude oceans, subtropical land regions, and high latitude oceans. The AOGCM simulations do not reproduce wetter environment in the subtropics relative to the present-day, which is suggested by terrestrial proxy data. The differences between the two types of AOGCM runs are little over the land but evident over the ocean particularly in the North Atlantic and polar regions.

Description: Phanerozoic environments of black shale deposition and the Wilson Cycle Solid Earth, 3, 29-42, 2012 Author(s): J. Trabucho-Alexandre, W. W. Hay, and P. L. de Boer The spatial and temporal distribution of black shales is related to the development of environments in which they accumulate and to a propitious combination of environmental variables. In recent years, much has been done to improve our understanding of the mechanisms behind the temporal distribution of black shales in the Phanerozoic and of the environmental variables that result in their deposition. However, the interpretation of ancient black shale depositional environments is dominated by an oversimplistic set of three depositional models that do not capture their complexity and dynamics. These three models, the restricted circulation, the (open) ocean oxygen minimum and the continental shelf models, are an oversimplification of the variety of black shale depositional environments that arise and coexist throughout the course of a basin's Wilson Cycle, i.e. the dynamic sequence of events and stages that characterise the evolution of an ocean basin, from the opening continental rift to the closing orogeny. We examine the spatial distribution of black shales in the context of the Wilson Cycle using examples from the Phanerozoic. It is shown that the geographical distribution of environments of black shale deposition and the position of black shales in the basin infill sequence strongly depend on basin evolution, which controls the development of sedimentary environments where black shales may be deposited. The nature of the black shales that are deposited, i.e. lithology and type of organic matter, also depends on basin evolution and palaeogeography. We propose that in studies of black shales more attention should be given to the sedimentary processes that have led to their formation and to the interpretation of their sedimentary environments.

Description: Land surface Verification Toolkit (LVT) – a generalized framework for land surface model evaluation Geoscientific Model Development Discussions, 5, 229-276, 2012 Author(s): S. V. Kumar, C. D. Peters-Lidard, J. Santanello, K. Harrison, Y. Liu, and M. Shaw Model evaluation and verification are key in improving the usage and applicability of simulation models for real-world applications. In this article, the development and capabilities of a formal system for land surface model evaluation called the Land surface Verification Toolkit (LVT) is described. LVT is designed to provide an integrated environment for systematic land model evaluation and facilitates a range of verification approaches and analysis capabilities. LVT operates across multiple temporal and spatial scales and employs a large suite of in-situ, remotely sensed and other model and reanalysis datasets in their native formats. In addition to the traditional accuracy-based measures, LVT also includes uncertainty and ensemble diagnostics, information theory measures, spatial similarity metrics and scale decomposition techniques that provide novel ways for performing diagnostic model evaluations. Though LVT was originally designed to support the land surface modeling and data assimilation framework known as the Land Information System (LIS), it supports hydrological data products from non-LIS environments as well. In addition, the analysis of diagnostics from various computational subsystems of LIS including data assimilation, optimization and uncertainty estimation are supported within LVT. Together, LIS and LVT provide a robust end-to-end environment for enabling the concepts of model data fusion for hydrological applications. The evolving capabilities of LVT framework are expected to facilitate rapid model evaluation efforts and aid the definition and refinement of formal evaluation procedures for the land surface modeling community.

Description: Water isotope variations in the global ocean model MPI-OM Geoscientific Model Development Discussions, 5, 277-307, 2012 Author(s): X. Xu, M. Werner, M. Butzin, and G. Lohmann The stable water isotopes H 2 18 O and HDO are incorporated as passive tracers into the oceanic general circulation model MPI-OM, and a control simulation under present-day climate conditions is analyzed in detail. Both δ 18 O and δD distributions at the ocean surface and deep ocean are generally consistent with available observations on the large scale. The modelled δD-δ 18 O relations in surface waters slightly deviates from the slope of the global meteoric water line in most basins, and a much steeper slope is detected in Arctic Oceans. The simulated deuterium excess of ocean surface waters shows small variations between 80° S and 55° N, and a strong decrease north of 55° N. The model is also able to capture the quasi-linear relationship between δ 18 O and salinity S , as well as δD and S , as seen in observational data. Both in the model results and observations, the surface δ– S relations show a steeper slope in extra-tropical regions than in tropical regions, which indicates relatively more addition of isotopically depleted water at high latitudes.

Description: Assimilation of OMI NO 2 retrievals into the limited-area chemical transport model DEHM (V2007.0) with a 2-D OI algorithm Geoscientific Model Development Discussions, 5, 309-346, 2012 Author(s): J. D. Silver, J. Brandt, M. Hvidberg, and J. Frydendall Data assimilation is the process of combining real-world observations with a modelled geophysical field. The increasing abundance of satellite retrievals of atmospheric trace gases makes chemical data assimilation a powerful tool for improving air quality forecasts. We implemented a two-dimensional optimal interpolation (OI) algorithm to assimilate satellite-derived estimates of tropospheric NO 2 column concentrations into the Danish Eulerian Hemispheric Model (DEHM, version V2007.0), a three-dimensional, European-scale, chemical transport model. In particular, we describe how we used observational data to estimate the background error covariance matrix, B . In the assimilation, the tropospheric column NO 2 field was adjusted and the modelled NO 2 profile was scaled accordingly; other species were only adjusted indirectly via changes to NO 2 concentrations. We ran a number of experiments to compare different parameterisations of B ; this involved varying the length scale used in B , the relative weighting of the background and observation errors, the errors assigned to observations and the influence of clustered observations. We assessed model performance by comparing the analysed fields to an independent set of observations: ground-based measurements of NO 2 concentrations. Ozonosonde profiles were also used for verification. The analysed NO 2 and O 3 concentrations were more accurate than those from a reference simulation without assimilation, with lower bias for both species and improved correlation for NO 2 . The experiments showed that appropriately chosen parameters for the B matrix, estimated using innovation statistics, yielded more accurate surface NO 2 concentrations. There was good agreement between the seasonally-averaged observed and modelled O 3 profiles. The simple OI scheme was effective and computationally feasible in this context, where only a single species was assimilated and only a two-dimensional field was adjusted. However there are certain limitations to using this assimilation scheme for more highly multi-dimensional problems. Although forecast accuracy was not examined here, we discuss the potential for improving NO 2 forecasts by using assimilation to generate initial conditions.

Description: Assessing climate model software quality: a defect density analysis of three models Geoscientific Model Development Discussions, 5, 347-382, 2012 Author(s): J. Pipitone and S. Easterbrook A climate model is an executable theory of the climate; the model encapsulates climatological theories in software so that they can be simulated and their implications investigated. Thus, in order to trust a climate model one must trust that the software it is built from is built correctly. Our study explores the nature of software quality in the context of climate modelling. We performed an analysis of defect reports and defect fixes in several versions of leading global climate models by collecting defect data from bug tracking systems and version control repository comments. We found that the climate models all have very low defect densities compared to well-known, similarly sized open-source projects. We discuss the implications of our findings for the assessment of climate model software trustworthiness.

Description: MAESPA: a model to study interactions between water limitation, environmental drivers and vegetation function at tree and stand levels, with an example application to [CO 2 ] × drought interactions Geoscientific Model Development Discussions, 5, 459-513, 2012 Author(s): R. A. Duursma and B.E. Medlyn Process-based models (PBMs) of vegetation function can be used to interpret and integrate experimental results. Water limitation to plant carbon uptake is a highly uncertain process in the context of environmental change, and many experiments have been carried out that study drought limitations to vegetation function at spatial scales from seedlings to entire canopies. What is lacking in the synthesis of these experiments is a quantitative tool that can be used to integrate and analyse experimental results at scales of both the whole-plant and the forest canopy, and that includes a detailed mechanistic representation of the water balance. To fill this gap, we developed an individual tree-based model (MAESPA), largely based on combining the well-known MAESTRA and SPA ecosystem models. The model includes a hydraulically-based model of stomatal conductance, root water uptake routines, drainage, infiltration, runoff and canopy interception, as well as detailed radiation interception and leaf physiology routines from the MAESTRA model. The model can be applied both to single plants of arbitrary size and shape, as well as stands of trees. The utility of this model is demonstrated by studying the interaction between elevated [CO 2 ] (e C a ) and drought. Based on theory, this interaction is generally expected to be positive, so that plants growing in e C a should be less susceptible to drought. Experimental results, however, are varied. We apply the model to a previously published experiment on droughted cherry, and show that changes in plant parameters due to long-term growth at e C a (acclimation) may strongly affect the outcome of C a × drought experiments. We discuss potential applications of MAESPA and some of the key uncertainties in process representation.

Description: Towards a public, standardized, diagnostic benchmarking system for land surface models Geoscientific Model Development Discussions, 5, 549-570, 2012 Author(s): G. Abramowitz We examine different conceptions of land surface model benchmarking and illustrate the importance of internationally standardized evaluation experiments that specify data sets, variables, metrics and model resolutions. We additionally show how essential the definition of a priori expectations of model performance can be, based on the complexity of a model and the amount of information being provided to it, and give an example of how these expectations might be quantified. Finally, we introduce the Protocol for the Analysis of Land Surface models (PALS), a free, online land surface model benchmarking application, and show how it is structured to meet both of these goals.

Description: A standard test case suite for two-dimensional linear transport on the sphere Geoscientific Model Development Discussions, 5, 189-228, 2012 Author(s): P. H. Lauritzen, W. C. Skamarock, M. J. Prather, and M. A. Taylor It is the purpose of this paper to propose a standard test case suite for two-dimensional transport schemes on the sphere intended to be used for model development and facilitating scheme intercomparison. The test cases are designed to assess important aspects of accuracy in geophysical fluid dynamics such as numerical order of convergence, "minimal" resolution, the ability of the transport scheme to preserve filaments, transport "rough" distributions, and to preserve pre-existing functional relations between species/tracers under challenging flow conditions. The experiments are designed to be easy to set up. They are specified in terms of two analytical wind fields (one non-divergent and one divergent) and four analytical initial conditions (varying from smooth to discontinuous). Both conventional error norms as well as novel mixing and filament preservation diagnostics are used that are easy to implement. The experiments pose different challenges for the range of transport approaches from Lagrangian to Eulerian. The mixing and filament preservation diagnostics do not require an analytical/reference solution which is in contrast to standard error norms where a "true" solution is needed. Results using the CSLAM (Conservative Semi-Lagrangian Multi-tracer) scheme on the cubed-sphere are presented for reference and illustrative purposes.

Description: Simulations over South Asia using the weather research and forecasting model with chemistry (WRF-Chem): chemistry evaluation and initial results Geoscientific Model Development Discussions, 5, 1-66, 2012 Author(s): R. Kumar, M. Naja, G. G. Pfister, M. C. Barth, C. Wiedinmyer, and G. P. Brasseur This study presents annual simulations of tropospheric ozone and related species made for the first time using the WRF-Chem model over South Asia for the year 2008. The model simulated ozone, CO, and NO x are evaluated against ground-based, balloon-borne and satellite-borne (TES, OMI and MOPITT) observations. The comparison of model results with surface ozone observations from seven sites and CO and NO x observations from three sites, indicate the model's ability in reproducing seasonal variations of ozone and CO, but show some differences in NO x . The modeled vertical ozone distribution agrees well with the ozone soundings data from two Indian sites. The vertical distributions of TES ozone and MOPITT CO are generally well reproduced, but the model underestimates TES ozone, OMI tropospheric column NO 2 and MOPITT total column CO retrievals during all the months except MOPITT retrievals during August–January. Largest differences between modeled and satellite retrieved quantities are found during spring when intense biomass burning activity occurs in this region. The evaluation results indicate large uncertainties in anthropogenic and biomass burning emission estimates, especially for NO x . The model results indicate clear regional differences in the seasonality of surface ozone over South Asia with estimated net ozone production during daytime (11:30–15:30 h) over inland regions of 0–5 ppbv h −1 during all seasons and of 0–2 ppbv h −1 over marine regions during outflow periods. The model results indicate that ozone production in this region is mostly NO x -limited. This study shows that WRF-Chem model captures many important features of the observations and gives confidence to using the model for understanding the spatio-temporal variability of ozone over South Asia. However, improvements of South Asian emission inventories and simulations at finer model resolution, especially over the complex Himalayan terrain in Northern India, are also essential for accurately simulating ozone in this region.

Description: Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area Solid Earth, 3, 43-51, 2012 Author(s): P. Tong, D. Zhao, and D. Yang High-resolution tomographic images of the crust and upper mantle in and around the area of the 2011 Iwaki earthquake (M 7.0) and the Fukushima nuclear power plant are determined by inverting a large number of high-quality arrival times with both the finite-frequency and ray tomography methods. The Iwaki earthquake and its aftershocks mainly occurred in a boundary zone with strong variations in seismic velocity and Poisson's ratio. Prominent low-velocity and high Poisson's ratio zones are revealed under the Iwaki source area and the Fukushima nuclear power plant, which may reflect fluids released from the dehydration of the subducting Pacific slab under Northeast Japan. The 2011 Tohoku-oki earthquake (Mw 9.0) caused static stress transfer in the overriding Okhotsk plate, resulting in the seismicity in the Iwaki source area that significantly increased immediately following the Tohoku-oki mainshock. Our results suggest that the Iwaki earthquake was triggered by the ascending fluids from the Pacific slab dehydration and the stress variation induced by the Tohoku-oki mainshock. The similar structures under the Iwaki source area and the Fukushima nuclear power plant suggest that the security of the nuclear power plant site should be strengthened to withstand potential large earthquakes in the future.

Description: Modelling the mid-Pliocene Warm Period climate with the IPSL coupled model and its atmospheric component LMDZ4 Geoscientific Model Development Discussions, 5, 515-548, 2012 Author(s): C. Contoux, G. Ramstein, and A. Jost This paper describes the experimental design and model results of the climate simulations of the mid-Pliocene Warm Period (mPWP, ca. 3.3–3 Ma) using the Institut Pierre Simon Laplace model (IPSLCM5A), in the framework of the Pliocene Model Intercomparison Project (PlioMIP). We use the IPSL atmosphere ocean general circulation model (AOGCM), and its atmospheric component alone, to simulate the climate of the mPWP. Boundary conditions such as sea surface temperatures (SSTs), topography, ice sheet extent and vegetation are derived from the ones imposed by the Pliocene Model Intercomparison Project (PlioMIP), described in Haywood et al. (2010, 2011). We first describe the IPSL model main features, and then give a full description of the boundary conditions used for atmospheric model and coupled model experiments. The climatic outputs of the mPWP simulations are detailed and compared to the corresponding control simulations. The simulated warming is 1.94 °C in the atmospheric and 1.83 °C in the coupled model experiments. In both experiments, warming is more important at high latitudes. Simulated precipitation has a different behaviour in the coupled model than in the atmospheric model alone, because of the reduced gradients in imposed SSTs, which impacts the Hadley and Walker circulations. In addition, a sensitivity test to the change of land-sea mask in the atmospheric model, representing a sea-level change from present-day to 25 m higher during the mid-Pliocene, is described. We find that surface temperature differences can be important (several degrees Celsius) but are restricted to the areas that were changed from ocean to land or vice versa. In terms of precipitation, there is no impact on polar regions although the change in land-sea mask is important in these areas.

Description: Using open sidewalls for modelling self-consistent lithosphere subduction dynamics Solid Earth, 3, 313-326, 2012 Author(s): M. V. Chertova, T. Geenen, A. van den Berg, and W. Spakman Subduction modelling in regional model domains, in 2-D or 3-D, is commonly performed using closed (impermeable) vertical boundaries. Here we investigate the merits of using open boundaries for 2-D modelling of lithosphere subduction. Our experiments are focused on using open and closed (free slip) sidewalls while comparing results for two model aspect ratios of 3:1 and 6:1. Slab buoyancy driven subduction with open boundaries and free plates immediately develops into strong rollback with high trench retreat velocities and predominantly laminar asthenospheric flow. In contrast, free-slip sidewalls prove highly restrictive on subduction rollback evolution, unless the lithosphere plates are allowed to move away from the sidewalls. This initiates return flows pushing both plates toward the subduction zone speeding up subduction. Increasing the aspect ratio to 6:1 does not change the overall flow pattern when using open sidewalls but only the flow magnitude. In contrast, for free-slip boundaries, the slab evolution does change with respect to the 3:1 aspect ratio model and slab evolution does not resemble the evolution obtained with open boundaries using 6:1 aspect ratio. For models with open side boundaries, we could develop a flow-speed scaling based on energy dissipation arguments to convert between flow fields of different model aspect ratios. We have also investigated incorporating the effect of far-field generated lithosphere stress in our open boundary models. By applying realistic normal stress conditions to the strong part of the overriding plate at the sidewalls, we can transfer intraplate stress to influence subduction dynamics varying from slab roll-back, stationary subduction, to advancing subduction. The relative independence of the flow field on model aspect ratio allows for a smaller modelling domain. Open boundaries allow for subduction to evolve freely and avoid the adverse effects (e.g. forced return flows) of free-slip boundaries. We conclude that open boundaries in combination with intraplate stress conditions are to be preferred for modelling subduction evolution (rollback, stationary or advancing) using regional model domains.

Description: Stable water isotopes in the coupled atmosphere–land surface model ECHAM5-JSBACH Geoscientific Model Development Discussions, 5, 3375-3418, 2012 Author(s): B. Haese, M. Werner, and G. Lohmann In this study we present first results of a new model development, ECHAM5-JSBACH-wiso, where we have incorporated the stable water isotopes H 2 18 O and HDO as tracers in the hydrological cycle of the coupled atmosphere–land surface model ECHAM5-JSBACH. The ECHAM5-JSBACH-wiso model was run under present-day climate conditions at two different resolutions (T31L19, T63L31). A comparison between ECHAM5-JSBACH-wiso and ECHAM5-wiso shows that the coupling has a strong impact on the simulated temperature and soil wetness. Caused by these changes of temperature and the hydrological cycle, the δ 18 O in precipitation also shows variations from −4.5‰ up to 4.5‰. One of the clearest anomalies is shown over North-East Asia where, depending on an increase of temperature, the δ 18 O in precipitation increases as well. In order to analyze the sensitivity of the fractionation processes over land, we compare a set of simulations with various implementations of water isotope fractionation processes over the land surface. The simulations allow us to distinguish between no fractionation, fractionation included in the evaporation flux (from bare soil) and also fractionation included in both evaporation and transpiration (from water transport through plants) fluxes. The simulated δ 18 O and δD in precipitation of these setups generally fit well with the observations and the best agreement between observation and simulation is given in the case where no fractionation over land surface is assumed.

Description: Evaluation of roadway Gaussian plume models with large-scale measurement campaigns Geoscientific Model Development Discussions, 5, 3343-3373, 2012 Author(s): R. Briant, C. Seigneur, M. Gadrat, and C. Bugajny Gaussian models are commonly used to simulate atmospheric pollutant dispersion near sources because they provide an efficient compromise between reasonable accuracy and manageable computational time. The Gaussian dispersion formula provides an exact solution to the atmospheric diffusion equation for the dispersion of a pollutant emitted from a point source. However, the Gaussian dispersion formula for a line source, which is convenient to model emissions from on-road traffic, is exact only when the wind is perpendicular to the line source. A novel approach that reduces the error in the line source formula when the wind direction is not perpendicular to the road was recently developed. This model, combined with a Romberg integration to account for the road width, is used to simulate NO x concentrations in a large case study (1371 road sections representing about 831 km). NO 2 , NO and O 3 concentrations are then computed using the photostationary-state approximation. NO 2 concentrations are compared with measurements made at 242 locations in the domain area. Model performance is satisfactory with errors ranging from 24% to 31%. Results obtained here are also compared with those obtained with a previous formulation and with a standard model used for regulatory applications, ADMS-Urban. Discrepancies among the results obtained with those models are discussed.

Description: The SOCOL version 3.0 chemistry-climate model: description, evaluation, and implications from an advanced transport algorithm Geoscientific Model Development Discussions, 5, 3419-3467, 2012 Author(s): A. Stenke, M. Schraner, E. Rozanov, T. Egorova, B. Luo, and T. Peter We present the third generation of the coupled chemistry-climate model (CCM) SOCOL (modeling tools for studies of SOlar Climate Ozone Links). The most notable modifications compared to the previous model version are: (1) the dynamical core has been updated with the fifth generation of the middle-atmosphere general circulation model MA-ECHAM, and (2) the advection of the chemical species is now calculated by a mass-conserving and shape-preserving flux-form transport scheme instead of the previously used hybrid advection scheme. The whole chemistry code has been rewritten according to the ECHAM5 infrastructure and transferred to Fortran95. In contrast to its predecessors, SOCOLvs3 is now fully parallelized. The performance of the new SOCOL version is evaluated on the basis of transient model simulations (1975–2004) with different horizontal (T31 and T42) resolutions, following the approach of the CCMVal-1 model validation activity. The advanced advection scheme significantly reduces the artificial loss and accumulation of tracer mass in regions with strong gradients that was observed in previous model versions. Compared to its predecessors, SOCOLvs3 generally shows more realistic distributions of chemical trace species, especially of total inorganic chlorine, in terms of the mean state, but also of the annual and interannual variability. Advancements with respect to model dynamics are for example a better representation of the stratospheric mean state in spring, especially in the Southern Hemisphere, and a slowdown of the upward propagation in the tropical lower stratosphere. Despite a large number of improvements model deficiencies still remain. Examples include a too fast vertical ascent and/or horizontal mixing in the tropical stratosphere, the cold temperature bias in the lowermost polar stratosphere, and the overestimation of polar total ozone loss during Antarctic springtime.