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.     

A Simple And General Subgrid Model Suitable Both For Surface Layer And Free-Stream Turbulence

A new and general approach is presented to allow standard subgrid schemes to besuitable both for surface layer and free-stream turbulence. Simple modificationsto subgrid schemes are proposed and derived for any vertical stabilityconditions. They are simple to implement in models and are suitable for morecomplicated simulations such as large-eddy simulation with inhomogeneoussurface conditions or complex topography. They are also applicable to mesoscaleand large-scale models. These modifications are physically justified by recentmeasurements of spectra close to the ground. The spectral analysis presentedshows how the energy deficit of blocked turbulence for a given dissipation(`anomalous dissipation') dramatically affects the coefficients to be used insubgrid schemes. As shown for neutral and convective cases with wind shear,these changes allow us to substantially improve the prediction of profiles for themean quantities in the surface layer. Agreement with similarity laws in the unstablecase is found up to about 0.2z_i, for simulated shear, stabilityprofiles and dissipation rates of turbulent kinetic energy.     

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.     

Dissolution of calcitic marble and dolomitic rock in high iron concentrated acid mine drainage: application to anoxic limestone drains

Oxidation of sulphide mining waste can generate acid mine drainage (AMD) that has the potential to seriously affect the ecosystems. Acid mine drainage is characterised by a high acidity, high concentrations of sulphates and metals. To reduce the environmental impacts due to AMD, neutralisation using limestone drains is an option proposed in the literature and used around the world. The present study focuses on the influence of the carbonate rock mineralogy and their particle size on the neutralising capacity. The tests were performed in two different anoxic conditions: in batch reactors, and in columns having a hydraulic retention time of 15?h. The results showed that the neutralisation capacity of calcite was more important than for dolomitic rock, and smaller particle size gave higher alkalinity production (fine calcite dissolved faster in contact with AMD). A characterization of metal precipitate in sludge and in limestone coating was performed and demonstrated that gypsum, lepidocrocite and goethite were the predominant secondary minerals to be formed. Finally, this study underlines that anoxic limestone drain cannot be used alone to treat high iron concentrated AMD.     

The role of temperature on treeline migration for an eastern African mountain during the Last Glacial Maximum

Paleo-data suggest that East African mountain treelines underwent an altitudinal shift during the Last Glacial Maximum (LGM). Understanding the ecological and physiological processes underlying treeline response to such past climate change will help to improve forecasts of treeline change under future global warming. In spite of significant improvements in paleoclimatic reconstruction, the climatic conditions explaining this migration are still debated and important factors such as atmospheric CO₂ concentration, the impact of lapse rate decreasing temperature along altitudinal gradients and rainfall modifications due to elevation have often been neglected or simplified. Here, we assess the effects of these different factors and estimate the influence of the most dominant factors controlling changes in past treeline position using a multi-proxy approach based on simulations from BIOME4, a coupled biogeography and biogeochemistry model, modified to account for the effect of elevation on vegetation, compared with pollen, and isotopic data. The results indicate a shift in mountain vegetation at the LGM was controlled by low pCO₂ and low temperatures promoting species morphologically and physiologically better adapted to LGM conditions than many trees composing the forest belt limit. Our estimate that the LGM climate was cooler than today’s by ?4.5 °C (range: ?4.3 to ?4.6 °C) at the upper limit of the treeline, whereas at 831 m it was cooler by ?1.4 °C (range: ?2.6 to ?0.6 °C), suggests that a possible lapse rate modification strongly constrained the upper limit of treeline, which may limit its potential extension under future global warming.     

Mean flow-transient perturbation interaction in the Southern Hemisphere: a simulation using a variable-resolution GCM

The ability of an atmospheric general circulation model to reproduce fundamental features of the wintertime extratropical Southern Hemisphere (SH) circulation is evaluated with emphasis on the daily variability of the SH mean flow and the mean flow-transient perturbations interaction. Two 10-year simulations using a new version of the LMDZ GCM with a stretched grid scheme centered at 45 °S and forced by climatological SST are performed: a high (144Ꮡ) and low (64Ꭹ) horizontal resolution runs. The performance of both simulations was determined by comparing several simulated fields (zonal wind, temperature, kinetic energy, transient eddy momentum and heat fluxes, Eliassen-Palm fluxes, Eady growth rate and baroclinic conversion term) against the European Centre for Medium Range Weather Forecast reanalyses (ERA). High and low-resolution simulations are similar in many respects; in particular, both experiments reproduce the main patterns of the southern extratropical large-scale circulation satisfactorily. Increasing resolution does not improve universally some spurious aspects of the low resolution simulation (e.g. the cold bias in the high polar troposphere, the debilitated subtropical jet, the low baroclinic conversion rate). Those aspects present little sensitivity to the model resolution. The interaction between transient eddies and zonal mean flow are examined. The low-resolution experiment is able to qualitatively represent the acceleration/deceleration of the mean flow by transient perturbations, south/north of 30 °S with an accuracy similar to that of the high-resolution experiment. Although both experiments represent the baroclinic structure of the mean flow satisfactorily, the model underestimates some transient properties due to the underestimation of the baroclinic conversion term in middle latitudes. Such misrepresentation does not improve with increasing resolution and is related to the relatively weak meridional temperature gradient and the inadequate geographical distribution of the eddy heat fluxes. In particular, the eddy kinetic energy is always underestimated. Eddy kinetic energy does not improve convincingly with increasing resolution, suggesting that the adequate representation of the storm tracks is highly influenced by the physical parametrizations.     

Impacts of Indian Ocean SST biases on the Indian Monsoon: as simulated in a global coupled model

In this study, the impact of the ocean–atmosphere coupling on the atmospheric mean state over the Indian Ocean and the Indian Summer Monsoon (ISM) is examined in the framework of the SINTEX-F2 coupled model through forced and coupled control simulations and several sensitivity coupled experiments. During boreal winter and spring, most of the Indian Ocean biases are common in forced and coupled simulations, suggesting that the errors originate from the atmospheric model, especially a dry islands bias in the Maritime Continent. During boreal summer, the air-sea coupling decreases the ISM rainfall over South India and the monsoon strength to realistic amplitude, but at the expense of important degradations of the rainfall and Sea Surface Temperature (SST) mean states in the Indian Ocean. Strong SST biases of opposite sign are observed over the western (WIO) and eastern (EIO) tropical Indian Ocean. Rainfall amounts over the ocean (land) are systematically higher (lower) in the northern hemisphere and the south equatorial Indian Ocean rainfall band is missing in the control coupled simulation. During boreal fall, positive dipole-like errors emerge in the mean state of the coupled model, with warm and wet (cold and dry) biases in the WIO (EIO), suggesting again a significant impact of the SST errors. The exact contributions and the distinct roles of these SST errors in the seasonal mean atmospheric state of the coupled model have been further assessed with two sensitivity coupled experiments, in which the SST biases are replaced by observed climatology either in the WIO (warm bias) or EIO (cold bias). The correction of the WIO warm bias leads to a global decrease of rainfall in the monsoon region, which confirms that the WIO is an important source of moisture for the ISM. On the other hand, the correction of the EIO cold bias leads to a global improvement of precipitation and circulation mean state during summer and fall. Nevertheless, all these improvements due to SST corrections seem drastically limited by the atmosphere intrinsic biases, including prominently the unimodal oceanic position of the ITCZ (Inter Tropical Convergence Zone) during summer and the enhanced westward wind stress along the equator during fall.     

North Atlantic atmospheric circulation and surface wind in the Northeast of the Iberian Peninsula: uncertainty and long term downscaled variability

The variability and predictability of the surface wind field at the regional scale is explored over a complex terrain region in the northeastern Iberian Peninsula by means of a downscaling technique based on Canonical Correlation Analysis. More than a decade of observations (1992–2005) allows for calibrating and validating a statistical method that elicits the main associations between the large scale atmospheric circulation over the North Atlantic and Mediterranean areas and the regional wind field. In an initial step the downscaling model is designed by selecting parameter values from practise. To a large extent, the variability of the wind at monthly timescales is found to be governed by the large scale circulation modulated by the particular orographic features of the area. The sensitivity of the downscaling methodology to the selection of the model parameter values is explored, in a second step, by performing a systematic sampling of the parameters space, avoiding a heuristic selection. This provides a metric for the uncertainty associated with the various possible model configurations. The uncertainties associated with the model configuration are considerably dependent on the spatial variability of the wind. While the sampling of the parameters space in the model set up moderately impact estimations during the calibration period, the regional wind variability is very sensitive to the parameters selection at longer timescales. This fact illustrates that downscaling exercises based on a single configuration of parameters should be interpreted with extreme caution. The downscaling model is used to extend the estimations several centuries to the past using long datasets of sea level pressure, thereby illustrating the large temporal variability of the regional wind field from interannual to multicentennial timescales. The analysis does not evidence long term trends throughout the twentieth century, however anomalous episodes of high/low wind speeds are identified.     

Future regional projections of extreme temperatures in Europe: a nonstationary seasonal approach

This paper analyzes changes of maximum temperatures in Europe, which are evaluated using two state-of-the-art regional climate models from the EU ENSEMBLES project. Extremes are expressed in terms of return values using a time-dependent generalized extreme value (GEV) model fitted to monthly maxima. Unlike the standard GEV method, this approach allows analyzing return periods at different time scales (monthly, seasonal, annual, etc). The study focuses on the end of the 20th century (1961?C2000), used as a calibration/validation period, and assesses the changes projected for the period 2061?C2100 considering the A1B emission scenario. The performance of the regional models is evaluated for each season of the calibration period against the high-resolution gridded E-OBS dataset, showing a similar South-North gradient with larger values over the Mediterranean basin. The inter-RCM changes in the bias pattern with respect to the E-OBS are larger than the bias resulting from a change in the boundary conditions from ERA-40 to ECHAM5 20c3m. The maximum temperature response to increased green house gases, as projected by the A1B scenario, is consistent for both RCMs. Under that scenario, results indicate that the increments for extremes (e.g. 40-year return values) will be two or three times higher than those for the mean seasonal temperatures, particularly during Spring and Summer in Southern Europe.