Impact of the LMDZ atmospheric grid configuration on the climate and sensitivity of the IPSL-CM5A coupled model

The IPSL-CM5A climate model was used to perform a large number of control, historical and climate change simulations in the frame of CMIP5. The refined horizontal and vertical grid of the atmospheric component, LMDZ, constitutes a major difference compared to the previous IPSL-CM4 version used for CMIP3. From imposed-SST (Sea Surface Temperature) and coupled numerical experiments, we systematically analyze the impact of the horizontal and vertical grid resolution on the simulated climate. The refinement of the horizontal grid results in a systematic reduction of major biases in the mean tropospheric structures and SST. The mid-latitude jets, located too close to the equator with the coarsest grids, move poleward. This robust feature, is accompanied by a drying at mid-latitudes and a reduction of cold biases in mid-latitudes relative to the equator. The model was also extended to the stratosphere by increasing the number of layers on the vertical from 19 to 39 (15 in the stratosphere) and adding relevant parameterizations. The 39-layer version captures the dominant modes of the stratospheric variability and exhibits stratospheric sudden warmings. Changing either the vertical or horizontal resolution modifies the global energy balance in imposed-SST simulations by typically several W/m² which translates in the coupled atmosphere-ocean simulations into a different global-mean SST. The sensitivity is of about 1.2 K per 1 W/m² when varying the horizontal grid. A re-tuning of model parameters was thus required to restore this energy balance in the imposed-SST simulations and reduce the biases in the simulated mean surface temperature and, to some extent, latitudinal SST variations in the coupled experiments for the modern climate. The tuning hardly compensates, however, for robust biases of the coupled model. Despite the wide range of grid configurations explored and their significant impact on the present-day climate, the climate sensitivity remains essentially unchanged.     

Effect of land model ensemble versus coupled model ensemble on the simulation of precipitation climatology and variability

Through a series of model simulations with an atmospheric general circulation model coupled to three different land surface models, this study investigates the impacts of land model ensembles and coupled model ensemble on precipitation simulation. It is found that coupling an ensemble of land models to an atmospheric model has a very minor impact on the improvement of precipitation climatology and variability, but a simple ensemble average of the precipitation from three individually coupled land-atmosphere models produces better results, especially for precipitation variability. The generally weak impact of land processes on precipitation should be the main reason that the land model ensembles do not improve precipitation simulation. However, if there are big biases in the land surface model or land surface data set, correcting them could improve the simulated climate, especially for well-constrained regional climate simulations.     

Surface energy balance complexity in GCM land surface models. Part II: coupled simulations

 Global coupled simulations with the Bureau of Meteorology Research Centre climate model and the CHAmeleon Surface Model (CHASM) are used to examine how four general extensions to the representation of the basic land surface energy balance affect simulated land-atmosphere interface variables: evaporation, precipitation, skin temperature and air temperature. The impacts of including separate surface energy balance calculations for: vegetated and non-vegetated portions of the land surface; an explicit parametrisation of canopy resistance; explicit bare ground evaporation; and explicit canopy interception are isolated and quantified. The hypothesis that these aspects of surface energy balance parametrisation do not contain substantial information at the monthly time scale (and are therefore not important to consider in a land surface model) is shown to be false. Considerable sensitivity to each of the four general surface energy balance extensions is identified in average pointwise monthly changes for important land-atmosphere interface variables. Average pointwise changes in monthly precipitation and land evaporation are equal to about 40 and 31–37% of the global-average precipitation and land evaporation respectively. Average pointwise changes for land surface skin temperature and lowest model layer air temperature are about 2 and 0.9 K respectively. The average pointwise change and average pointwise biases are statistically significant at 95% in all cases. Substantial changes to zonally average variables are also identified. We demonstrate how the globally averaged surface resistance parameter can vary from 150 to 25 s/m depending on which aspects of the surface energy balance are treated implicitly. We also show that if interception is treated implicitly, the effective surface resistance must vary geographically in order to capture the behaviour of a model which treats this process explicitly. The implication of these results for the design of land surface models is discussed. Received: 8 July 1999 / Accepted: 1 September 2000     

