An Overview of BCC Climate System Model Development and Application for Climate Change Studies

This paper reviews recent progress in the development of the Beijing Climate Center Climate System Model(BCC-CSM) and its four component models(atmosphere,land surface,ocean,and sea ice).Two recent versions are described:BCC-CSM1.1 with coarse resolution(approximately 2.8125°×2.8125°) and BCC-CSM1.1(m) with moderate resolution(approximately 1.125°×1.125°).Both versions are fully coupled climate-carbon cycle models that simulate the global terrestrial and oceanic carbon cycles and include dynamic vegetation.Both models well simulate the concentration and temporal evolution of atmospheric CO_2 during the 20th century with anthropogenic CO2 emissions prescribed.Simulations using these two versions of the BCC-CSM model have been contributed to the Coupled Model Intercomparison Project phase five(CMIP5) in support of the Intergovernmental Panel on Climate Change(IPCC) Fifth Assessment Report(AR5).These simulations are available for use by both national and international communities for investigating global climate change and for future climate projections.Simulations of the 20th century climate using BCC-CSMl.l and BCC-CSMl.l(m) are presented and validated,with particular focus on the spatial pattern and seasonal evolution of precipitation and surface air temperature on global and continental scales.Simulations of climate during the last millennium and projections of climate change during the next century are also presented and discussed.Both BCC-CSMl.l and BCC-CSMl.l(m) perform well when compared with other CMIP5 models.Preliminary analyses indicate that the higher resolution in BCC-CSM1.1(m) improves the simulation of mean climate relative to BCC-CSMl.l,particularly on regional scales.     

Ocean heat transport into the Arctic in the twentieth and twenty-first century in EC-Earth

The ocean heat transport into the Arctic and the heat budget of the Barents Sea are analyzed in an ensemble of historical and future climate simulations performed with the global coupled climate model EC-Earth. The zonally integrated northward heat flux in the ocean at 70°N is strongly enhanced and compensates for a reduction of its atmospheric counterpart in the twenty first century. Although an increase in the northward heat transport occurs through all of Fram Strait, Canadian Archipelago, Bering Strait and Barents Sea Opening, it is the latter which dominates the increase in ocean heat transport into the Arctic. Increased temperature of the northward transported Atlantic water masses are the main reason for the enhancement of the ocean heat transport. The natural variability in the heat transport into the Barents Sea is caused to the same extent by variations in temperature and volume transport. Large ocean heat transports lead to reduced ice and higher atmospheric temperature in the Barents Sea area and are related to the positive phase of the North Atlantic Oscillation. The net ocean heat transport into the Barents Sea grows until about year 2050. Thereafter, both heat and volume fluxes out of the Barents Sea through the section between Franz Josef Land and Novaya Zemlya are strongly enhanced and compensate for all further increase in the inflow through the Barents Sea Opening. Most of the heat transported by the ocean into the Barents Sea is passed to the atmosphere and contributes to warming of the atmosphere and Arctic temperature amplification. Latent and sensible heat fluxes are enhanced. Net surface long-wave and solar radiation are enhanced upward and downward, respectively and are almost compensating each other. We find that the changes in the surface heat fluxes are mainly caused by the vanishing sea ice in the twenty first century. The increasing ocean heat transport leads to enhanced bottom ice melt and to an extension of the area with bottom ice melt further northward. However, no indication for a substantial impact of the increased heat transport on ice melt in the Central Arctic is found. Most of the heat that is not passed to the atmosphere in the Barents Sea is stored in the Arctic intermediate layer of Atlantic water, which is increasingly pronounced in the twenty first century.     

Impact of climate change on the hydrogeology of two basins in northern France