LMDZ5B: the atmospheric component of the IPSL climate model with revisited parameterizations for clouds and convection

Based on a decade of research on cloud processes, a new version of the LMDZ atmospheric general circulation model has been developed that corresponds to a complete recasting of the parameterization of turbulence, convection and clouds. This LMDZ5B version includes a mass-flux representation of the thermal plumes or rolls of the convective boundary layer, coupled to a bi-Gaussian statistical cloud scheme, as well as a parameterization of the cold pools generated below cumulonimbus by re-evaporation of convective precipitation. The triggering and closure of deep convection are now controlled by lifting processes in the sub-cloud layer. An available lifting energy and lifting power are provided both by the thermal plumes and by the spread of cold pools. The individual parameterizations were carefully validated against the results of explicit high resolution simulations. Here we present the work done to go from those new concepts and developments to a full 3D atmospheric model, used in particular for climate change projections with the IPSL-CM5B coupled model. Based on a series of sensitivity experiments, we document the differences with the previous LMDZ5A version distinguishing the role of parameterization changes from that of model tuning. Improvements found previously in single-column simulations of case studies are confirmed in the 3D model: (1) the convective boundary layer and cumulus clouds are better represented and (2) the diurnal cycle of convective rainfall over continents is delayed by several hours, solving a longstanding problem in climate modeling. The variability of tropical rainfall is also larger in LMDZ5B at intraseasonal time-scales. Significant biases of the LMDZ5A model however remain, or are even sometimes amplified. The paper emphasizes the importance of parameterization improvements and model tuning in the frame of climate change studies as well as the new paradigm that represents the improvement of 3D climate models under the control of single-column case studies simulations.     

Long-term ice sheet–climate interactions under anthropogenic greenhouse forcing simulated with a complex Earth System Model

Several multi-century and multi-millennia simulations have been performed with a complex Earth System Model (ESM) for different anthropogenic climate change scenarios in order to study the long-term evolution of sea level and the impact of ice sheet changes on the climate system. The core of the ESM is a coupled coarse-resolution Atmosphere–Ocean General Circulation Model (AOGCM). Ocean biogeochemistry, land vegetation and ice sheets are included as components of the ESM. The Greenland Ice Sheet (GrIS) decays in all simulations, while the Antarctic ice sheet contributes negatively to sea level rise, due to enhanced storage of water caused by larger snowfall rates. Freshwater flux increases from Greenland are one order of magnitude smaller than total freshwater flux increases into the North Atlantic basin (the sum of the contribution from changes in precipitation, evaporation, run-off and Greenland meltwater) and do not play an important role in changes in the strength of the North Atlantic Meridional Overturning Circulation (NAMOC). The regional climate change associated with weakening/collapse of the NAMOC drastically reduces the decay rate of the GrIS. The dynamical changes due to GrIS topography modification driven by mass balance changes act first as a negative feedback for the decay of the ice sheet, but accelerate the decay at a later stage. The increase of surface temperature due to reduced topographic heights causes a strong acceleration of the decay of the ice sheet in the long term. Other feedbacks between ice sheet and atmosphere are not important for the mass balance of the GrIS until it is reduced to 3/4 of the original size. From then, the reduction in the albedo of Greenland strongly accelerates the decay of the ice sheet.     

Comparison of the climate simulated by the CCM3 coupled to two different land-surface models

We present results from a coupled atmosphere-biosphere model CCM3/IBIS (the Community Climate Model coupled to the Integrated BIosphere Simulator), which is designed to study the dynamic interactions between climate and vegetation and the global carbon cycle. We analyze the climate simulated by CCM3/IBIS with fixed vegetation conditions and we compare it to the climate simulated by the standard CCM3, which includes the LSM (land surface model) land-surface package. Important differences between the two models include simple parametrizations of lakes, wetlands and crops in CCM3/LSM not taken into account in CCM3/IBIS. CCM3/IBIS and CCM3/LSM share common biases (compared to observations) in the temperature field in boreal winter and in the precipitation field annually, making the atmospheric model the most probable cause of those biases. The models differ in the temperature field and surface energy balance in the Sahara annually and in the mid-to high latitudes from spring through fall. CCM3/IBIS simulates global annual air temperatures that are on average 1.7 °C higher than CCM3/LSM and 0.5 °C higher than the observed climatology. Differences in albedo and/or snow parametrization explain most of the Sahara and high-latitude temperature disagreement. Our sensitivity study with CCM3/LSM shows that the presence of lakes and wetlands in CCM3/LSM can account for about half of the difference in temperature in summer over the lake and wetland regions of the mid-latitudes. A second sensitivity study shows that higher surface roughness length in CCM3/IBIS can also explain part of the difference in summer surface temperature in the mid-latitudes. Surface roughness length affects the surface temperature through a feedback mechanism linking surface wind speed, planetary boundary layer height, low level cloudiness and radiation     

An Updated Coupled Model for Land-Atmosphere Interaction. Part I： Simulations of Physical Processes

A new two-way land-atmosphere interaction model （R42_AVIM） is fulfilled by coupling the spectral atmospheric model （SAMIL_R42L9） developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences （LASG/IAP/CAS） with the land surface model, Atmosphere-Vegetation-Interaction-Model （AVIM）. In this coupled model, physical and biological components of AVIM are both included. Climate base state and land surface physical fluxes simulated by R42_AVIM are analyzed and compared with the results of R42_SSIB [which is coupled by SAMIL_R42L9 and Simplified Simple Biosphere （SSIB） models]. The results show the performance of the new model is closer to the observations. It can basically guarantee that the land surface energy budget is balanced, and can simulate June-July-August （JJA） and December-January- February （DJF） land surface air temperature, sensible heat flux, latent heat flux, precipitation, sea level pressure and other variables reasonably well. Compared with R42_SSIB, there are obvious improvements in the JJA simulations of surface air temperature and surface fluxes. Thus, this land-atmosphere coupled model will offer a good experiment platform for land-atmosphere interaction research.     

A Regional Air--Sea Coupled Model and Its Application over East Asia in the Summer of 2000

A regional air-sea coupled model,comprising the Regional Integrated Environment Model System (RIEMS)and the Princeton Ocean Model(POM)was developed to simulate summer climate features over East Asia in 2000.The sensitivity of the model's behavior to the coupling time interval(CTI),the causes of the sea surface temperature(SST)biases,and the role of air-sea interaction in the simulation of precipitation over China are investigated.Results show that the coupled model can basically produce the spatial pattern of SST,precipitation,and surface air temperature(SAT)with five different CTIs respectively.Also,using a CTI of 3,6 or 12 hours tended to produce more successful simulations than if using 1 and 24 hours.Further analysis indicates that both a higher and lower coupling frequency result in larger model biases in air-sea heat flux exchanges,which might be responsible for the sensitivity of the coupled model's behavior to the CTI. Sensitivity experiments indicate that SST biases between the coupled and uncoupled POM occurring over the China coastal waters were due to the mismatch of the surface heat fluxes produced by the RIEMS with those required by the POM.In the coupled run,the air-sea feedbacks reduced the biases in surface heat fluxes,compared with the uncoupled RIEMS,consequently resulted in changes in thermal contrast over land and sea and led to a precipitation increase over South China and a decrease over North China.These results agree well observations in the summer of 2000.     

An Updated Coupled Model for Land-Atmosphere Interaction. Part Ⅰ: Simulations of Physical Processes

A new two-way land-atmosphere interaction model (R42_AVIM) is fulfilled by coupling the spectral at- mospheric model (SAMIL_R42L9) developed at the State Key Laboratory of Numerical Modeling for Atmo- spheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sci- ences (LASG/IAP/CAS) with the land surface model, Atmosphere-Vegetation-Interaction-Model (AVIM). In this coupled model, physical and biological components of AVIM are both included. Climate base state and land surface physical fluxes simulated by R42_AVIM are analyzed and compared with the results of R42_SSIB [which is coupled by SAMIL_R42L9 and Simplified Simple Biosphere (SSIB) models]. The results show the performance of the new model is closer to the observations. It can basically guarantee that the land surface energy budget is balanced, and can simulate June-July-August (JJA) and December-January- February (DJF) land surface air temperature, sensible heat flux, latent heat flux, precipitation, sea level pressure and other variables reasonably well. Compared with R42_SSIB, there are obvious improvements in the JJA simulations of surface air temperature and surface fluxes. Thus, this land-atmosphere coupled model will offer a good experiment platform for land-atmosphere interaction research.     