This study presents an analysis of climate-change impacts on the water resources of two basins located in northern France, by integrating four sources of uncertainty: climate modelling, hydrological modelling, downscaling methods, and emission scenarios. The analysis focused on the evolution of the water budget, the river discharges and piezometric heads. Seven hydrological models were used, from lumped rainfall-discharge to distributed hydrogeological models, and led to quite different estimates of the water-balance components. One of the hydrological models, CLSM, was found to be unable to simulate the increased water stress and was, thus, considered as an outlier even though it gave fair results for the present day compared to observations. Although there were large differences in the results between the models, there was a marked tendency towards a decrease of the water resource in the rivers and aquifers (on average in 2050 about ?14 % and ?2.5 m, respectively), associated with global warming and a reduction in annual precipitation (on average in 2050 +2.1 K and ?3 %, respectively). The uncertainty associated to climate models was shown to clearly dominate, while the three others were about the same order of magnitude and 3–4 times lower. In terms of impact, the results found in this work are rather different from those obtained in a previous study, even though two of the hydrological models and one of the climate models were used in both studies. This emphasizes the need for a survey of the climatic-change impact on the water resource.     

Projection and possible causes of summer precipitation in eastern China using self-organizing map

Self-organizing map (SOM) is used to simulate summer daily precipitation over the Yangtze–Huaihe river basin in Eastern China, including future projections. SOM shows good behaviors in terms of probability distribution of daily rainfall and spatial distribution of rainfall indices, as well as consistency of multi-model simulations. Under RCP4.5 Scenario, daily rainfall at most sites (63%) is projected to shift towards larger values. For the early 21st century (2016–2035), precipitation in the central basin increases, yet decreases occur over the middle reaches of the Yangtze River as well as a part of its southeast area. For the late 21st century (2081–2100), the mean precipitation and extreme indices experience an overall increase except for a few southeast stations. The total precipitation in the lower reaches of the Yangtze River and in its south area is projected to increase from 7% at 1.5 °C global warming to 11% at 2 °C, while the intensity enhancement is more significant in southern and western sites of the domain. A clustering allows to regroup all SOM nodes into four distinct regimes. Such regional synoptic regimes show remarkable stability for future climate. The overall intensification of precipitation in future climate is linked to the occurrence-frequency rise of a wet regime which brings longitudinally closer the South Asia High (eastward extended) and the Western Pacific Subtropical High (westward extended), as well as the reduction of a dry pattern which makes the two atmospheric centers of action move away from each other.     

The stratospheric version of LMDz: dynamical climatologies, arctic oscillation, and impact on the surface climate

A climatology of the stratosphere is determined from a 20-year integration with the stratospheric version of the Atmospheric General Circulation Model LMDz. The model has an upper boundary at near 65 km, uses a Doppler spread non-orographic gravity waves drag parameterization and a subgrid-scale orography parameterization. It also has a Rayleigh damping layer for resolved waves only (not the zonal mean flow) over the top 5 km. This paper describes the basic features of the model and some aspects of its radiative-dynamical climatology. Standard first order diagnostics are presented but some emphasis is given to the model’s ability to reproduce the low frequency variability of the stratosphere in the winter northern hemisphere. In this model, the stratospheric variability is dominated at each altitudes by patterns which have some similarities with the arctic oscillation (AO). For those patterns, the signal sometimes descends from the stratosphere to the troposphere. In an experiment where the parameterized orographic gravity waves that reach the stratosphere are exaggerated, the model stratosphere in the NH presents much less variability. Although the stratospheric variability is still dominated by patterns that resemble to the AO, the downward influence of the stratosphere along these patterns is near entirely lost. In the same time, the persistence of the surface AO decreases, which is consistent with the picture that this persistence is linked to the descent of the AO signal from the stratosphere to the troposphere. A comparison between the stratospheric version of the model, and its routinely used tropospheric version is also done. It shows that the introduction of the stratosphere in a model that already has a realistic AO persistence can lead to overestimate the actual influence of the stratospheric dynamics onto the surface AO. Although this result is certainly model dependent, it suggests that the introduction of the stratosphere in a GCM also call for a new adjustment of the model parameters that affect the tropospheric variability.     