陆面过程模式CoLM和NCAR_CLM3.0对中国典型森林生态系统陆气相互作用的模拟 I. 不同模式模拟结果的初步分析

CoLM(Common Land Model)和NCAR_CLM3.0(NCAR Community Lanol Model 3.0)是目前国际上广为应用的两个发展比较完善的陆面过程模式.本研究利用中国陆地生态系统通量观测研究网络(Chi-naFLUX)在长白山温带混交林和千烟洲亚热带人工针叶林观测站点的长期连续强化观测资料,对这两个模式在上述地区的模拟性能进行了初步评估.同观测资料的对比表明,两个模式均能较好地模拟出观测站点地表能量和水分平衡的基本特征,其中,CoLM对潜热通量的模拟性能更好.以对潜热通量为期1年的日均值的模拟为例,COLM和CLM3.0在长白山观测站模拟值和观测值时间序列的相关系数分别为0.80和0.65,在千烟洲站分别为0.69和0.64,均通过了0.01的信度检验;两个模式对全年平均的模拟与观测日平均值的比值在长白山分别为1.21和0.86,在千烟洲分别为0.83和0.60.研究结果表明,这两个陆面过程模式可以作为研究这两种典型森林生态系统陆气交换的基本工具.同时,对模式模拟性能差异的深入分析将有助于进一步改进陆面模式的参数化过程,为相关研究奠定更坚实的基础.     

全球模式NCAR CESM和CAS ESM对亚洲东部夏季气候的模拟性能评估:气候平均态和降水日变化分析

主要评估了美国国家大气研究中心的NCAR CESM（Community Earth System Model,NCAR）和中国科学院的CAS ESM（Earth System Model,Chinese Academy of Sciences）两个地球系统模式对亚洲东部夏季气候态的模拟性能。使用NCAR CESM和CAS ESM各两种不同的水平分辨率,一共进行了4组长达19年（1998~2016年）的数值积分试验,并通过对2 m气温、降水强度和降水日变化等的分析,比较了这两个模式在亚洲东部的模拟性能。结果表明,CAS ESM和NCAR CESM均能模拟出夏季2 m气温和降水强度的大尺度分布特征,但整体上模拟得到的地表面气温偏暖、降水强度偏弱。对于降水日变化而言,观测的日降水峰值在陆地上主要发生在下午到傍晚时段,在海洋上则出现在午夜到凌晨时段。两组低分辨率试验模拟的陆地降水峰值出现过早,且无法模拟出四川盆地的夜间降水峰值和部分海洋地区凌晨或上午的降水峰值。提高分辨率对模式的模拟性能有显著的提升作用。高分辨率下,NCAR CESM和CAS ESM对陆地和海洋的降水日变化模拟性能都明显提高。对降水日变化的定量化分析表明,高分辨率CAS ESM模式对整个亚洲东部降水日变化的模拟最优。目前模式对海陆风的模拟还不太理想,未来要进一步提高模式模拟性能,需要重点完善与气温、降水过程相关的物理参数化方案。     

Dynamical downscaling of IPSL-CM5 CMIP5 historical simulations over the Mediterranean: benefits on the representation of regional surface winds and cyclogenesis

The Mediterranean region is identified as one of the two main hot-spots of climate change and also known to have the highest concentration of cyclones in the world. These atmospheric features contribute significantly to the regional climate but they are not reproduced by the Atmosphere–Ocean General Circulation Models (AOGCM), due to their coarse horizontal resolution, which have recently been run in the frame of the 5th Climate Model Intercomparison Project. This article investigates the benefit of dynamically downscaling the Institut Pierre Simon Laplace (IPSL) AOGCM (IPSL-CM5) historical simulation by the weather and research forecasting (WRF) for the representation of the Mediterranean surface winds and cyclonic activity. Indeed, when considering IPSL-CM5 atmospheric fields, the dramatic underestimation of the cyclonic activity in the most cyclogenetic region of the world jeopardizes our ability to investigate in-depth the Mediterranean regional climate and trend in the context of global change. The WRF model shows remarkable skill to reproduce regional cyclogenesis. Indeed, cyclones occurrence is quasi-absent in IPSL-CM5 data but when applying dynamical downscaling their spatial–temporal variability is very close to the re-analysis. This is a clear benefit of dynamical downscaling in regions of strong topographic forcing. This “steady” source of forcing allows the production of lee cyclogenesis and the development of strong cyclones, whatever the quality of the large-scale circulation provided at the WRF’s boundaries by IPSL-CM5. However, dynamical downscaling still presents disadvantages as for instance the fact that large-scale inaccurate features of the IPSL-CM5 regional circulation are replicated by WRF due to the boundary controlled (small domain) simulation. The advantages and disadvantages of dynamical downscaling are thoroughly discussed in this paper revealing its importance for climate research, especially in the context of future scenarios and wind impacts.     

陆面模式CLSM的设计及性能检验Ⅰ.模式设计

在现有陆面模式的基础上,详细考虑了地气系统中的积雪、土壤水热传输、植被及湍流边界层中的各种物理过程,发展了一个既能反映积雪变化、干旱/半干旱区地气交换过程,同时又能描述不同陆面状况物理过程的陆面模式(Comprehensive Land Surface Model,简称CLSM),改善了陆面模式对全球范围内不同下垫面条件下的陆面过程及地-气交换过程的模拟能力.     

Overview of the CMIP6 Historical Experiment Datasets with the Climate System Model CAS FGOALS-f3-L