The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Crack nucleation in saturated porous media

Crack nucleation has been the subject of important contributions in the last two last decades. However, it seems that few attention has been granted to the case of saturated porous media. This is the question addressed in the present paper which is devoted to nucleation in traction mode. From a physical point of view, nucleation is a sudden phenomenon, so that the material response is both adiabatic and undrained. In the spirit of the variational approach, the nucleated crack is viewed as the final state of a region of space in which the material undergoes a full damage process. In traction mode, the opening of a saturated crack in undrained condition induces a drop of fluid pressure. In case of low fluid compressibility, the presence of the fluid delays the brittle failure usually associated with nucleation, as long as the fluid pressure remains above the saturation vapor pressure. Nucleation is therefore possible only if a partial vaporization of the fluid takes place.     

Magnitude of latent heat in thermally loaded clays

Temperature changes are known to induce specific couplings in clay, in particular, an anomalously high thermal pressurization in undrained conditions or a thermal compaction in drained conditions, both of which are potential threats for the mechanical stability and sealing capacity of the geomaterials. Thermodynamical analysis of those peculiar thermomechanical couplings points to a potentially important latent energy, which in turn could limit the temperature change upon heating or cooling. The direct measurement of latent energy developed during a laboratory geomechanical test is challenging. Instead, proper identification of thermal hardening in conventional experiments with temperature changes provides an alternative route to estimate latent energy. In this work, existing laboratory thermomechanical tests of clays are analyzed with a rigorous thermodynamic framework to quantify the magnitude of latent energy in thermomechanically loaded clays. A thermodynamically consistent constitutive model for fully saturated clays that combines two key features, (a) the temperature dependence of the blocked energy and (b) the framework of bounding plasticity, is proposed. The performance of the model is validated by reproducing results obtained in laboratory tests for Boom and Opalinus clays. The thermomechanical loads considered to validate the model performance were then used to estimate the percentage of work that remains latent in the clayey material during plastic yielding. We find that the magnitude of latent energy is quite significant, typically a few tens of percent of the total dissipated energy, and increases significantly with temperature. Accordingly, it is expected to play an important role in the thermomechanical response of clays.     

Provenance discrimination of fine sediments by mid-infrared spectroscopy: Calibration and application to fluvial palaeo-environmental reconstruction

Determining sediment provenance allows a better understanding of fluvial palaeo-dynamics, and identifying involved watersheds, at broad spatio-temporal scales. Conventional approaches for source identification are usually based on the physical, mineralogical, geochemical, magnetic or isotopic properties of sediments. Rapid, non-destructive and, in well-established contexts, highly accurate, mid-infrared spectroscopy is an alternative method for investigating sediment sources. The present research objectives are: (i) to use the mid-infrared spectroscopy method to discriminate the provenance of fine sediments, by applying discriminant analysis on a large set of reference samples from three different watersheds in the Upper Rhine area (associated with the Rhine, Ill and Vosges tributaries); (ii) to clarify whether the provenance spectra signatures are influenced by riverine depositional contexts (bars versus banks) and, to some extent, by grain size and/or high organic matter content; and (iii) to apply the mid-infrared spectroscopy – discriminant analysis method to a study of fluvial palaeo-dynamics and determine the provenance of palaeo-channel infillings. Three main sedimentary sources, divided into eight sub-categories, have been characterized by 196 modern reference samples from 78 collecting sites. Discriminant analysis displayed a strong separating power by classifying correctly the origin of samples without any inter-group overlap, independently from the geomorphological context (bar or bank) and associated slight changes in organic matter contents or grain size. Mid-infrared spectroscopy – discriminant analysis investigations of the palaeo-channel infill, complemented by radiocarbon dates and mineralogical data, allowed reconstructing general trends for the local morpho-sedimentary dynamics over the last ca 12 millennia.