The three-member historical simulations by the Chinese Academy of Sciences Flexible Global Ocean–Atmosphere–Land System model, version f3-L(CAS FGOALS-f3-L), which is contributing to phase 6 of the Coupled Model Intercomparison Project(CMIP6), are described in this study. The details of the CAS FGOALS-f3-L model, experiment settings and output datasets are briefly introduced. The datasets include monthly and daily outputs from the atmospheric, oceanic, land and sea-ice component models of CAS FGOALS-f3-L, and all these data have been published online in the Earth System Grid Federation(ESGF, https://esgf-node.llnl.gov/projects/cmip6/). The three ensembles are initialized from the 600th, 650th and 700th model year of the preindustrial experiment(piControl) and forced by the same historical forcing provided by CMIP6 from 1850 to 2014. The performance of the coupled model is validated in comparison with some recent observed atmospheric and oceanic datasets. It is shown that CAS FGOALS-f3-L is able to reproduce the main features of the modern climate, including the climatology of air surface temperature and precipitation,the long-term changes in global mean surface air temperature, ocean heat content and sea surface steric height, and the horizontal and vertical distribution of temperature in the ocean and atmosphere. Meanwhile, like other state-of-the-art coupled GCMs, there are still some obvious biases in the historical simulations, which are also illustrated. This paper can help users to better understand the advantages and biases of the model and the datasets.     

Mid-Century Ensemble Regional Climate Change Scenarios for the Western United States

To study the impacts of climate change on water resources in the western U.S., global climate simulations were produced using the National Center for Atmospheric Research/Department of Energy (NCAR/DOE) Parallel Climate Model (PCM). The Penn State/NCAR Mesoscale Model (MM5) was used to downscale the PCM control (20 years) and three future(2040–2060) climate simulations to yield ensemble regional climate simulations at 40 km spatial resolution for the western U.S. This paper describes the regional simulations and focuses on the hydroclimate conditions in the Columbia River Basin (CRB) and Sacramento-San Joaquin River (SSJ) Basin. Results based on global and regional simulations show that by mid-century, the average regional warming of 1 to 2.5 °C strongly affects snowpack in the western U.S. Along coastal mountains, reduction in annual snowpack was about70% as indicated by the regional simulations. Besides changes in mean temperature, precipitation, and snowpack, cold season extreme daily precipitation increased by 5 to 15 mm/day (15–20%) along theCascades and the Sierra. The warming resulted in increased rainfall at the expense of reduced snowfall, and reduced snow accumulation (or earlier snowmelt) during the cold season. In the CRB, these changes were accompanied by more frequent rain-on-snow events. Overall, they induced higher likelihood of wintertime flooding and reduced runoff and soil moisture in the summer. Changes in surface water and energy budgets in the CRB and SSJ basin were affected mainly by changes in surface temperature, which were statistically significant at the 0.95 confidence level. Changes in precipitation, while spatially incoherent, were not statistically significant except for the drying trend during summer. Because snow and runoff are highly sensitive tospatial distributions of temperature and precipitation, this study shows that (1) downscaling provides more realistic estimates of hydrologic impacts in mountainous regions such as the western U.S., and (2) despite relatively small changes in temperature and precipitation, changes in snowpack and runoff can be much larger on monthly to seasonal time scales because the effects of temperature and precipitation are integrated over time and space through various surface hydrological and land-atmosphere feedback processes. Although the results reported in this study were derived from an ensemble of regional climate simulations driven by a global climate model that displays low climate sensitivity compared with most other models, climate change was found to significantly affect water resources in the western U.S. by the mid twenty-first century.     

Global validation of the ISBA sub-grid hydrology

Over recent years, many numerical studies have suggested that the land surface hydrology contributes to atmospheric variability and predictability on a wide range of scales. Conversely, land surface models (LSMs) have been also used to study the hydrological impacts of seasonal climate anomalies and of global warming. Validating these models at the global scale is therefore a crucial task, which requires off-line simulations driven by realistic atmospheric fluxes to avoid the systematic biases commonly found in the atmospheric models. The present study is aimed at validating a new land surface hydrology within the ISBA LSM. Global simulations are conducted at a 1° by 1° horizontal resolution using 3-hourly atmospheric forcings provided by the Global Soil Wetness Project. Compared to the original scheme, the new hydrology includes a comprehensive and consistent set of sub-grid parametrizations in order to account for spatial heterogeneities of topography, vegetation, and precipitation within each grid cell. The simulated runoff is converted into river discharge using the total runoff integrating pathways (TRIP) river routing model (RRM), and compared with available monthly observations at 80 gauging stations distributed over the world’s largest river basins. The simulated discharges are also compared with parallel global simulations from five alternative LSMs. Globally, the new sub-grid hydrology performs better than the original ISBA scheme. Nevertheless, the improvement is not so clear in the high-latitude river basins (i.e. Ob, MacKenzie), which can be explained by a too late snow melt in the ISBA model. Over specific basins (i.e. Parana, Niger), the quality of the simulated discharge is also limited by the TRIP RRM, which does not account for the occurrence of seasonal floodplains and for their significant impact on the basin-scale water budget.     

Boreal forest and tundra ecosystems as components of the climate system

The effects of terrestrial ecosystems on the climate system have received most attention in the tropics, where extensive deforestation and burning has altered atmospheric chemistry and land surface climatology. In this paper we examine the biophysical and biogeochemical effects of boreal forest and tundra ecosystems on atmospheric processes. Boreal forests and tundra have an important role in the global budgets of atmospheric CO₂ and CH₄. However, these biogeochemical interactions are climatically important only at long temporal scales, when terrestrial vegetation undergoes large geographic redistribution in response to climate change. In contrast, by masking the high albedo of snow and through the partitioning of net radiation into sensible and latent heat, boreal forests have a significant impact on the seasonal and annual climatology of much of the Northern Hemisphere. Experiments with the LSX land surface model and the GENESIS climate model show that the boreal forest decreases land surface albedo in the winter, warms surface air temperatures at all times of the year, and increases latent heat flux and atmospheric moisture at all times of the year compared to simulations in which the boreal forest is replaced with bare ground or tundra. These effects are greatest in arctic and sub-arctic regions, but extend to the tropics. This paper shows that land-atmosphere interactions are especially important in arctic and sub-arctic regions, resulting in a coupled system in which the geographic distribution of vegetation affects climate and vice versa. This coupling is most important over long time periods, when changes in the abundance and distribution of boreal forest and tundra ecosystems in response to climatic change influence climate through their carbon storage, albedo, and hydrologic feedbacks.